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Abstract 

Correlations in the nuclear wave-function beyond the mean-field or Hartree-Fock ap- 
proximation are very important to describe basic properties of nuclear structure. Various 
approaches to account for such correlations are described and compared to each other. 
This includes the hole-line expansion, the coupled cluster or "exponential S" approach, 
the self-consistent evaluation of Greens functions, variational approaches using correlated 
basis functions and recent developments employing quantum Monte-Carlo techniques. 
Details of these correlations are explored and their sensitivity to the underlying nucleon- 
nucleon interaction. Special attention is paid to the attempts to investigate these correla- 
tions in exclusive nucleon knock-out experiments induced by electron scattering. Another 
important issue of nuclear structure physics is the role of relativistic effects as contained 
in phenomenological mean field models. The sensitivity of various nuclear structure ob- 
servables on these relativistic features are investigated. The report includes the discussion 
of nuclear matter as well as finite nuclei. 

1 Introduction 



One of the central challenges of theoretical nuclear physics is the attempt to describe the basic 
properties of nuclear systems in terms of a realistic nucleon-nucleon (NN) interaction. Such 
an attempt typically contains two major steps. In the first step one has to consider a specific 
model for the NN interaction. This could be a model which is inspired by the quantum- 
chromo-dynamics|]l]], a meson-exchange or One-Boson- Exchange model0, ^ or a purely phe- 
nomenological ansatz in terms of two-body spin-isospin operators multiplied by local potential 
functions [|, |[. Such models are considered as a realistic description of the NN interaction, if 
the adjustment of parameters within the model yields a good fit to the NN scattering data at 
energies below the threshold for pion production as well as energy and other observables of the 
deuteron. 

After the definition of the nuclear hamiltonian, the second step implies the solution of 
the many-body problem of A nucleons interacting in terms of such a realistic two-body NN 
interaction. The simplest approach to this many-body problem of interacting fermions one 
could think of would be the mean field or Hartree-Fock approximation. This procedure yields 
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very good results for the bulk properties of nuclei, binding energies and radii, if one employs 
simple plienomenological NN forces like e.g. the Skyrme forces, which are adjusted to describe 
such nuclear structure data[0. However, employing realistic NN interactions the Hartree-Fock 
approximation fails very badly: it leads to unbound nuclei [0]. 

This failure of the Hartree-Fock approximation is a consequence of the strong short-range 
components of a realistic NN interaction, which are necessary to describe the NN data. The 
Hartree-Fock wavefunction describes the nucleus as a system of nucleons moving independent 
from each other in a mean field derived from the average interaction with all other nucleons. 
This implies that the wavefunction contains large amplitudes of configurations, in which two 
nucleons are so close to each other, that they are exposed to the very repulsive components of 
the NN interaction at short distances. The hard-core potentials [§], which were very popular 
in the sixties, describe these components in terms of an infinitely repulsive core for relative 
distances smaller than the radius of this hard core of about 0.4 fm. In this case a Hartree- 
Fock calculation would even yield an infinite repulsive energy. Modern models for the NN 
interaction, in particular the meson-exchange models which lead to non-local NN interactions, 
contain softer cores. Nevertheless, a careful treatment of two-body short range correlations 
beyond the mean field approximation is indispensable to describe the structure of nuclei in 
terms of realistic NN interactions. The same is true for the correlations which are induced 
by the strong tensor components in the NN interaction, which mainly originate from the pion 
exchange contributions. 

Various different tools have been developed to account for such correlation effects. Below we 



will describe some of these methods including the Brueckner hole-line expansionP, 10, 11], the 



coupled cluster or "exponential S" approachO, O], the self-consistent evaluation of Greens 



functions[Q, variational approaches using correlated basis functions [p!5| , |T6|, and recent 



developments employing quantum Monte-Carlo techniques ||18|, |IS 



Using such methods one obtains correlated many-body wave functions, which are rather 
sensitive to the NN interaction under consideration. Therefore the question arises which kind 
of experimental observables might be considered to investigate details of these correlations. 
The hope is that different predictions derived from the various model of the NN interaction 
will allow to distinguish between the various model for the NN interaction at short distances. 
Therefore, finally such nuclear structure studies will help to explore the details of the strong 
baryon baryon interaction at short distances. 

What is the effect of correlations on the nuclear wave function? In order to discuss this 
question let us consider for the moment a system of infinite nuclear matter. The mean field or 
Hartree-Fock wave function of nuclear matter corresponds to the Slater determinant of plane 
wave states, in which all single-particle states with momenta below the Fermi momentum kp are 
occupied. Correlations on top of this Hartree-Fock state will yield a depletion of single-particle 
states with momenta k below kp and a non-vanishing occupation of states with high momentum. 
From this consideration one may think that the study of high-momentum components in the 
single-particle wave functions of nuclear states might be an ideal tool to explore correlation 
effects in the nuclear wave function. Therefore it has been suggested to study these high- 
momentum components by means of exclusive single nucleon knock-out experiments like (e, e'A^) 
with the {A — 1) nucleus remaining in the ground state or other well defined statep^, pl[| . 
However, a detailed analysis of the spectral function for such knock-out experiments shows that 
the expected high-momentum components in the nuclear wavefunction could only be observed 
at high missing energies, i.e. with an excitation energy of the residual target nucleus well above 
the threshold for the emission of a secondary nucleon [E^, |2^, p4|. Therefore one expects exclusive 



two-nucleon knock-out experiments like (7, A^A^) or (e, e'NN) ||25|, |26[| to be more sensitive to the 



effects of correlations. Details for the analysis of such nucleon knock-out experiments induced 
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by inelastic electron scattering will be given below. 

Up to this point we merely discussed the nuclear system within the framework of the non- 
relativistic quantum many-body theory. At first sight this seems to be well justified since the 
binding energies of nuclei with single-particle potentials of -50 MeV or so, are very small on 
the scale of the nucleon rest mass. It has already been discussed, however, that the short-range 
components of the NN interaction are very strong with large repulsive components, which 
are compensated by strong attractive ones. Within the meson-exchange model for the NN 
interaction the repulsive components are generated by the exchange of the cu, a vector meson, 
while the attractive parts are described in terms of the exchange of a scalar particle, the o 
meson 0. This o meson does not represent a meson in the usual sense but its exchange is used 
to describe the correlated two-pion exchange without and with intermediate excitation of the 



interacting nucleonsp7|, p8| . Calculating the nucleon self energies from such a meson exchange 
model within a Hartree approximation, one finds that the io exchange yields a component So, 
which transforms under a Lorentz transformation like the time-like component of a vector, while 
the scalar meson exchange yields a contribution S^, which transforms like a scalar. Inserting 
this self energy into the Dirac equation for a nucleon in the medium of nuclear matter leads to 
binding effects which are as small as the -50 MeV discussed above. This small binding effect, 
however, results from a strong cancellation between the repulsive Sq and the attractive 
component. The attractive scalar component leads to Dirac spinors for the nucleons in the 
nuclear medium, which contain a small component significantly enhanced as compared to the 
Dirac spinor of a free nucleon. This effect is often described in terms of an effective Dirac mass 
m* for the nucleon, which can be of the order of 600 MeV in nuclear matter around saturation 
density. 

This modification of the Dirac spinors in the nuclear medium requires a self-consistent eval- 
uation of the matrix elements of the meson exchange interaction. Within the phenomenological 
0" — or Walecka model of nuclear physics p9|, pO|] it is this self-consistency, which provides 
the saturation of nuclear matter. This relativistic effect, however, is also observed in Dirac- 



Brueckner-Hartree-Fock studies, which are based on realistic NN interactions fHT 



Including these relativistic features improves the predictions for the saturation point of nuclear 
matter as well as the bulk properties of finite nuclei significantly The question 



is: Are there other observables, which are sensitive to the enhancement of the small component 
of the nucleon Dirac spinors in the nuclear medium? 

The calculation scheme discussed so far, determine the interaction of two nucleons in the 
vacuum in a first step and then solve the many-body problem of nucleons interacting by such 
realistic potentials in a second step, is of course based on the picture that nucleons are ele- 
mentary particles with properties, which are not affected by the presence of other nucleons in 
the nuclear medium. One knows, of course, that this is a rather simplified picture: nucleons 
are built out of quarks and their properties might very well be influenced by the surrounding 
medium. A cartoon of this feature is displayed in Fig. ^ So the question is, how important are 
the sub-nucleonic degrees of freedom, must we expect a change of the nucleon properties in the 
nuclear medium? At low energies it might be sufficient to account for the internal structure 
of the nucleons by considering the possibility that the nucleon may get excited to the A(3,3) 
excitation at 1232 MeV. It has been demonstrated |3^ that the process of two interacting nu- 
cleons, which polarize each other to such isobar excitations are an important ingredient to the 
medium range attraction of the NN interaction. Attempts have been made to account for such 
isobar excitations also in microscopic studies of nuclear structure EDI, ET . 



Furthermore, not only the basic properties of the nucleons, their mass or electromagnetic 
form-factor might be affected by the nuclear medium, also the masses and thereby the propa- 
gation of the mesons might be modified as well[H2], H3. This will have consequences for the 
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Figure 1: Cartoon of a nucleus, displaying the size of the nucleons as compared to the 
typical distance to nearest neighbors. Also indicated are the internal structure of nucleons 
and mesons. 




Figure 2: Two nucleons interacting in a nuclear medium, see discussion in the text. 



meson-exchange model of the NN interaction. For example, a lowering of the meson masses in 
the nuclear medium as compared to the vacuum would lead to a larger range of the correspond- 
ing meson exchange interaction terms, which should lead to different results in the nuclear 
structure calculation. 

In studying such features of sub-nucleonic degrees of freedom, however, one should be careful 
to avoid the mixing of different points of view on the same process. As an example, let us 
consider a process like the one indicated by the diagram in Fig. 0: Two nucleons, represented 
by the two upward going lines, interact with each other by the exchange of e.g. a pion, which is 
represented by the dashed line. Propagating in a nuclear medium between these two nucleons, 
the pion may interact with a third nucleon leading to an intermediate particle-hole excitation of 
the nuclear system. Such processes might be considered as a modification of the pion propagator 
in the nuclear medium, which could be characterized by a modification of the effective pion 
mass inside a nucleus. The process displayed in Fig. |, however, will also be included in any 
many-body calculation based on two-nucleon interactions, which accounts for three-nucleon 
correlations or includes effects of ring diagrams. Therefore effects which one may identify with 
a modification of meson propagators are accounted for in an approach, which does not even 
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mention any mesonic degrees of freedom. 

The pion re-scattering process displayed in Fig. ^ may also lead to an excitation of the 
intermediate nucleon to e.g. the A resonance. Again such terms might be considered as a 
modification of the meson propagator. Such terms would also be accounted for in a many-body 
theory based on baryon - baryon interactions, like e.g. the Argonne V28|^5[, which includes 
the A degrees of freedom explicitly. Within the framework of a many-body theory which does 
not account for such A excitations, such a process might be taken into account by means of a 
three-nucleon interaction. This example has been given to demonstrate the model dependence 
of three-body forces or statements about the relevance of mesonic degrees of freedom. The 
approaches which we will discuss throughout this review, are all based on realistic NN inter- 
actions and try to solve the many-body problem as good as possible. It is the aim to study 
how far such an approach, which ignores all sub-nucleonic degrees of freedom, can predict the 
properties of nuclear systems at zero temperature and without external pressure. Remaining 
discrepancies between theoretical prediction and empirical data can then bet attributed to the 
necessity to account for sub-nucleonic degrees of freedom, explicitly. 

Such many-body calculations can then also be used to obtain predictions for nuclear matter 
under extreme conditions like in neutron stars, supernova explosions |[46|| or in central heavy 
ion collisions. They yield a prediction for the equation of state of nuclear matter at high 
densities and/or high temperatures which is free of any parameters. The success or failure of 
such theoretical calculations in predicting the properties of normal nuclear matter from the 
NN interactions of two nucleons in the vacuum can be used to judge the reliability of these 
predictions for nuclear matter at extreme conditions. 

Nuclear systems provide an intriguing and challenging object for the development of quan- 
tum many-body theories. The underlying NN interaction is non-trivial: It contains a rich 
operator structure with strong tensor components and is in general non-local. The interacting 
particles must be considered as quasiparticle and the inclusion of internal degrees of freedom 
might get important. An extension of the many-body theory to account for relativistic effects 
might be necessary. Empirical data are available and should be reproduced for finite systems 
with particle numbers ranging from A = 2 to infinite nuclear matter. For other systems of 
condensed matter or Fermi liquids, like electron gas or liquid Helium 3, reliable data on finite 
samples are difficult to obtain but of high interest. Therefore nuclei are an ideal testing ground 
for many-body theory for finite systems in particular. The experience collected here should be 
rather useful in the study of other systems. 

After this introduction, we will present an outline of some many-body approaches which are 
used in nuclear physics. This includes the Brueckner - Bethe hole-line expansion, the coupled 
cluster or "exponential S" approach, the self-consistent evaluation of Greens functions, varia- 
tional approaches using correlated Jastrow-type basis functions and variational and quantum 
Monte Carlo techniques. A short review on the status of realistic NN interactions and on the 
differences between those models will be given in the first part of section 3. This section 3 
also contains the discussion of results for bulk properties of infinite nuclear matter and finite 
nuclei. As an example on experimental efforts to explore the effects of correlations in detail, we 
will report on the analysis of inelastic electron scattering in section 4. The section 5 describes 
attempts to extend the many-body theory to account for relativistic features as well as the 
sensitivity of some observables to these relativistic effects. 



2 Many-Body Approaches 
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Figure 3: Schematic picture of a NN interaction with hard core and its effect on the correlated 
NN wave function \l/(r). 



2.1 Hole - Line Expansion 

As it has been discussed aheady above one problem of nuclear structure calculations based on 
realistic NN interactions is to deal with the strong short-range components contained in all 
such interactions. This problem is evident in particular when so-called hard-core potentials 
are employed, which are infinite for relative distances smaller than the radius of the hard core 
Tc- The matrix elements of such a potential V evaluated for an uncorrelated two-body wave 
function $(r) diverges since $(r) is different from zero also for relative distances r smaller 
than the hard-core radius Vc (see the schematic picture in Fig. |^. A way out of this problem 
is to account for the two-body correlations induced by the NN interaction in the correlated 
wave function \E'(r) or by defining an effective operator, which acting on the uncorrelated wave 
function $(r) yields the same result as the bare interaction V acting on \E'(r). This concept is 
well known for example in dealing with the scattering matrix T, which is defined by 

< $|T|$ >=< > . (1) 

As it is indicated in the schematic Fig. ^ the correlations tend to enhance the amplitude of 
the correlated wave function relative to the uncorrelated one at distances r for which the 
interaction is attractive. A reduction of the amplitude is to be expected for small distances for 
which V{r) is repulsive. From this discussion we see that the correlation effects tend to make 
the matrix elements of T more attractive than those of the bare potential V. For two nucleons 
in the vacuum the T matrix can be determined by solving a Lippmann-Schwinger equation 

T|<l>> = vl\<^>+ V|^>1 

I u — Ho + ie J 

= \v + V -t1|$>. (2) 

In these equations the starting energy u stands for the energy of the two interacting nucleons 
and Hq represents the operator for the nucleons in the intermediate state without residual 
interaction, i.e. the kinetic energy. 

This concept of an effective operator accounting for correlation effects in the wave function is 
one way to present the Brueckner-Bethe-Goldstone (BBG) approach to the many-body problem. 
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To show the features of the BBG approach, we consider the nuclear hamiltonian, a sum of the 
kinetic energy tkin and the two-body interaction V, and introduce an appropriate single-particle 
potential U 

H = Ho + Hi, Ho = hra + U, Hi = V-U. (3) 

The eigenstates of Hq for A particles, i.e. Slater determinants build from the corresponding 
single-particle wavefunctions, provide a basis of the Hilbert-space for the A nucleon problem 

HoH^,...,A) = E^<i>,{l,...,A), (4) 

with the sum of the single particle energies. This basis <I>j is used to split the Hilbert-space 
into a model space and the rest. If one is interested in the ground-state properties of closed 
shell nuclei like ^^O this model space could be defined by a single Slater determinant. If one is 
interested in the states of open shell nuclei at low energies, this model space could consist out of 
those functions $i which represent a typical set of basis states for a shell-model configuration 



mixing calculation^^. The concept of a model space, however, can also be considered for 
the system of infinite nuclear matter. One may choose the model space to contain all Slater 
determinants, which can be constructed by considering occupied states with momenta below 
a model space momentum kM- If this momentum kM is equal to the Fermi momentum of a 
free Fermi gas kp of the density considered, the model space reduces just to this model wave 



function P8||. With this choice for a model space one can define projection operators V and Q, 
which project onto the model space and on the rest of the A-nucleon Hilbert space, respectively. 
It is the aim to derive an effective hamiltonian He/f, which is defined within the model space, 
with eigenvalues Ei identical to the eigenvalues of the original hamiltonian H and eigenvectors, 
which are just the projection of the exact eigenvectors of H on the model space. This implies 
that 

VH,jfV^, = E,V^,. (5) 

With this separation of the Hilbert-space into a model space and a rest also the treatment of 
correlations is separated: There are correlations which can be represented by degrees of freedom 
within this model space. They are treated explicitly and accounted for by the diagonalization 
of the effective hamiltonian Heff (see Eq. (^)) within the model-space. This will typically 
be correlations which are described in terms of single-particle excitations around the Fermi 
energy with low excitation energy, small momentum transfers and therefore of longer range. 
These long range correlations are quite sensitive to the shell structure of finite systems and 
will therefore be different for nuclear matter and specific finite nuclei. The correlations, which 
cannot be represented by the mixing of configurations within the model space are treated by 
means of the effective operator. Such correlations include excitations up to states high above 
the Fermi surface. This implies large excitation energies and high momentum transfers, which 
means correlations of short range. 

Several expansions have been formulated for this effective hamiltonian. There are the 
energy-dependent expansions of Feshbach and of Bloch and Horowitz . They correspond 
to the Brillouin - Wigner perturbation expansion and therefore yield an effective hamiltonian 
which depends on the exact energy to be calculated. To get rid of this energy-dependence in 
the case of multi-dimensional degenerate model-spaces one has to consider the so-called folded 



diagram expansion, which has been formulated e.g. by Brandow |51] and Kuo, Lee and Ratcliff 



52| . Using this folded- diagram formulation, the effective hamiltonian could be written as 

Q 

—Upff + toiaea aiagrams } 

linked 

with E^ referring to the eigenvalue of Hq for the state in the model space considered. 



Heft = Ho + Hi + \Hi ^ Heff + folded diagrams) , (6) 

I E^ - Ha J linked 
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Figure 4: Example of an unlinked diagram, which should be ignored in the expansion of 
H,ff in Eq. (g) 



Note that in this expansion one only needs to consider the contribution of the so-called linked 
diagrams. This implies that contributions like the one characterized by the diagram in Fig. ^ 
can be ignored. Such unlinked contributions contain products of matrix elements of Hi (the 
two-body parts are represented by the dashed lines in Fig. which are completely unlinked in 
the sense that the summation over single-particle quantum numbers of the states attached to 
the matrix elements of Hi are completely disconnected. Translated into the representation of 
diagrams this means that the corresponding diagram can be separated into two pieces without 
cutting a solid or dashed line, representing the propagation of a nucleon and an interaction, 
respectively. This restriction to linked diagram contribution is of particular importance for 
systems with many particles, like e.g. infinite nuclear matter. Unlinked contributions, which 
represent a combination of processes which are independent from each other, diverge in the 
limit of particle number to infinity. Details on the definition and on techniques to evaluate the 
contribution of folded diagrams can be found in |5l|, |5^, Folded diagrams do not 



occur if the model space is of dimension one, an assumption, which is most commonly made in 
calculating ground-state properties of closed shell nuclei or infinite nuclear matter. 

From the expansion in Eq. (^) it is obvious that H^ff will not only contain one-body and 
two-body operators but, in general, also operators involving three and more nucleons. If, 
however, we restrict the discussion to the evaluation of ground-state properties and consider 
only the effective two-body operators in Heff, the expansion is reduced to the summation of 
ladder diagrams and one obtains the Bethe-Goldstone equation 

Gico) = V + V-^G, (7) 

UJ - hi2 

which defines an effective NN interactions of two nucleons. Note that the projection operator 
on many-body states Q outside the model space in (|]) has been replaced by the Pauli operator 
Q, which is a two-body operator defined by 

Q\^j^j I > if i and j are single-particle unoccupied in $ ^g-j 
I else 

with $ referring to the Slater-determinant defining the model-space. The energy denominator 
in (0) corresponds to the excitation energy of the intermediate two-particle two-hole state, 
i.e. it is the sum of the single-particle energies for the interacting nucleons in states below 
the Fermi energy, expressed by the starting energy uj minus the energy of the two nucleons in 
the intermediate states above the Fermi energy, denoted by the operator hi2. (See also the 
graphical representation in Fig. ^) 
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The Bethe-Goldstone equation (|^) is quite similar to the Lippmann-Schwinger Eq. for the 
scattering matrix T. The only differences are the Pauli operator Q and the energy denominator 
which in is defined in terms of single-particle energies of the many-body system instead of 
the kinetic energies in Eq. (|^). Therefore the solution of the Bethe-Goldstone equation, the 
G-matrix, corresponds to an effective interaction between two nucleons, which accounts for 
correlation effects in the nuclear medium. Similar to Eq. (0) one can define for each product 
wave function of two uncorrelated single-particle wave functions > a correlated wave 

function lipap > by 

G\(j)at3 > = V\iJaf3 > 

= v{l + —^G]\<l>.f,>, (9) 
which implies that the correlated wave function \ipai3 > can be identified with 

>= |0a/3 > + ^7— > (10) 

to - hi2 

The second term on the right hand side of this equation, the difference between the correlated 
and uncorrelated wave function, is called the defect function. If the uncorrelated state refers to 
two single-particle states a and (3 below the Fermi surface, the difference between the starting 
energy u = ea + {ea denoting the single-particle energies of Hq) and the eigenvalue of hi2 is 
negative for all intermediate states, which are restricted to two-particle states above the Fermi 
surface. This means that the summation or integration over intermediate particle states does 
not meet any pole in the propagators of (|^) or (|T0|). This implies that the matrix elements 
of G are real. There exist not phase shifts between the uncorrelated and the correlated wave 
function, the defect function vanishes for large relative distances. This vanishing of the defect 
function at large r is called the healing property, the correlated wave functions "heals" to the 
uncorrelated one at large r. 

One of the main points of the Brueckner-Bethe-Goldstone approach to the many-body 
system is to evaluate the contributions to the effective hamiltonian -ffe// in Eq- ® not in a 
perturbation expansion in which the contributions are ordered with respect to the numbers of 
bare interaction V terms, but in terms of the G-matrix. 

Up to this point we have not specified yet, how to choose the unperturbed hamiltonian Hq or 
the single-particle potential U in Eq. (H). The aim is of course to choose Hq such that the Slater 
determinant $ defining the model space is close to the exact wavefunction for the ground-state. 
Therefore it seems quite natural to define the single-particle potential U in analogy to the 
Hartree-Fock definition with the bare interaction V replaced by the corresponding G-matrix. 
To be more precise, the Brueckner-Hartree-Fock (BHF) definition of U is given by 

{J2u<F < {G{oJau) + G{upy)) \ j3v >, if a and {3 < F 
T.u<F <a^\G{uJau)\pv>: li a<F axid (3> F . (11) 

if a and P > F, 

In this definition a < F refers to single-particle states a below the Fermi surface and the 
starting energies in calculating G are defined by Uah = ^a + ^h using the single-particle energies 

ea =< a\tkin + U\a> . (12) 

For matrix elements of U involving hole states states, i.e. < a|f/|/5 > with a and/or l3 < F, 
it can be shown by a theorem of Bethe, Brandow and Petschek (BBP)[Q, that the on-shell 
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Figure 5: The upper half of this figure displays Goldstone diagrams, which are are included 
in the calculation of the ground-state energy within the BHF approximation (|13|). Diagram 
a) corresponds to the direct, b) to the Fock-exchange contribution. The wiggly lines repre- 
sent the G matrix and the lower half of the figure visualizes the Bethe-Goldstone Eq. (0) 
demonstrating that the G interaction line represents the sum of all particle-particle ladder 
diagrams. This implies that e.g. the Goldstone diagram displayed in c) is redundant as it is 
contained already in a). 



definition of the starting energy in the BHF choice for U (|TT]) yields an exact cancellation of 
many diagrams of higher order in G. The BBP theorem cannot be applied to the particle- 
particle matrix elements of U {a and /3 > F). Therefore the choice = for particle states, 
the so-called conventional choice, has been favored in many BHF calculations. 

Looking at the Eqs. (|^), (|TT|) and (0), which define the BHF approach, one can see that 
these equations request the solution of a self-consistency problem: in order to solve the Bethe- 
Goldstone equation one has to know the single-particle states \a > and single-particle ener- 
gies ea, to define the Pauli operator Q and starting energies, respectively. On the other hand, 
one should know the G-matrix already to determine the single-particle states and energies by 
diagonalizing tkin + U . This self-consistency requirement goes beyond the usual self-consistency 
requirement, which request the knowledge of the hole states to define U . A self-consistent 
solution of the BHF equations can be obtained by solving the Bethe-Goldstone equation and 
the Hartree-Fock equation in an iteration scheme. After the self-consistency has been achieved, 
it is easy to calculate the total energy by 

Ebhf = t: ^ tu + , (13) 

the sum of single-particle energies e^, and the expectation value of the kinetic energy t^, for 
the single-particle wave functions of the hole states. Goldstone diagrams representing this 
approximation for the energy are displayed in Fig. ^. Results for other observables like e.g. the 
radius of the mass or charge-distribution are usually determined by simply calculating the 
expectation value of the corresponding operator for the model space Slater determinant $. 
Strictly speaking one should derive an effective operator also for such observables, which account 
for the restriction of the complete Hilbert space to the model space. 



10 



The optimal choice of Hq to be used for the particle states in the BHF definition of f/, 
and consequently also the hu in the Bethe-Goldstone equation, has widely been discussed in 
particular for the case of nuclear matter calculations [0, Good arguments have been 



presented to favor a single-particle spectrum which is continuous at the Fermi surface and 
therefore contains an attractive potential for the low-lying particle states p8| or all particle 



states 153] • The answer to the questions, which is the optimal choice, can only be obtained by 
evaluating the contributions of higher order correction and their sensitivity to the choice of U. 

So before we can answer the question, how to choose U in an optimal way, we should discuss, 
how to go beyond the BHF approximation, which is also often called the lowest order Brueckner 
theory. The total energy is calculated in the BHF approximation (see Eq. (pTSD) by taking into 
account the contribution of all ladder diagrams with any number of intermediate two-particle 
states, which are summed up in the G-matrix. This means that all diagrams with two hole lines 
are taken into account. The first step beyond this two- hole line approach would be to include 
the contributions of all linked diagrams with three hole lines. This ordering with respect to the 
number of hole lines, which is the basic assumption of the hole line expansion, can be justified 
with the following argument: Linked diagrams including n hole lines describe processes, in 
which n nucleons are interacting in a coherent way. This means that all n nucleons should be 
within a volume that they can all interact with each other. If the range of the strong interaction 
rv is smaller than the average distance to the next neighbor d^, the probability that n + 1 
nucleons are found within a volume of mutual interaction is smaller than the corresponding 
probability for n nucleons by a factor {ry/d]\fY. In nuclear matter around saturation density, 
the average distance to the next neighbor is around 1.8 fm, which is larger than range of the 
strong short-range components of the NN interaction (the radius of a hard core is around 0.4 fm, 
the range for the exchange of an uj meson, h/nc, corresponds to 0.26 fm), however, comparable 
to the range of the pion exchange (1.45 fm). It is obvious that this hole line expansion should 
work at small densities, but may fail at large densities, at which the assumption, ry small 
compared to dj^, is not justified. The convergence of the hole-line expansion will be discussed 
in section 3 below. 

The inclusion of all two-hole line contributions requires the solution of the two-body problem 
in the nuclear medium, as expressed by the Bethe-Goldstone equation. For the inclusion of 
all three-hole lines one has to solve the three-body problem in the nuclear medium, which 
corresponds to solving the Bethe-Fadeev equations ||58| (see also diagrams in Fig. H). Techniques 
and details how to solve the Bethe-Fadeev equations have been described by Day||S^. Numerical 
solutions for the three-hole line contributions in nuclear matter have recently been presented 
by Song et al.0, 0. 

Various techniques have been developed to solve the BHF equations, i.e. the Bethe-Goldstone 
equation in particular for nuclear matter and finite nuclei. As an example for a standard so- 
lution in infinite nuclear matter we would like to refer to the work of Haftel and Tabakin 



A description, how to solve the Bethe-Goldstone equation for finite systems, including a FOR- 
TRAN program, can be found in [B3]. 



2.2 Coupled Cluster Method 



A very detailed description of the coupled cluster Method (COM), which has been introduced 
more than 40 years ago by Coester and Kiimmel jH^, |65[| can be found in [|l^. More recently a 
review on this approach has been given by Bishop |13|. Also in the COM or "Exponential S" 



approach one starts assuming an appropriate Slater determinant $ as a first approximation for 
the exact eigenstate of \1/ for the A-particle system. If the overlap < \E'|$ > is different from 
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Figure 6: Three-particle ladder diagrams to be included in the Bethe-Fadeev equation. 



zero, one can always consider the "Exponential S" ansatz 



with S an operator of the form 



(14) 



(15) 



n=l 



where Sn is an n-particle operator which can be written 

1 



E < pi . . .p„|5'„|z/i . . . z/„ > a 



(16) 



Pl-Pn 



Here and in the following a]^ (oq,) stand for creation (annihilation) operators for nucleons in a 
state a with a a single-particle state of the basis, which also defines the Slater determinant 
We will use labels i/i . . . z/„ to refer to hole states, i.e. single particle states which are occupied 
in the model state $ and indices pi . . . p„ to denote particle states, i.e. states above the Fermi 
surface of $. Expanding the exponential in (0) one obtains 



^ = {l + Si + ^Sf + S2 + ^Sf + ^2^1 + ^3 • • • ) $ • 



(17) 



So exact eigenstates of the hamiltonian is written as a sum of the reference state one- 
particle on-hole (5*1) excitations relative to $, two-particle two-hole (1/25"^ -|- 5*2) and so on up 
to A-particle A-hole excitations. Therefore it is obvious that the ansatz (|14D is sufficient, but 
one may ask the question, why we are using the exponential form of the ansatz and not simply 
write 

= (1 + Fi + F2 + F3...)$ 

= {l + F)<l>, (18) 
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where F„ contains all ra-particle n-hole contributions, which implies 

Fi = Si 

= \sl + S2Si + S^ (19) 

If the two expansions are treated including all terms up to n = A both expansions lead to the 
exact result. The question is, which approach is more appropriate if one has to truncate the 
expansion and consider terms or Si only up to an order i < n with n smaller than the total 
particle number A. 

A first answer on the question why the "exponential S" ansatz (|I4D is preferable to the 
parametrisation (|TB|) can be given by Thouless theorem ||66|]: If we restrict the "exponential 
S" ansatz including only terms up to = 1, the Thouless theorem says that the ansatz ([T^ ) 
includes all Slater determinants which are not orthogonal to $. Therefore the solution of the 
CCM method in the Si approximation would correspond to the Hartree-Fock approach. On 
the other hand, however, the Fi approach, i.e. ignore all contributions of Fi with i = 2 ... A in 
([T8|), is more restrictive and more sensitive to the choice of the initial state $. 

A more convincing argument, however, can be found by considering the fact that e.g. F2 
will contain two-particle two-hole excitation which are unrelated and completely independent 
from each other. Just from statistical arguments it is clear that such unlinked contributions will 
become more and more important if the size of the system, or the particle number gets large. 
This problem can be seen by inspecting the set of equations, which determine the amplitudes 
Fi. To derive these equations one considers the Schrodinger equation 

if|^>=E|^> (20) 

assumes the ansatz ([T8|) and projects from the left with < $|, < $|aj,ap, < ^\al_^al^ap^ap-^ and 
so on. This yields 

E = < <l>|if(l + F)|<1> > 

E<p\Fi\u> = < <l>\aiapH{l + F)\<^ > 
E < piP2\F2\uiU2> = <^\alap,ap,H{l + F)\<^> (21) 

One finds that all left-hand sides of these equations contain the energy, which is an extensive 
quantity, i.e. it grows proportional to the particle number A. Therefore also the right-hand 
side of these equations must be extensive, which is a reflection of the fact that these equations 
contain unlinked terms. 

The "exponential S" ansatz provides a way out of this problem. This can be seen already 
from Eqs. ([T9|) , in which e.g. F2 is rewritten in terms of unlinked contributions Sf and a linked 
term 5*2. To demonstrate this feature we consider the equations which determine the amplitudes 
Sn- To do this we consider the Schrodinger equation in the form 

e-^He^\'i> >= e~^E|^ >= E|$ > (22) 

and project it from the left with < $|, which yields 

E = < $|e-^i/e^|$ >=< (!>\He^\<^ > 

= < <l>|if(^l + Si + ^52 + ^2) 1$ > (23) 
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Note that in the first hne we have used the fact >= 0, which is obvious from the definitions 
of the Sn in (0). The other equations are obtained by muhiplying from the left with 
< $|aj,ap, < ^\al_^al^ap^ap-^ etc. which yields 

< <l>|4ape"^iJe^|$ > = 
< ^\al^al^ap^ap^e~^He'^\^ > = (24) 

and so on. Note that these equations do not contain the extensive quantity E as the corre- 
sponding Eqs. ( |2lD for the Fi. 

The fact, that the Eqs. (|2^) do not contain extensive quantities, growing with A, is only a 
hint that the matrix elements contained in these equations do not contain terms, which corre- 
spond to unlinked diagrams. In order to proof this feature, one has to consider the evaluation 
of the matrix elements more in detail. For that purpose we consider the expansion 

e-'He' = H+[H,S] + ^^ [[H, S] , S] + [[[H, S] , S] , S] + . . . (25) 

As an example for such commutator expression let us consider the commutator between the 
operator of kinetic energy with 5*1 

[tkin,Si] = ^ < a\tkin\P >< p\Si\l^ > 
afipv 

= H ^ H < «|4m|p >< p\Si\v > ala^ \ '^\tkin\P >< p\Si\iy > aja/? [> 

Each commutator removes one and one a operator from the operator product, and links 
single-particle labels of two amplitudes together. Translated into the language of diagrams 
such a link means that there is an internal line, a single-particle propagator, connecting the 
operator symbols of t^m and Si. The components of S commute with each other since they 
all contain creation operator for particle states and annihilation operators for hole states in 
Therefore links produced by the commutators in (^) can only occur between the hamiltonian 
H and S. This implies that the expansion in (|25| ) is finite: if e.g. the hamiltonian contains 
one- and two-body operators only, there will only be up to four nested commutators on the 
right hand side. Furthermore we see that the left-hand sides of Eqs. correspond to matrix 
elements of a linked operator product between the unperturbed model ground-state $ as ket- 
and n-particle n-hole states relative to $ as bra-state. 

In order to obtain explicit equations for the amplitudes Sn one introduces the so-called 
n-particle subsystem amplitudes defined by 

< rci .. .x„|^„|i/i ... z/„ >=< ... 4^ aa;„ ... a^.J^ > . (27) 

So this subsystem amplitudes correspond to the overlap of the exact A-particle wave function 
^ with an n-particle n-hole excitation of the reference state $. Note, however, that the label 
Xi does not necessarily refer to quantum numbers of a single-particle state in $, it may also 
refer to coordinates of a nucleon in the usual space or in momentum space. As an example we 
consider explicitly 

< Xi|\l'i|z/i > = < a;i|z/i > + < xil^ilz/i > 
< XiX2|^2kl«^2 > = ^{< xi|^i|z/i >< X2|^l|z/2 >}+ < a;iX2|S'2|z/iZ/2 > (28) 

If we consider xi to represent coordinate, we see that the single-particle amplitude contains 
the uncorrelated wave function of the hole state < xi\vi> and possible corrections defined by 



(26) 
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5*1. In a similar way \l/2 contains the anti-symmetrized [A represents the operator of antisym- 
metrisation of the indices z/j) product of the plus corrections due to 5*2 This would imply 
that < xiX2|5'2|i^iz/2 > plays a role similar to the defect function in ([101) . For the three- and 
more-particle amplitudes it is convenient to introduce also[p7[ 

< XiX2p\xs Wi^2^3 > = < XiX2p\^3Wi^2^3 > -A{< pl^^il^s >< a;ia;2 1 5*2 1 Z/1Z/2 >} 

= ^{< X2|^i|z/2 >< a;ip|5'2|z/iz/3 >} 

+^{< Xi\'^i\Ui >< X2P\S2\V2V^ >} 

-I- < XiX2p\S-i\viV2l''i > . (29) 

Furthermore we can define a single-particle potential for hole states u by 

^x\U\u>= <xv'\V^2W > ■ (30) 

u'<F 

If we identify \E'2 with a correlated two-particle wave function like in ([10|) , the product V'^2 plays 
the same role than the G matrix in the hole-line expansion and the single-particle potential of 
([30| ) can directly be compared with the definition of the single-particle potential ( p!TD in the 
BHF approach. 

Using all these definitions the one-particle, Hartree-Fock like equation, determining Si or 
can be written 

< Xi\tkin,l^ > + X! < xW\tkin,2S2\Vlv' > 
y'<F 

+ <Xi\U\ui> + < XiVv'\V2zxf\iyiiyp' > = ^ <a;i|^i|z/>/i,,, (31) 

uv'<F u<F 

with h,yui the matrix elements for single-particle energy between hole states 

=< y\tkin,i'^ i^i > + < y\U\vi > (32) 

The indices 1, 2 in these equation shall denote that operators of kinetic energy and the two-body 
interaction V act on the corresponding single-particle states in the ket-vector. 

The corresponding equation for the two-body amplitudes, determining 5*2 or \E'2 can be 
written (omitting some corrections proportional to Si) 

< xiX2\Q{tkin,i + tkin,2)S2\i^ii^2 > " XI (< 2:1X2 1 ^2 1 z/i^2 > K^^+ < a;ia;2 1 5*2 | Z^lZ/ > K^^) 

v<F 

< XxX2\QVi2S2V'\V2> = - < a;ia;2|Vl2|z^iz^2 > - < a;iX2|5'2PVi2\E'2|'^i'^2 > 

- X < x,X2u\QVnX^'^ + QV23xt^Wi^2iy > 

u<F 

~\ < xiX2iyiy'\QV3iX4'^^Wi^2^^' > ■ (33) 

Here Q refers again to the Pauli operator for two-particle states, which we have defined already 
in (^. From these two equations (^Tj) and ( P^D we can observe some general features of the 
problem to determine the amphtudes Sn- The different equations are coupled. If we restrict our 
considerations to two-body interaction terms only, the n-body equation, which should provide 
Sn requires the knowledge of Sn+i and Sn+2- If the hamiltonian would also contain a three-body 
interaction, Sn+3 would be needed as well. In order to solve the one-body equation (|3TD which 
will determine 5*1 or \E'i, one should know the amplitudes S2 and 6*3, which are hidden in U, "$2 
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Figure 7: Hole-hole ladder diagram included in CCM SUB2 approximation 



and Xs, respectively. For the two-body equation (|33D the information on 5*3 and 5*4 (contained 
in Xi) is needed. 

This implies that one has to truncate the hierarchy of this set of equations by assuming 
in the so-called "SUBn" approximation that all amplitudes Si with i > n are assumed to be 
zero. As a first example we will consider the SUB2 approximation and compare it to the BHF 
approach discussed in the preceeding subsection. The SUB2 assumption implies that we ignore 
the last term on the left-hand side of Eq. (|3TD (xs) and the last two terms in Eq. (|33D. If for 
the moment we furthermore assume that 

hu,, = eJu^ (34) 
is diagonal and neglect the matrix elements of 5'2PVi2^E'2 we can rewrite Eq. ( |33D into 

Q (e^^ + 6^,2 - {tkin,l + tfcm,2)) 32^1^2 >= QV'^2Wl^2 > (35) 

which corresponds to 

5'2|z/iZ/2 >= — ^ — -¥-^2^11^2 > ■ (36) 

Comparing with Eq. (|T0|) we may identify 

!^i!^2 > ^ ^ V'^2\^l^2 > 

— ^— (710^,^2 > < — > 82^11^2 >, (37) 
a; - hi2 

where the last line identifies the so-called defect functions in the two approaches. With this 
identification we also see that the single-particle equation (^) corresponds to the single-particle 
equation if we ignore the term tkin,2S2 in that equation. 

The main difference between the Coupled Cluster Method (CCM) SUB2 approach and the 
BHF approximation is the inclusion of the 5*2^1^12^2 in the two-body Eq. (|33D as compared to 
the Bethe-Goldstone equation (|^. In terms of diagrams this means that the energy calculated 
in the CCM SUB2 approach also includes hole-hole scattering diagrams like the one displayed 
m Fig. 0. This means that the hole-hole ladders are treated on the same level as the particle- 
particle ladders. This should be of particular importance if the phase space of hole-hole states, 
which can be reached by the two-body interaction is as large as the corresponding space of 
particle-particle configurations. As we will see in the discussion of results in chapter 3, for 
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realistic interactions in nuclear physics, the particle-particle ladders are dominant as compared 
to the hole-hole ladders. Therefore the differences between the results of CCM SUB2 and 
BHF calculations are usually small. Note that particle-particle and hole-hole ladders are also 
often referred to as particle-particle hole-hole ring diagrams, which are included in a pphh RPA 



calculation p8|, |69| , |70| 

Within the framework of the CCM approach one does not have any choice for the auxiliary 
potential U as it was the case in the BHF approach. From the comparison of the two approaches 
discussed above (see Eq. (|37|) ) one finds that the CCM SUB2 approach is rather close to the 
BHF assuming the conventional choice for the particle state-spectrum in the Bethe-Goldstone 
equation, which means that hi2 is replace by the sum of the kinetic energies. All corrections to 
the CCM SUB2 approach occur due to the inclusion of the three-body and higher order terms. 

In order to discuss the effects of the three-body amplitude on the calculated energies, we 
write the equation for Xs^l 



+ < X2X3\S2\iy2l^3 >< a;i|\£'l|z/l >} 

+ J2 < XlX2X3\ — \y2y3yi >< l/l2/2l/3|G'i2X3^^Vl'^2i^3 > 

+ ^ < xiX2X3\ — \y2y3yi >< ymy2\Gi2xT\^i^2T^^ > (38) 

In this equation we have identified ¥"$2 with G according to (|37|). Q3 denotes the Pauli operator 
for three-particle states and the energy denominator 63 is defined by 

63 = ^ui + ^U2 + ^u-i — (tkinl + tkin2 + tkins) ■ (39) 



The resulting amplitudes X3 multiplied with V enter into the one-body ( PT| ) as well as two-body 
Eq. (pSl). The contribution to the binding energy E originating from these contributions is of 
third and higher order in G. The Goldstone diagrams representing the contributions of third 
order (omitting exchange diagrams) are displayed in Fig. ^. The diagram of Fig. ^a is a particle- 
hole ring diagram. The summation of all ph ring diagram contributions can be obtained by 
calculating the correlation energy arising from particle-hole RPA calculations. The diagram 
displayed in Fig. §3 is a contribution, which within the BHF approach one would try to cancel 
by an appropriate definition of the single-particle spectrum U for particle states. This indicates 
again that U should be chosen to minimize the effects of three- and four-body terms. 

If one would like to evaluate the amplitudes < Xi . . . Xn\Sn\i^i . . . Un > in the real space 
representation, which means that Xi refer to the usual coordinates of the nucleons, the basic 



equations (|3l|) and (^) must be considered as differential equations. This is particularly useful 



if the NN interaction V is described in terms of local potential terms p7|, [7I| . For non-local 
potentials it is me more convenient to use momentum space representation. In this case the 
equations lead to inhomogeneous integral equations|]72|. The CCM equations have also been 



solved in a Hilbert space, which is spanned by a basis of appropriate oscillator functions ||73 
which leads to a solution of a coupled system of nonlinear equations. 

If the amplitude Si have been determined, it is easy to evaluate the energy according to 
([23|). Knowing Si we also have the information on the exact wave function (not normalized). 
An efficient evaluation of observables different from the energy can be done using the following 
considerations. As an example we consider the single-particle density matrix 

dpa = ^^rr 40 
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a) b) 



Figure 8: Lowest order contributions originating from the inclusion of three-body terms in 
the CCM. The diagram shown in the left part (a) is a particle-hole ring diagram, while the 
one displayed in (b) contains a bubble insertion into a particle line, which one would like to 
compensate by an appropriate choice of the single-particle potential for particle states. 



One can rewrite this expression into 
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<*|e-«aXe-|4>+i:^ ^° (41) 



The last line has been obtained by inserting the unity operator 

\<^><^+Y,Xi\^>< Phi\Xn (42) 

n 

with the n-particle n-hole operator 

^n = 7^ E (43) 
V'--^ p^>F,u,<F 

Since commutes with S, the expression (|4l|) can be rewritten into 

d^a =< $|e-^4a„e^|<l> > +J2 < ^l^^e-^^a^e^l^ > ^ ^1^^^ ^ • (44) 

71 < W I W > 

The matrix elements of operators like e~ be calculated employing an expansion 

similar to the expansion in (pSl) . The factors of the form < \E'|X^|\E' > / < \E'|\E' >, on the other 



hand correspond to the matrix elements of the n-body density operator, which brings us back 
to the starting point in (^Ol). This means that the final matrix elements can be calculated in 
an iterative scheme |75|, 

2.3 Many-Body Theory in Terms of Green's Functions 

The two-body approaches discussed so far, the hole-line expansion as well as the CCM, are es- 
sentially restricted to the evaluation of ground-state properties. The Green's function approach. 
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which will shortly be introduced in this section also yields results for dynamic properties like 
e.g. the single-particle spectral function which is closely related to the cross section of particle 
knock-out and pick-up reactions. It is based on the time-dependent perturbation expansion 
and also assumes a separation of the total hamiltonian into an single-particle part Hq and a 
perturbation Hi as introduced in (H). A more detailed description can be found e.g. in the 
textbooks of Fetter and Walecka [[75[] , Negele and Orland|7^ or in the book by MattuckfT? 



which particularly provides a rather intuitive interpretation of the Feynman diagrams. A very 
comprehensive description of the main features has been presented by Mahaux and Sartor [[78[|. 
Introductions are also given in various review articles . 
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The expectation value of any operator O is calculated in the so-called interaction picture as 

<^>\t)\0,{m\t)>. (45) 

where the time-dependent operator Oi{t) in the interaction picture is related to the operator 
O in the usual Schrodinger picture by 

Oi(t) = e'^^'Oe-'^^' . (46) 

and the time-dependence of the state |\I'^(t) > can be derived from 

l^^(t) >=Hn{m\t) > . (47) 

If we introduce the evolution operator in the interaction scheme Ui{t,to) by 

\¥{t) >=Ui{t,to)\^\to)>, (48) 

the equation of motion for the state |\E'^ > ( ^T] ) can be rewritten into a differential equation for 
the evolution operator 

t^^Ui{t,to) = Hnim{t,to) . (49) 
This equation can be transformed into an integral equation 

r* 

'to 



Ui(t, to) = Ui{to, to) - % t dti Hn{ti)Ui{ti, to) . (50) 

Jtn 



Using the fact that Ui{to,to) is identical to the unit operator 1, we can iterate this integral 
equation to a perturbative expansion in powers of Hn: 



Ui{t,to) = i + H) I dtiHniti) 

J to 

+ {-if f dti ifn(ti) r dt2 Hn{t2) + ... (51) 

J to J to 



The integration variables used on the left-hand side are nested in a rather inconvenient way. 
Therefore one rewrites the term of order n in this expansion using 

ft ptn-l \ ft ft 

/ dti Hniti) ... dtn Hn{tn) = dti . . . t^T (Hniti) . . . HM) (52) 
Jto Jto n\ Jto Jto 

with the time ordering or chronological operator T, which is defined for two operators by 

r {A{ti)B{t2)) = I ^^l]^B{t2)A{t,), Otherwise. ^^^^ 
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Here m is the number of exchanges of fermion creation and annihilation operators contained in 
A and B, which are needed to bring A and B into chronological order. Note that for our present 
purpose (0) the factor (— 1)"^ is always equal to 1, as the number of fermion operators defining 
Hii is even. The definition of T in ( ^^D for two operators is easily extended to n operators. 
Applying ( |52D to the expansion of the time evolution operator in (|51|), one gets 



Ui{t,to) = J2^^ / dh... dt^r{Hn{h)...Hn{t^)) (54) 
^ nl Jto Jto 



In order to arrive at a perturbation expansion for the calculation of matrix elements, one 
assumes that the perturbation Hi is "switched off" at times t = —oo and t = +oo and can be 
switched on in an adiabatic way for times t ~ 0. This can be achieved by a time- dependent 
hamiltonian of the form 

Hait) = Ho + e-^^'^Hi (55) 

where a is a small positive number that becomes infinitesimal in the adiabatic limit. This 
procedure implies that the eigenstates of the hamiltonian are identical to the eigenstates 
of the unperturbed hamiltonian Hq at times |t| = oo and should evolve to the corresponding 
eigenstates of the total if at t ~ if we use the time evolution operator Uia{t, t') for the hamil- 
tonian Ha- If we denote the (nondegenerate) ground state of Hq by $0, which is independent 
of time in the interaction picture, this means that we obtain an eigenstate of the exact H at 
time t = by 

|^o(t = 0)>=lim[/i„(0,-oo)|<l>o> . (56) 



It has been shown by Gell-Mann and Low[^ that \E'o is indeed an exact eigenstate of H if the 
perturbation expansion converges. 

In order to calculate matrix elements we now consider the Heisenberg scheme. In this 
representation, which corresponds to the interaction scheme with Hq = H, the wave function 
is time-independent and identical to "^oit = 0) and the time-dependence is completely assigned 
to the operators Onit) which are defined as in (^Bj), replacing Ho by H. This means that a 
matrix element of an operator O can be calculated as 

lim ^ '^ol^i^(QO^O)QH(t)^ia(0,-oo)|$o > 
< $o|f/i„(oo,0)?7i„(0,-oo)|<I)o > 
<%\Uio,{oo,t)Oi{t)Uio,{t,-oo)\%> 

< <I)o|f/la(oO,-Oo)|$o > 



< ^0 




< ^0 


^0 > 



Using the explicit representation of the time evolution operator in (|54|) one can furthermore 
show (see, e.g.|]75[) that the matrix element for any time-ordered product of two Heisenberg 



operators can be calculated as 

<*o|r(AH(t)5H(t'))l*o> 



lim 



< ^ol^o > ^0 ^- -^-^ 



< ^ol^o > 

X < <I>o|T (Vl(ti) . . . Vi(t„)^i(t)5i(0) \% > 



(5J 



which means that one has to evaluate matrix elements of time-ordered products of operators 
in the interaction scheme for the ground state of the unperturbed hamiltonian $o- 

The operators in these matrix elements are time-ordered. Therefore we can apply Wick's 
theorem to evaluate them. In the following we will shortly review this scheme and illustrate 
how the various contributions are visualized in terms of Feynman diagrams. 
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First we recall that the operator for the residual interaction Hi as well as any other operator 
can be expressed in terms of the basic single-particle creation (a|) and annihilation operators 
{tti). It is easy to verify that these operators in the interaction scheme are given by 

aij(t) = ajexp{—iejt) , 

a\j(t) = a] exp (+zejt) , (59) 

with ej the single-particle energies defining the unperturbed hamiltonian Hq The ground state 
for this unperturbed hamiltonian, $05 is given by a Slater determinant in which all single- 
particle states i with an energy below the Fermi energy ep are occupied. The Fermi energy 
separates hole-states {e^ < ep, occupied in the unperturbed ground state) from particle states 
(ep > ep, unoccupied in the unperturbed ground state). If M, N,0, P, . . . represent creation 
or annihilation operators, in the Schrodinger or in the interaction picture, one can define the 
normal product of such operators by 

Af{MNOP...) = {-iyOP...MN... (60) 

where the sequence of operators on the right-hand side of this equation is such that all creation 
operators for particle states (aj,) and annihilation operators for hole states (a^) are moved to 
the left {0,P), whereas all creation operators for hole states (aj,) and annihilation operators 
for particle states (op) are moved to the right (M, A^) and 7 counts the number of exchanges 
of these operators required to obtain the normal ordering. This normal ordering guarantees 
that the matrix element of such an ordered product calculated for the unperturbed state $0 
vanishes: 

< $o|Ar(MiVOP...) |$o >= 0. (61) 

Furthermore we define a "contraction" of two operators, using the chronological operator of 
( |53| ) for two such operators M, N in the interaction scheme, as 

MiV = < $o|r (A'/A^) -Ar(MA^) |$o > 

= <<!>o\T{MN)\%>, (62) 

which is just a complex number given by 

'^ifcC"*'^^ if j is a particle state and t > t', 

aij{t)al^^{t') = { -S-f^e-'^At-t') if j is a hole state and t' > t, (63) 
otherwise, 

al{t')ai,{t) = -a,,it)alit') . (64) 



Using Wick's theorem it is easy to show that matrix elements of time-ordered products calcu- 
lated for the unperturbed ground state $0 as in ( p^ can be calculated as a sum of all terms 
in which the single-particle operators are pairwise contracted, which means that each pair is 
replaced by the corresponding value for the contraction (pH]). 

The use of Feynman diagrams provides an easy control of all the contributions. To demon- 
strate this, we consider as an example the matrix element occurring in the first-order term of 

m 

- ^ < ij\V\kl >< <l>o|T [aUti)a\j{ti)aii{ti)aik{ti)aia{t)alp{t')) |$o > , (65) 

i,j,k,l 

where we have made the substitution Ai aia, Bi a|^, and used the explicit representation 
of V. For this example we will furthermore assume that t' <ti < t. As before the operator of 
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d) 



e) 



Figure 9: Building blocks for Feynman diagrams (a) and completely contracted contributions 
as discussed in the text. 



the residual interaction is represented in the Feynman diagram by a horizontal dashed line, with 
an outgoing and incoming arrow at each end. The outgoing arrows refer to creation operators 
contained in V (ajj, a^) in our notation, and the incoming ones refer to annihilation operators. 
The external operators are represented by a dot with an ingoing and an outgoing arrow for the 
annihilation and creation operator, respectively. These objects are displayed in Fig. With 
the assumption that the vertical axis represents a time axis the objects are ordered according 
to the choice of our example: t' <ti < t. 

Any contraction implies that two of the lines with arrows must be paired. In the graphical 
representation this is achieved by connecting them. Looking at the contractions which yield 
results different from zero (see (|63|)-(|6^), one finds that only those pairs of lines in which 
the arrows point into the same direction must be connected. Furthermore we can distinguish 
connected lines with an arrow pointing upwards, which refers to a particle state, i.e., the 
corresponding summation indices in (|65| ) can be restricted to particle states, whereas connected 



lines with an arrow pointing downwards refer to hole lines. 

All diagrams representing the completely contracted terms of the expression shown in (|65| ) 
are displayed in Fig. |^)-e). Note that we show only those diagrams that are topologically 
distinct in the sense that a diagram obtained from another one by just mirroring the ends 
of an interaction line is not displayed again. Recalling Wick's theorem, it is evident that all 
non- vanishing contributions to the perturbation calculation of the nth-order term in the matrix 



element of a time-ordered product of operators in (|58|) can be obtained in the following way. 



Draw n interaction lines (dashed lines in Fig. ^) and mark on the creation and annihilation 
operators for the external operators to be calculated (dots in Fig. H). 

Construct all diagrams by connecting the "arrows" linked to these basic building blocks 
according to the rules given in the example of Fig. ^ 

Keep in mind that all possible time orderings of the interaction vertices relative to the 
external operators must be considered (see time integrations in (pSI)). 



Each of the resulting diagrams represents a non-vanishing contribution to the evaluation 



of (|58D and there exist well-established Feynman rules that translate the contribution of 



the diagram into a calculable expression (see, e.g.|77[]). 
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• Consider only the contribution of linked diagrams. The linked cluster theorem, which 
has already been discussed in sections ^TT| and also holds for this case (see e.g.|]T5|). 
For the example displayed in Fig. ^ this implies that only the contributions of the linked 
diagrams displayed in Fig. ^) and e) have to be taken into account. 

Up to this point the diagrams have been used only as a kind of book-keeping tool to identify 
all non-vanishing contributions in the perturbation expansion. However, we have seen already 
that each connecting line in those diagrams represents a contraction, and a line with an arrow 
pointing upwards stands for [see (|6^)] 

aij{t)al{t'\ =< %\T {ai,{t)al{t')) \% > , (66) 

with the creation of a particle taking place before the annihilation {t > t'). This means that 
we can ignore the operator T and rewrite this contraction as 

<^o\aij{t)a\^{t')\^o> = E < '^o|ay(t)|/5></5|alfc(t')|$o > 

= 5jfce""-'^*"*'^ for j a particle state. (67) 

In the first line of this equation we have inserted a summation over a complete set of states 
1/3 > with one particle in addition to the number of fermions in |$o >? in order to show that 
this contraction describes the product of a probability amplitude to create a particle at a time 
t', producing a state (3 and the probability that it is annihilated at the later time t reproducing 
the unperturbed ground state $o- Iii the second line of this equation we have copied the result 
for the contraction from (|63|) that such a propagation of a particle on top of the unperturbed 
state is only possible if j = k refers to a state above the Fermi energy, in order not to violate 
the Pauli principle. 

In a similar way one can convince oneself that a line with an arrow pointing down represents 
the propagation of a hole state, i.e., a particle must be removed first from a state h below the 
Fermi energy before it is put back at a later time. With this interpretation of the contractions 
visualized in the diagrams one can easily interpret the Feynman diagrams in terms of time- 
dependent processes. 

The single-particle Green's function can be considered as a special example of an expectation 
value for the time-ordered product of two operators calculated for the exact ground-state in 
(|5§). It is defined by 

ig{at,(3t') =< ^o\T [ana{t)aif,{t')) |^o > ■ (68) 

Note that here and in the following we have dropped the denominator < ^1/01^1/0 >, assuming 
that the exact ground state is properly normalized. The creation and annihilation operators in 
the Heisenberg representation are defined in an appropriate basis, characterized by quantum 
numbers a and (3. If we assume that a refers to a position r' and /3 to a position r of the 
considered fermion in r space, the single-particle Green's function ig[r't',ft) describes the 
propagation of this fermion from the space-time point [ft) to [r't'). In contrast to the discussion 
of the contractions in the previous section, in this case the propagation is with respect to the 
exact ground state and the complete hamiltonian. 

For a system that is invariant under translation, such as the infinite nuclear or neutron 
matter, which we want to consider, the appropriate set of basis states is that of the momentum 
eigenstates; the Green's function is diagonal in this representation. Rewriting the effect of the 
chronological operator T defined in (^) in terms of step functions B(t — t'), we find that the 
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Green's function is given as 



ig{k,t-t') = Q{t-t') < ^o|aHfc(t)aL(Ol*o > -e(t' - 1) < ^o|4fc(OaHfc(t)|^o > 



0(t-t')Ee"'^''' 

7 

-e(t'-t)Ee 



< ^ 



(A-1) 



(69) 



To arrive at the second part of this equation we have inserted a complete set of eigenstates 
for the system with A + 1 particles {'^^~^^'^) and A — 1 particles {"^s^ ^^), as appropriate, 
and replaced the hamiltonian in the exponential functions of the definition for the Heisenberg 
operators by the corresponding eigenvalues. This means that the energies Eq, and 

refer to the exact energies for the ground state of our reference system, and the exact 



(A-l) 



E, 

energies of eigenstates with A + 1 and A — 1 particles, respectively. Note that the step function 
can be represented by its integral form: 



Q{t) = - lim — 
v^o 2Tci 



duo 

^oo to + ir] 



(70) 



Thus the Fourier transformed Green's function, transforming the time difference t 
energy variable oo, can be written as 
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(71) 



This is the so-called spectral or Lehmann representation of the single-particle Green's function 
Inspecting this equation one finds that the single-particle Green's function is represented in 
terms of quantities that are measurable. It shows poles at energies that correspond to energies 
of the system with one particle added (A + 1) and one particle removed {A — 1) relative to the 
energy of the ground state for the reference system. The residua of these poles are given by the 
spectroscopic factors, i.e., the probabilities of adding and removing one particle with momen- 
tum k to produce the specific state 7 (5) of the residual system. The infinitesimal quantity rj 
shifts those poles below the Fermi energy (the states of the A — 1 system) to slightly above the 
real axis and those above the Fermi energy (the states of the A + 1 system) to slightly below 
the real axis. 

In our spectral representation of the single-particle Green's function (^) we assumed that 
the spectra of the many-body system are defined in terms of discrete energies E!:^"^~^^\ This 
is true of course only for systems confined to a finite space. For infinite systems, the energy 
spectra are continuous and it is more appropriate to introduce the so-called hole and particle 
spectral functions defined by 
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-lmg{k,oo), for < ep 



'£\<^^,''-'^\a,\^o>\5{oo-iE^-E!^^ 
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(72) 
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Figure 10: Classifying the diagrams for the single-particle Green's function in terms of the 
mass operator S and graphical representation of the Dyson equation(|78]). The thick line 
represents the "dressed" single-particle Green's function while the thin lines stand for the 
unperturbed Green's function g^. 



where we have made reference to the case of discrete spectra in the second and fourth lines. 
These definitions imply that the single-particle Green's function is given in terms of the spectral 
functions by 



g(k,uj) = lim / diu' '-"^•"'"-^ + / du' -P'-]- ^ . (73) 

V^O \J-oo UJ — OJ' — IT] Jep OJ — OJ' + IT} j 

The single-particle Green's function or the spectral functions defining the single-particle Green's 
function can be used to evaluate observables of the system. As a first example we mention the 
momentum distribution or momentum density 

du;S^{k,u). (74) 

-OO 

The single-particle Green's function allows the evaluation of the expectation value for any 
single-particle operator O. If we return to the non-diagonal nomenclature used, e.g., in (|68| ) 
and generalize the definition of spectral functions in a corresponding way, such a calculation is 
performed via 

< ^olOl^o >= r '^^ -^hM, oj) < a\0\(3 > , (75) 

with < a|0|/3 > denoting the single-particle matrix element calculated in the basis a, /?.... 
Furthermore, if the interaction between the fermions is a two-body interaction, one can even 
calculate the energy of the correlated ground state via 

1 ^ rep 

= oXl / du Sh{aP,u) {< a\T]^in\P > +ujSa,p) , (76) 

Z r) J — OO 

with Tkin representing the single-particle operator for the kinetic energy. 

What remains to be discussed are the tools and approximations used to determine the single- 
particle Green's function. For that purpose we consider the single-particle Green's function for 
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Figure 11: Diagrams contributing to the self-energy E. Note that the thick connecting hues 
represent "dressed" single-particle Green's functions. This means that expanding e.g. the 
rightmost diagram in the upper line yields among other the terms displayed in the lower 
part of the figure. 



the unperturbed hamiltonian 



e(e„-eF) , e(eF 



^o(a/3, uj) = Sa(3 < — + ^ \ . (77) 

and classify all linked diagrams contributing to the expansion for the final, "dressed" single- 
particle Green's function following the diagrammatic representation in the upper part of Fig. ^ 
In this figure we group all Feynman diagrams into parts, which are irreducible with respect to 
the single-particle Green's function (represented by the ellipses S) and single-particle Green's 
function qq connecting these irreducible parts. The irreducible self-energy or mass operator S 
contains all diagrams of any order in the interaction Hi or V, without an intermediate state, 
which is just a single-particle Green's function. Examples of such contributions are displayed 
in Fig. |Tl]. Note that here we do not distinguish between particle- and hole-propagation, since 
the unperturbed Green's function go (see Eq. (^)) as well as the resulting Green's function 
contains both contributions. 



The series displayed in the upper part of Fig. |T0| represents the Dyson equation and can be 
written 

+ Zlfi'o(^Ai'^)^M^(^)fi'o('^/?>^) + • • • 
= go{a(3,u) + J29o{a'y,uj)i:^5{uj)g{5(3,uj). (78) 

Examples of diagrams to be included in the definition of the self-energy are displayed in 
Fig. |Tl|. The diagram of first order in the interaction V is again the Hartree-Fock contribution. 
Note that it is defined in terms of the dressed Green's function, which is expressed by the 
fact that the loop connected to the V interaction line is drawn as a thick line representing the 
dressed Green's function. This reflects the problem of finding a self-consistent solution of the 
Hartree-Fock equations. The structure of the single-particle Green's function is identical to the 



unperturbed one, defined in ([77|) . The only difference is that Hartree-Fock wave functions and 
energies should be used. 
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If, however, one goes beyond the Hartree-Fock approximation by including diagrams of 
higher order in the definition of the self-energy E, like e.g. the second term displayed in Fig. ^ 
the problem of a self-consistent determination of the Green's function gets much more involved. 
To demonstrate this feature let us consider as a first step this diagram of second order in 
V, replacing the dressed single-particle Green's functions by the corresponding Hartree-Fock 
Green's function qq. This means that the intermediate states in that diagram are given in terms 
of intermediate 2-particle 1-hole (2plh) and 2-hole 1-particle (2hlp) states. The corresponding 
expressions for the self-energy are given by 
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+ irj 



(79) 



and 



p>F ui,V2<F 



UJ — (e 



+ ez.2 - Cp) ir] 



(80) 



If we insert the sum of these two contributions into the Dyson equation (ffSD, we obtain a 
single-particle Green's function, which exhibits a much richer pole-structure. In fact assuming 
a discretized space of single-particle states (as we have done also in Eqs. (^) and (pOD) one can 



rewrite the Dyson equation into an eigenvalue problem ||82 
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The matrix to be diagonalized contains the Hartree-Fock single-particle energies and the cou- 
pling to the P different 2plh configurations, which is described in terms of 

aai =< au\V\pip2 > (82) 

and Q different 2hlp configurations, for which we have introduced the abbreviation 

Aa^ =< ap\V\uiU2 > . (83) 

As long as we are still ignoring any residual interaction between the various 2plh and 2hlp 
configurations the corresponding parts of the matrix in (^) are diagonal with elements defined 
by Ci (Ej) for 2plh {2hlp) 



Ej 6jyj^ ~l~ 6^2 1 

where as before the indices pi and Ui refer to particle and hole states, respectively. 
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mean field approach 



general structure 



Figure 12: Schematic picture of the spectral function in the mean field approach (upper 
part) and the general case. 



Solving the eigenvalue problem (Eq. (^TJ)) one gets as a result the single-particle Green's 
function in the Lehmann representation in the discrete basis of the box defined as in ([7lD. The 
eigenvalues Un define the position of the poles of the Green's function 



Eq^ for Un > Ef 
[e^-EI^-'^) for uJn<EF (85) 



with Ep the Fermi energy for the A-nucleon system. The corresponding spectroscopic ampli- 
tudes are given by 

<<|afcJv^^i>=Xo",^ for cu„ > 

<<|4jvl>^i>=Xo",^ far u.^<Ef (86) 

This multiplicity of the poles is a consequence of going beyond the mean-field approximation. 
Within the Hartree-Fock or independent particle approach, one can either remove a particle 
from the orbit j (if this is an orbit below the Fermi energy) or add a particle into this orbit (if 
j is a state above the Fermi energy). The spectroscopic factor is always one. This is reflected 
in the mean field approximation to the single-particle Green's function by the feature, that it 
exhibits for each orbit j exactly one pole with a spectroscopic factor of one. In the approach 
which we are discussing now this spectroscopic strength gets redistributed as indicated in 
the schematic picture of Fig. The single-particle peak with strength one in the mean field 
approach gets reduced to a quasi-particle peak in the spectral function with a strength Z, which 
is smaller than one. Part of the strength gets redistributed into removal or hole-strength, which 
is predominantly due to the coupling to 2hlp configurations. Another part, however, is shifted 
to energies above Ep and therefore represents a spectroscopic amplitude for adding a particle. 
This implies of course that occupation numbers for the various orbits (see ([7^)) will in general 
be different from zero and one. Also it is not guaranteed any longer that the total particle 
number will be conserved, i.e. identical to A the one of the unperturbed state 

? 

dkn{k) =A (87) 
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In order to obtain an approach, which is particle-number conserving, one has to go beyond 
the approximation which we have discussed sofar and calculate the Green's function in a self- 
consistent way. General rules for calculational schemes of Green's functions, which conserve 
symmetries like the property of a fixed particle number have been developed by Baym and 
Kadanov[p3|]. For our present purpose it is sufficient to say that one has to evaluate the 
self-energy S, which enters the Dyson equation (|78D in a self-consistent way in sofar that 
the resulting Green's function should be used in the calculation of S. This self-consistency 



requirement is displayed in terms of diagrams in the lower part of Fig. [11 . 

The realization of this self-consistency requirement has to face a serious problem: Calcu- 
lating the self-energy S in terms of Hartree-Fock Green's functions led to a dressed Green's 
function which exhibits L + P + Q poles, where L stands for the number of Hartree-Fock states 
of a given symmetry and P and Q refer to the number of 2plh and 2hlp states of this symmetry, 
respectively. In a next step towards a self-consistent solution, one may consider these Green's 
functions with L + P + Q poles already in the calculation of the self-energy. This means that 
one also accounts for coupling to 3p2/i, 3/i2p, 4p3/i etc. configurations. This implies that the 
number of poles in the Green's function resulting from this next iteration step will dramatically 
increase. This "infiation of poles" in the single-particle Green's function can be handled very 



easily by means of the so-called "BAsis GEnerated by Lanczos" (BAGEL) scheme p^, p^ , p5 
In this scheme on represents the Green's function in terms of a few "characteristic" poles in the 
Lehmann representation. The number of these poles can be kept fixed in the iteration scheme, 
leading to self-consistency. These pole are determined by solving the eigenvalue Eq. (|81|) by 
means of the Lanczos scheme, starting with the single-particle states as appropriate initial 
vectors. Another possibility is to keep track of the distribution of the spectral strength in a 
purely numerical way. Corresponding calculations have been performed by van Neck et al. p6| 
assuming a model-space of finite dimension. 

Up to this point we only discussed the determination of the single-particle Green's function 
assuming an expansion for the self-energy, which is self-consistent in terms of the Green's 
function but perturbative in terms of the interaction V . From the discussion in the preceeding 
section we know already that two-body correlations should be included in a non-perturbative 
way, which means that one should consider all ladder diagrams, when one is calculating e.g. the 
self-energy S. How can this be achieved within the framework of the Green's function approach? 

For that purpose one should solve a Dyson equation for the two-particle Green's function. 
The Lehmann representation of the two-particle Green's function is given in terms of energies 
and states of the systems with A and A ±2 particles 



E 



n - {E^+^ - E^) + IT, 



The result for the noninteracting product of dressed propagators, including the exchange con- 
tribution reads 

(a/?,7(5;fi) = i j {g{a,^]u)g{l3,5]VL- u) - g{a,5;u)g{l3,-i]VL- u)} 
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-(7 < > 5) 



(89) 



The integration in the first hne of this equation can be performed by employing the Lehmann 
representation for the single-particle Green's functions. The ladder approximation to the two- 
particle propagator (R3) is then given by: 



1 



(90) 



eriec 



This ladder approximation for the two-particle Green's function can then be used to define 
the self-energy S in a non-perturbative way|14]. The Dyson equations ( pO|) and ( [78D for the 
two-body and one-body Green's function have to be solved in a self-consistent manner, this 
procedure has been named Self-Consistent Green Function formalism (SCGF) and has been 
extensively discussed in Ref. 



14 . 



2.4 Variational Method and Correlated Basis Function 

An efficient way to handle the correlations induced by the NN-interaction is to embody them, 
from the very beginning, in a trial wave function "^t, which describes the system of A nucleons 

^t(1, A) = F(l, A)^mf{1, A), (91) 

where is a mean field wave function corresponding to the uncorrelated system and the 
operator F is intended to take care of the dynamical correlations. 

Once a trial wave function is defined, the variational principle ensures that if we are capable 
to calculate the expectation value of the nuclear hamiltonian 

Et, (92) 





H 






T 





then Et will be an upperbound to the ground state energy, Eq. Parameters in the variational 
wave function are varied to minimize Et and the best \E't can then be used to evaluate other 
observables of interest. Obviously, for the method to be efficient, the trial wave function 
must give a good representation of the real ground state many-body wave function. Although 
conceptually it looks very simple, the evaluation of the expectation value is by no means an 
easy task and very sophisticated algorithms which require large computer capabilities have been 
devised during the last years. 

Therefore the ingredients for a variational calculation are the hamiltonian and the wave 
function. In addition, one requires an efficient machinery to evaluate the expectation value. By 
the definition of the trial wave function, which is given in configuration space, the variational 
method is tailored for a non-relativistic framework where the only constituents are the nucleons, 
interacting through a local NN interaction with no energy or momentum dependence. A realistic 
interaction based on the exchange of mesons between the nucleons would be naturally given in 
momentum space, being non-local and energy dependent p|, There are, however, realistic 
interactions which are specially suitable for variational calculations [^, 

Following these guidelines, one can write a realistic nuclear hamiltonian in the form: 

^ = -|^Ev? + E^.- (93) 

I Kj 
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where the two-body potential Vij has a local operatorial structure, i.e. is given by the sum 
of functions depending on the relative distance between the nucleons and spin-isospin operators 
build according to invariance requirements. Although we devote a full section to discuss different 
modern realistic potentials, it is convenient at this point to make some comments on the 
potentials usually employed in variational calculations. For instance, the Argonne Vi4 0] NN 
potential is given by the sum of 14 isoscalar terms 
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p=l,14 
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Sij — 3((Tj ■ rij){dj ■ Vij) 



(94) 



[l,r,-r,], (95) 



(96) 

and S is the total 



is the usual tensor operator, L is the relative orbital angular momentum 
spin of the pair. 

The radial components of the potential contain a long range part which is given by a static 
nonrelativistic reduction of the one pion exchange potential (OPE) which contributes only to 
the ((Jj ■ cr,)(rj ■ Tj) and Sij^Ti ■ Tj). The intermediate and short range parts are given in terms 
of a physically plausible parameterization with a reasonable number of adjustable parameters. 

However, these accurate NN potentials do not satisfactorily reproduce the nuclear matter 
saturation point and underbind nuclei with A > 2 |T^. As discussed in the introduction, a 
possible reason for that is the existence of three body forces which can have different origins. 
We will come back to the problem of three body forces, but here it is important to keep 
in mind that the nuclear hamiltonian can certainly have an additional term with three body 
forces. Usually the recent realistic variational calculations include also a three body force in the 
nuclear hamiltonian, which is constructed in largely phenomenological fashion with parameters 
determined to reproduce the saturation of nuclear matter and the correct binding energy for 
A = 3,4 nuclei. 1^ 

There is also an updated version of the Argonne potential which breaks charge indepen- 
dence and charge symmetry and contains some additional isotensor and isovector components, 
responsible for the breaking of the isospin symmetries. This version of the two-body Argonne 
potential contains in total 43 parameters and an operatorial structure up to 18 operators, it is 
known as the Argonne fig NN potential |^ . 



Once the realistic hamiltonian has been defined, which for the time being we suppose to 
have only two body forces, we need an ansatz for the variational wave function. At this point, 
we should appeal to the physical intuition and sense in order to choose an appropriate trial wave 
function for the system under consideration. In the nuclear case, i.e. for nuclei and for nuclear 
matter around saturation, the density is small enough to assume that two body correlations 
will be the most relevant ones. A form of F{1, 2, A) that has shown to be suitable for nuclear 
systems is 



F{1,...,A) = S 



i<j 



(97) 



i.e. a symmetrized product of two-body correlation operators, F2{i,j). The natural choice is 
to allow the two-body correlation operator to have a similar structure as the two-body NN 
interaction. In recent calculations, the ansatz for the two body correlation operator F contains 
the same set of operators as in Eq. ( P3| ) except the quadratic terms in the relative angular 
momentum, 

F2{^,J)= E f^"'\n^)0l7\ (98) 
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where the sum runs up to the spin-orbit components in Eq. (p5|). 

When the components of the correlation operator with m > 2 are disregarded, one recovers 
the well known Jastrow correlation function |Q , which has been largely used in the context 



of quantum liquids E^ and also in nuclear physics for simple semi-realistic interactions 

The mean field wave function ••,^) is a Slater determinant of single particle wave 

functions, fai'i), where subscript a stands for the set of quantum numbers characterizing the 
single particle state, and (i) indicates the spatial and spin-isospin variables of particle (i). 
Those single particle wave functions are obtained by some mean field (MF) potential, which 
can be completely arbitrary. In the particular case of symmetric and spin and isospin saturated 
nuclear matter, the mean field wave function is a Slater determinant of plane waves: 

<I>mf(1, 2, A) = DeU^ [exp{ik, ■ r,)xf (j)] , (99) 

with all single particle momentum states kj occupied up to the Fermi momentum kp = 
(67r^p/(i2)^/^, where d is the spin-isospin degeneracy of each particle level, d = A for nuclear 
matter and d = 2 for neutron matter, xf (j) cind xIU) ci-re the spin and isospin functions. The 
mean field wave function ($mf) is intended to carry the correct quantum statistics and any of 
the required symmetries of the system. As the mean field wave function has already the correct 
antisymmetric character, the correlation operator requires the presence of a symmetrizer opera- 
tor S (Eq. (^)) to preserve the correct symmetry behavior of the trial wave function. Obviously 
in the case that the different two-body correlations commute with each other, as it would be 
the case for a simple Jastrow correlation, the presence of the symmetrizer is superfluous. 

Besides being symmetric in all the variables, the correlation operator possesses the cluster 
decomposition property, namely that upon separating one subgroup of particles (1,2, ...,n) far 
from the rest {n + l,n + 2, A), the operator F(l, 2, A) decomposes into a product 

F(l, A) = F(")(l, n)F^-"(n + 1, A). (100) 

According to the clusterization property, the scalar component f^^^^\r) in (|98|) heals to unity 
when r is large whereas f^P^^^{r) — > 0. 

A realistic trial wave function is expected to give an upperbound very near to the ground 
state energy and a wave function close to the ground state wave function. The radial correlation 
functions /('")(r) are determined by minimizing the energy, Et- However, the minimization of 
the expectation value of the hamiltonian respect to arbitrary variations of /("*)(r) 

(101) 

is a highly prohibitive job in the nuclear case. In practice, the correlations are either param- 
eterized and the parameters fixed by minimization or they are generated by solving a set of 
coupled differential equations obtained by minimizing the energy obtained in a two-body cluster 
expansion. 

Once the hamiltonian and the trial wave function are defined, one should be able to calculate 
the expectation value. The most efficient techniques are Fermi-Hyper-Netted- Chain (FHNC) 
theory and Variational Monte Carlo (VMC) method (to be discussed in the next sub-section). 
FHNC is an integral equation method that sums up series of clusters diagrams associated with 
the distribution functions of the many-body wave function. The method is suitable for infinite 
homogeneous systems as well as for finite, inhomogeneous systems. In fact, the expectation 
value of the hamiltonian (or any other operator) is written in terms of n-body densities, 

p, = £./(■■-■■■) I ^^>, (102) 



5 (^'t 


H 






T 





32 



(103) 



The densities are then cluster expanded in terms of dynamical correlations h(r) = /'•^•'(r 

products f^^\r)f^P-'^\r), and statistical correlations, Poihj) = Y.a 0l(^)0a(j)) associated to the 
exchanges in the mean field wave function. The terms of the resulting expansions are usually 
called cluster terms and are characterized by integrals containing a given number of correlations, 
and exchange functions joining the correlated particles. In the early calculations, the expansion 
was stopped at low orders | 



94, 95, m 



A great progress in the summation of the cluster 
expansion was achieved by realizing that the cluster terms are conveniently represented by 
cluster diagrams 



97, 



] and by paying attention to the enormous amount of cancelations 
between the different terms of the cluster expansion. A crucial step towards the summation 
was the proper classification of the different cluster terms, i.e. diagrams. Finally, a close set 
of equations for Fermi systems represented by a wave function containing two-body Jastrow 
factor were derived p9| , |100|| . In this way, the uncertainties associated to a low order expansion 
were removed to a large extent. The restriction of the FHNC equations to Bose systems is 
equivalent to the earlier HNC equations used to calculate distribution functions in the context 
of classical theory of liquids ||101|| . 

Rules for constructing the different topological diagrams and derive the FHNC equations 
to sum up the different sets of diagrams have been given several times in the literature, here 
we refer the reader to two recent reviews of Fabrocini and Fantoni where they give in a very 
clear and pedagogical way the derivation of the FHNC equations for both, state independent 
and state depedent correlations, homogeneous and inhomogeneous systems] T7, |^. A review 
which contains an extensive analysis on the different low order cluster expansions and a very 
detailed derivation of the FHNC equations with special attention to the different treatment 
of the exchange terms |^9|, |100|| is the one of Clark in an earlier volume of this series ||92|. In 
that review one can also find results for nuclear matter with simple semi-phenomenological 
potentials that do not require the full operatorial structure of the two-body correlations. 

For brevity, here we only illustrate the main ideas of the HNC summations by considering 
the case of a system of bosons, interacting through a potential that depends only on the 
distance between the particles. This example corresponds to a realistic situation in the context 
of Condensed Matter, i.e. to liquid ^He at zero temperature. In that case the inter-atomic 
potential is much simpler than the NN interaction, in the sense that it depends only of the 
distance between the atoms. However, as in the NN interaction, it has a strong repulsion at 
short distances and weak atraction at medium and large distances. A good representation of 
the interatomic potential is given by the Lennard- Jones potential. 



V{r) = 4e 



12 



(104) 



where e = 10.22 Kelvin (K) ( leV ~ llOOOi^') and a = 2.556 Angstr0m define the energy 
and length scales, respectively. In that case, the non interacting system would be described 
by a wave function with all atoms in the state of zero momentum, i.e. there is a macroscopic 
occupation of a single particle quantum state. This phenomenon is known as Bose-Einstein 
condensation. When the interactions is turned on, the correlations induced by the interaction 
lead to atoms in states with momentum different from zero. However, the liquid ^He still shows 
a condensate fraction, i.e. a ratio between the number of atoms in the zero momentum state 
and the total number of atoms, of around 10%. 
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Figure 13: Different HNC diagrams which appear in the calculation of the distribution 
function, g^ru) 



A convenient trial wave function to describe this system is given by a Jastrow wave function: 

vi>(l, 2, ...A) = F(l, 2, A) = n fir^J), (105) 

i<j 

a product of two body correlations depending only on the distance between the particles. We 
work in the thermodynamical limit, allowing N ^ oo and the volume Q ^ oo but keeping the 
density p = N/Q constant. In this limit, the energy per particle is given by: 



e(p) = I d^rg{r) 



V{r) 



2m. 



-VMn /(r) 



(106) 



where g{r) is the distribution function: 

:i,2) 



9[r) 











2 d^ri...d^rA 



(107) 



which gives the probability to find two particles separated by a distance r, and is normalized 
such that 

p [ d^r[g{r) - 1] = - 1. (108) 



This condition can be used as an accuracy check of the approximations implemented in the cal- 
culation of g{r). In this way, the problem to calculate the energy per particle is now translated 
in into the calculation of the distribution function. Now, the numerator and the denominator 
of Eq. ( |1071 ) are expanded in powers of the correlation interaction, defined as h{r) = [/(r)]^ — 1, 



9[r) 



A{A - 1) I d^r^...d^TA (1 + E Kn^) + j:hh + ... 



J d^ri...d^rA (1 + E h{ri 



(109) 



The terms of the resulting expansion are characterized by integrals containing a given num- 
ber of functions h. Now it is very useful to introduce a diagrammatic notation. This allows to 
classify the integrals according to their diagrammatic representation. The cluster diagrams are 
built employing the following simple rules: Their basic blocks are points and solid lines. Points 
(vertices) represent the coordinate rj of a generic i-particle. Solid points (internal points) imply 
integration over the coordinates times a factor p, circles (external points) refer to particles la- 
beled 1 and 2, and we do not integrate over their coordinates. Solid lines are correlation factors 
and can not be superimposed. Notice that in the diagrams of the denominator, all vertices are 
solid points. 
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In Fig. |T3| we give several examples of diagrams appearing in the numerator of Eq. ( |109| ) 



for which the expressions are given by 

Diagram{a) = P J d^fzh{ri3)h{r2'i) , (HO) 



Diagram{b) = P J d 'r^d r^d rc^h{ri3)h{ri^)h(r2-i)h(r2i)^ (HI) 

and 

Diagram{c) = p^ d^r3d^rih{ri-i)h{rii)h{r3^)h{r2i)h{r23)- (112) 



Cluster diagrams may be linked or unlinked. Unlinked diagrams have at least two parts with 
no common points. They can not be drawn without separating the pen from the paper. The 
linked diagrams can be classified in reducible and irreducible diagrams. The reducible diagrams 
are those which can be factorized in a product of two or more irreducible diagrams such that one 
of them contains the two external points. Diagram (d) of Fig. |13| is a reducible diagram. In the 
infinite systems the reducibility is closely connected to the translationally invariance character 
of the correlations. The irreducible diagrams are then classified into: Nodal (N), Composite 
(X) and the remaining ones called elementary. The nodal diagrams have at least one node, 
that is an internal point such that all ways of going from one 1 to 2 (the two external points) 
should go through it. The composite diagrams are those having two or more (12)-sub-diagrams. 
Where and (ij)-sub-diagram is a part of the diagram connected to the rest only through the 
points i and j. Diagrams a,b,c and d of Fig. O are examples of nodal, composite, elementary 



and reducible, respectively. There are two mathematical operations related to the construction 
of diagrams: the convolution product linked to the construction of nodal diagrams 

(a(rii) I b{ri2)) = p d^ria{rii)b{ri2) (113) 



and the algebraic product {a{rij)h{rij)) for the composite diagrams. 

In evaluating the distribution function Eq. ( p.09| ) one can take advantage of a lot of cancel- 
lations between the different terms. The expansion is linked and in addition in the thermody- 
namical limit, for the specific case of bosons, is irreducible up to terms of the order 1/A^. For 
homogeneous Fermi systems with Jastrow correlations, the expansion is fully irreducible. How- 
ever, for state dependent correlations the cluster expansion is not irreducible, i.e. one should 
consider vertex corrections. In case of finite systems, both for fermions and bosons the cluster 



expansion is not irreducible [93 



The irreducible diagrams are summed up by means of the HNC equation, which allows for 
an iterative process to sum the nodal {N{r)) and composite (X(r)) diagrams once the sum of 
the elementary diagrams {E{r)) is given. The HNC equation is a non-linear integral equation 
relating the function N{r) and X{r) 

{X{ru) I iV(r,2)) = N{n2) - {X{ru) \ X(r,2)), (114) 

i.e., doing the convolution product ( associated to the construction of nodal diagrams) of X(r) 
and N(r) we get all nodal diagrams except the ones corresponding to (X(rij) | X{ri2)). This 
equation can be written also in momentum space as 

iV(A;) = ^iil^ (115) 

where the Fourier transform of a given function a(r) is defined as 

a{r) = p f d\e''' ''a{r). (116) 
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The composite diagrams are given by 



X(r) = f(r)e^ir)+E(r) _ ^ _ ^^^^^ ^^^^^ 

while the distribution function is expressed as 

g{r) = f\r)e^ir)+E(r) = I + x{r) + N{r) (118) 

The HNC integral equation (Eq. (|115|) ) suggests an iterative process to calculate g{r). At 
the first iteration , X{r) = h{r), then the HNC equation is used to construct the first chain 
of diagrams. After a Fourier transformation to r-space, the resulting function N{r) is used 



to define the new X(r) (Eq. (|117D ). The number of diagrams summed in this way grows 
tremendously and in a few iterations one reaches convergence. However, there is a problem 
with the elementary diagrams, entering in the definition of X{r), i.e. the HNC iterative process 
sums up all nodal and composite diagrams , once the sum of the elementary diagrams is known. 
In this sense, the function E{r) is an input for solving the HNC equation. There is no exact 
method to compute this function and approximations are necessary. The simplest option, is 
to take E{r) = 0, known as HNC/0. Actually, this approximation is appropriate for nuclear 
systems, where the density is not very high. 

The Fourier transform of the radial distribution function defines the static structure function 

S{k) = 1 + p f d^r[g{r) - 1] (119) 



which in the case of quantum liquids is experimentally accessible by means of elastic neutron 
scattering against the liquids ||102|| . The condition of Eq. ( p.08| ), implies 5'(0+) = at zero 



temperature. When k oo then S{k) 1. 

A two body correlation appropriate for the Lennard- Jones potential which has been largely 
used in the literature of quantum liquids is the McMillan form 



1 r b<7 



f(r) = e-^i-\ . (120) 

This correlation function together with the corresponding distribution function are shown in 
Fig. |1^ at the ^He saturation density po = 0.365cr~^ in the HNC/0 approximation (i.e. by 



neglecting the elementary diagrams). The correlation function (Eq. ( |120| )) avoids the configu- 
rations where the two particles are close enough to feel the strong repulsion of the potential. 
Finally, the distribution function has a hole at the origin, and the probability to find two 
particles at short distances is strongly depleted. The system minimizes its energy by favoring 
distances between the particles (maximum of g(r)) near to the maximum depth of the potential. 
At the same time, the "holes" in the wave function yield an increment of the kinetic energy, 
which was zero in the uncorrelated system. At the end, there is a delicate balance between 
the potential and the kinetic energy which results in the binding energy of the "^He atoms to 
be —T.nK at the saturation density, while the kinetic energy and the potential energy are 
T ~ 147^^ and V ~ 21K. These arguments of the short range correlations apply equally well 
to the nuclear case, and as we will discuss below, the final binding energy is also a delicate 
balance between potential and kinetic energy. The short range correlations are responsible 
for an increment of the kinetic energy with respect to the kinetic energy of the uncorrelated 
system, which in the case of nuclear matter is given by 3h'^kl/10m, the average kinetic energy 
of the underlying free Fermi sea. Before leaving this example completely it is worth to mention 
that the HNC theory can treat at the same time short- and long-range correlations consistently 
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Figure 14: Two body correlation function (dot-dashed line) and the corresponding g{r) 
calculated in HNC/0 at the experimental saturation density of liquid ^He 



103|, when instead of just using a short-range correlation function, one considers the optimal 



correlation, solution of 
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This optimal correlation has a long range behavior , that translates in the right limit of S{k) 
at low k, S{k) ~ hk/2mvs, being m the mass of the atoms and Vg the speed of sound. This 
behavior reflects the phonon nature of the low energy spectrum of an interacting boson system. 

In the case of Fermi systems, the anti-symmetrization implied by the Pauli principle, in- 
corporates statistical correlations to the cluster expansions, i.e. to the diagrams. In the case 
of homogeneous systems, the statistical correlations are accounted for by the Slater function 
Kkprij), 

3 

Poihj) = = kprij) = — [sin(x) — x cos(x)] (122) 
This Slater function, defines the distribution function for a free Fermi sea, 

g{r) = 1 - -/{kpr). (123) 

To take care of all topological different diagrams, we need to introduce new nodal and 
composite functions, which are related by a system of coupled non-linear integral equations 
(FHNC) ||99|, |100|| . Once the FHNC equations are solved one can easily calculate the energy. 



The presence of the symmetrizer in the correlation operator and the non-commutativity of the 
operators involved in the correlations makes it very difficult to sum up all the contributions 
in the case of state dependent correlations. Actually, no scheme has been devised so far, 
which keeps track of all possible terms. A complete FHNC summation is possible only for the 
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Jastrow component. However, partial classes of diagrams containing operatorial correlations 
may exactly be summed by the Single Operator Chain approximation (FHNC/SOC)|TD^. The 
FHNC/SOC integral equations sum all the nodal diagrams containing only one operatorial 
correlation per internal side, in addition to all the Jastrow correlated clusters. By computing 
its leading corrections, the estimated accuracy of the FHNC/SOC approximation has been set 
to ~ 1 MeV/A for nuclear matter at saturation density ||105|| . The FHNC/SOC has been applied 
to study the equation of state of nuclear and neutron matter [jl5[, including three body terms 
in the hamiltonian. It is very difficult to apply the variational formalism with state dependent 
correlations to study asymmetric nuclear matter with an arbitrary asymmetry x = pp/ p and all 
the calculations, which are then later on used to study /?— stable matter, are based in quadratic 
parameterizations of the correlation energy 



E{p, x) = Tf{p, x) + Voip) + (1 - 2x) V2(p) 



(124) 



where 



Tf(p, x) 



X 



5/3 



+ (1 



X 



,5/3 



(125) 



is the kinetic energy of the asymmetric Fermi sea and Vq and V2 are related to the correlation 
energy of neutron and nuclear matter. In this parameterization, the symmetry energy of nuclear 
matter is 

ld'^E{p,x) 5, 

U=l/2 = 



Es{p) 



-Tf{pM2) + V2{p). 



(126) 



\ dx '^-^'^ 9" 

The extension of the FHNC/SOC to finite nuclei is an important issue in order to have 
a systematic study with the same basic tools for nuclei and nuclear matter. After the first 
efforts based on low order cluster expansions ||106| , |107|| , the FHNC formalism was extended to 
doubly closed shell nuclei, both in Is and jj coupling schemes, using Jastrow correlated wave 
functions ||108| , |109| , |110|| and semi-realistic interactions ||111|| . Very recently, the FHNC/SOC 
has been generalized to deal with realistic potentials and state dependent correlations, with 
non-central components for double closed shell nuclei as ^^O and ^°Ca ||112| , |113|| . The most 
complete calculations ||113|| deal with the Argonne Vi4,m supplemented with the Urbana VII 



model of the three-body interaction |S^. The two-body correlation function includes operators 
up to tensor components in Eq. (|98D and the single particle wave functions were taken from a 
Woods-Saxon mean field potential. The results are rather promising and one can conclude that 
now a days the calculations in finite nuclei with realistic interactions have achieved the same 
degree of accuracy than in nuclear matter. 

As we have already mentioned, a full minimization of the variational energy respect to 
arbitrary variations of /^^■'(r) is impossible. The correlations, in the most sophisticated cal- 
culations, both in nuclear matter and finite nuclei, are generated by solving a set of coupled 
Euler-Lagrange equations obtained by minimizing the second order cluster expansion of the 
energy. The differential equations are solved with the boundary conditions: 



and 



f{p>i)^r > d(P>i)) = 0, 



(127) 



:i28) 



where S^"' are the healing distances in each correlation channel. To ensure the continuity of 
the first derivative of the correlations 



dfip) 
Or 



=dip)'- 







(129) 
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one introduces Lagrange multipliers A^^^^ which are varied to satisfy the condition of the above 
equation. 

Some additional variational parameters are introduced because the second order minimiza- 
tion is actually performed by using a "quenched" potential 



Y,a^P\^P\rij)d^P\ (130) 



These quenching parameters a'^^^ are interpreted as representing medium effects. In this way 
the two body correlations contain several variational parameters 

Fv= E f^'\r,r.d^''\a^''^)Of^ (131) 

p=l,8 

which will be determined by energy minimization. The number of parameters is reduced in a 
drastic way by taking all the healing distances equal S-^'^ = dc, except those associated to the 
tensor channels dt- The quenching factors are typically assumed to be 1 for the scalar channels 
and those containing L^, while all other channels are quenched by the same amount a^^^ = a 

151, im. 



A main drawback of a variational calculation is that even if the expectation value of the 
hamiltonian is evaluated exactly for a given trial wave function, one has only un upper-bound 
on the exact energy and an approximation to the exact many body wave function. The way 
to proceed further is provided by the Correlated Basis Function (CBF) Method, which was 
introduced nearly forty years ago by Feenberg and his collaborators ||91], |115|| . Several people 
have contributed later on to make the method efficient and competitive, pushing its applications 



to very different physical systems p7| , p.l6|| . The CBF method provides the general frame to 
systematically improve on a trial wave function by introducing a complete set of correlated basis 
functions. The idea is to perform a perturbation theory built on such correlated basis, with the 
hope that the correlated basis states are close enough to the eigenstates of the hamiltonian that 
a lower order in the perturbation theory will be sufficient to take care of these small differences. 

The set of correlated basis states \l/„, is built by applying the many-body correlation oper- 
ator, F{1,2, A), usually determined in the variational calculation, to model basis functions, 

(t>n, 

^„(1, 2, A) = F(l, 2, A)M^, 2, A). (132) 

The zero order of the CBF theory, i.e. taking only one basis function, reduces to the 
variational procedure. Therefore it is crucial to have a good variational description of the 
ground state in order to have a fast convergence of the perturbative series. The model function 
00 is also used to define the Fermi level and as before particle states above and hole states 
below this Fermi level. The model states (pn refer then to some particle-hole configuration with 
respect to 0o- 

In order to perform perturbation theory we need to decompose the hamiltonian in an un- 
perturbed hamiltonian Hq and a perturbation Hj, 

H = Ho + Hi. (133) 

Here they are defined in terms of matrix elements of H in the correlated basis | i). The diagonal 
matrix elements 



define the unperturbed hamiltonian: 

{z\Ho\ j) = 5,,,El (135) 
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while Hj is given by 

{i\Hj\ j) = (1 - s,,) '"i:; (136) 





H 


j) 




j) 



By construction the diagonal matrix elements of Hj are zero and therefore there are no first 
order terms in the perturbation expansion. 

A main characteristic of the CBF theory is the non-orthogonality of the basis, 



(^1 


J) 






1 j)y^' 



= ./.i;/ I ^ ^^^r (137) 



Therefore, it is necessary to use non-orthogonal perturbation theory ||117|| . The perturbative 
correction to Eq, in a non-orthogonal basis is given by ||117|| 

^ (J5„-b;)(b„-b;)+-- 

Due to the non-orthogonality of the basis, the interaction Wij{0) contains the overlap matrix 
element Nij, 

WijiO) = Hij - EoNi, (139) 
and carries an additional energy dependence, besides the one coming from the denominators. 



To eliminate the unknown Eq from the perturbative series, one expands Eq. ( p.38| ) around 
Eq = Eq. In this way, 

W,,{Q) = {H,, - EIN,^) - AEqN,^ = W^^iO) - AEqN.j. (140) 

Notice that, VFjj(O) depends on the variational energy Eq. The expansion for the energy 
denominators gives 



Eq-e; e-q-e; + aEq E-Q-E;f^Q\ e-q~e; 

Actually, the analysis of the series is easier by expressing, 

W:,{Q) = W:^{i) + {E^-El)N,,, (142) 
with W^jii) = Hij — E^Nij. The series can then be written in terms of Wij{i) 

^Eq = E(W)-^^oA^op)f^S^-ivJ 

. ZJW,A0) - A^o^.) - N„) - N,) . ... ,143, 

Due to the correlations introduced in the wave functions, the calculation of matrix elements 
is much more involved than in normal perturbation theory. The CB matrix elements are 
constructed from the matrix elements of F~^F and F^HF on Fermi gas states, characterized by 
the number of particles and holes. The matrix elements are calculated by cluster expansions. 
After expanding the energy denominators of the series (Eq. (|143|) ), a careful analysis of the 
cluster expansion of the matrix elements and of the energy dependence of the different terms 
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leads to the conclusion that AEq is given by the sum of all linked Goldstone type diagrams 
(see section built up with W^j{i) and Nij interaction boxes. The interaction boxes (or 
equivalently the matrix elements) are the sums of only linked cluster diagrams. Therefore, the 
CBF perturbation series is well behaved in the thermodynamic limit ||118 . 



In dealing with CBF perturbation theory we are confronted with two types of diagrams, the 
Goldstone type diagrams representing the terms of the series and the diagrams corresponding 
to the cluster expansion of the matrix elements, which are summed up by means of FHNC 
techniques. 

Perturbative CBF calculations for the binding energy of nuclear matter have been performed 
up to second order. In this case, the perturbative correction to the ground state energy is 
expressed as: 

and has been calculated within the set of 2p—2h correlated basis states. In the thermodynamical 
limit, the energy denominators can be expressed in terms of the variational estimate of the single 
particle energies,e''(p), 

Eo" - E^ip,, p2, hi, h2) = e^ih) + e-{h2) - e^p^) - e^p^). (145) 

As expected, the second order correction to the ground state energy, increases the binding. 
In the first calculations, the correlation operator was taken of Jastrow type and in this case 
the perturbative second order correction was large ||119|| . On the other hand, the corrections 
are much smaller, when the correlation with the full operatorial structure is used. Therefore, a 
significant improvement in the convergence of the series is expected. Typical values of second 
order corrections to the binding energy at saturation density of nuclear matter range from ~ — 1 
to ~ —3 MeV. A very recent estimation for the Argonne Vi^ and the Urbana IX model for the 
three-body force, is —1.89 MeV at saturation density [|T6[. 



The CBF perturbation theory has also been implemented to study other observables, besides 
the hamiltonian, like the self-energy, the momentum distribution, the single particle spectral 
function, the different nuclear response functions, etc. In the previous section ^]3| we have 
considered the self-energy of a nucleon as the key quantity to calculate the single-particle 
Green's function. Here one can calculate the imaginary part and split it similar to (|79[) and 
([5D| ) into two pieces, the correlation term, i.e. coupling to 2hlp correlated configurations and 
the polarization term corresponding to the coupling to 2plh configurations. The imaginary 
part of T?^^'P{k,uj) for a particle state, in second order perturbation theory is defined as 



ImL^^^^ik, uo < e^ikp)) = I ^ I (k I i7 - i?^' I Pi; hi, h^) ? 6{uj + e^(pi) - e^{h,) - e''{h,)) 

(146) 

while 

ImT?P^\k, u > e^kp)) = I ^ I (k I if - I pi, P2; h^) \' 5{u; + e^h^) - e^(pi) - e^p^)^ 

(147) 

Similar expressions hold for the hole states. Using dispersion relations one can calculate the 
corrections to the variational estimate of the single particle potential 

AeW^.W-."w4££^<^. (148) 
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In the case of the single particle momentum distributions, the second order perturbative 
correction to n^{k) has been obtained by calculating the expectation value of the occupation 
number operator aj^-Ok on the perturbed ground state | Oper) containing two particle two-hole 
correlated states: 

I Oper) =1 0^) + ^ H «(Pi5 P2; hi, ha) I P1P2; hihs), (149) 

where 

( u u \ (OJ if-E^ I Pi,P2;hi,h2) , - 

tt Pi, P2; hi, h2 = , . . — ^ ^r^' ^^"^ 

and keeping the terms of order ||120| . 



However, there are still some problems with orthogonality corrections when the series is 
truncated at a given order. Moreover, it is not straightforward to extract the orthogonalized 
eigenvectors in a non-orthogonal CBF theory. Therefore, even if the applications of the non- 
orthogonal perturbation theory are rather successful, it is desirable to have also a method 
available in a orthogonal scheme. In this way, the Nij boxes would be eliminated from the 
diagrams and the parallelism with normal perturbation theory becomes more transparent. To 
this end, one can try a Schmidt orthogonalization process, however it happens that the vari- 
ational estimates of the hamiltonian on the Schmid orthogonalized states are worse than the 
ones calculated with the non-orthogonal basis. Very recently, a new orthogonalization proce- 
dure has been devised ||121|| , in which the diagonal matrix elements of the hamiltonian and of 



the identity operator are preserved. It is a two step process which combines the Schmid orthog- 



onalization with the Lowdin transformation [ p.22|| . Standard perturbation theory may be used 



in this new set of states. Although both schemes should be equivalent, in practice orthogonal 



CBF theory is more efficient particularly in calculating the spectral functions ||123| , p.24|| or the 
nuclear response functions ||125| , |126| , |127| . 



The one-body Green's function has also been studied in the framework of CBF. Instead of 
computing the Green's function from a Dyson equation, as in the perturbative scheme discussed 
in the previous section, one directly calculates the single particle spectral functions. The 
starting point is the variational estimate, which is equivalent to approximate the true ground 
state by the variational ground state wave function and take 1 particle and 2 particle-1 hole 
basis correlated states for A-l-1 particles in the case of Sp{k,uj) and the 1 hole and 2 hole -1 
particle correlated states for A-1 particles for Sh{k,uj). The projection into 1 hole (1 particle) 
state is the dominant one and is given by 

S%,){k,uj) =\ ^t^'^ r Siuj - e\k))e [±{kF - k)] , (151) 

where the upper (lower) sign in the argument of the B-function applies to the hole (particle) 
spectral function, e^{k) is the variational single particle energy and the overlap matrix element 
cp'^^ is defined by 

^'^={0,\ai\h). (152) 

Where | h) is the basis correlated state of A-1 particle build on a Slater determinant where the 
state with momentum k is missing. In a similar way, 

= {Oy I Ok I p). (153) 

where in this case the state | p) represents the correlated state of A-l-1 particles build on a 
Slater determinant that besides having all the states occupied up to the Fermi level it has an 
additional particle on the state k above the Fermi level. Actually these type of intermediate 
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states are the only non-vanishing contributions in any uncorrelated ground-state. In such a 
case, the spectral functions Sunc,h(p){k,uj) are just ^-functions with strength unity: 

SuncMp) = - 7^)0 [Hkp - A;)] . (154) 

As correlations are present in the CBF states also other intermediate states contribute to the 
zero order approximation to the spectral function. The contribution from 2h-lp (2p-lh) corre- 
lated states is spread out in energy and can be interpreted as a background contribution to be 
added to the quasi-particle part (5-function) of the spectral function. This spreading produces 
a considerable quench of the quasi-particle peak. On top of this variational estimate one can 
build perturbation corrections, which take into account the admixture of more complicated 
configurations in the states considered in the definition of the spectral function. For instance, 
the 2p-2h admixture in the ground state, and the 2h-lp (2p-lh) in the one-hole (particle) have 
been studied ||123| , |124|| . As a consequence of these corrections, the 5-peak of the quasiparticle. 



acquires a width (except for k = kp), which is related to the imaginary part of the self-energy 
at the quasiparticle-energy. For k = kp, the on-shell imaginary part of the self-energy is zero 
and the quasiparticle peak is a 5-function which strength {Zkp) which defines the discontinuity 
of the momentum distribution at k = kp. 

We can conclude that CBF has reach its maturity and systematic studies of different observ- 
ables, for both nuclear matter and nuclei are available. The starting wave function is realistic 
enough that the perturbative corrections to the different quantities can be safely calculated up 
to second order. The formal properties of the perturbative series have been established and the 
efforts now can be addressed to improve the evaluation of the matrix elements or to study also 
higher orders in the perturbative expansion. 



2.5 Variational and Quantum Monte Carlo Techniques 

There are various possibilities to apply the techniques of Monte Carlo (MC) sampling in cal- 
culations of quantum many-body systems. For example one can use MC based algorithms to 
sample the multidimensional integral occurring in the evaluation of the expectation value of 
the hamiltonian, or any other operator, on a given trial function (Eq. (p2|)). In fact for simple 
interactions and central correlations the procedure is well established and it has been applied 
both for bosons and fermions, for systems with finite and infinite particle numbers. This type 
of calculations are known as Variational Monte Carlo (VMC) methods. However, there are also 
methods based on stochastic algorithms, which aim at an exact integration of the many-body 
Schrodinger equation. Those are generically known as Quantum Monte Carlo (QMC) meth- 
ods. VMC and QMC calculations in systems of liquid helium, both "^He and ^He, which follow 
Bose and Fermi statistics respectively, have motivated an enormous progress in Monte Carlo 
Methods. However, the nuclear many-body hamiltonian is much more complicated and requires 
to introduce strong state dependence correlations in the nuclear wave function. Therefore the 
progress in accurate nuclear ground state calculations based on MC techniques has been rather 
slow. It is only recently, that the drastic improvement of the computational resources, i.e. 
parallel computers, allows to increase the number of nucleons involved in the calculations. 

Once the trial wave function is chosen, the variational Monte Carlo method provides an 
exact - at least in the statistical sense - evaluation of the expectation value, which is free 
of the approximations involved in the integral equation methods discussed in the previous 
section. However, we should not consider these methods as competing ones but rather as a 
complementary ones. For instance, it is clear that Monte Carlo methods have difficulties in 
considering the long-range behavior of the distribution function of an infinite homogeneous 
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system because of the finite size of the simulating box (see below). On the other hand one 
can very well determine this long-range behavior in the context of HNC methods. In the 
case of finite nuclei, Monte Carlo methods have enormous problems in increasing the number 
of nucleons of the nucleus under consideration. On the other hand recent calculations with 
FHNC techniques can easily be applied to nuclei with A=40. For the systems, where both 
methods can be applied, they can help each other in looking for the appropriate correlation, in 
understanding the physical requirements of the wave function and in evaluating the accuracy 
of the approximations involved in the FHNC method. 

Recent reviews on Quantum Monte Carlo methods can be found in Refs. |P!H|, |128| , |129|| . 
Another review-like book with presentation of the methodology and physical applications is 



130|| . Two introductions into this field and rather comprehensive reviews by R. Guardiola on 



both the Variational and the Quantum Monte Carlo method have been published in Refs. [131 



132|] . More oriented to nuclear physics is the paper of Carlson and Wiringa on variational Monte 



Carlo in finite nuclei |133|. A recent overview including also Green Function Monte Carlo for 



nuclear systems is due to S. Pieper ||134|| . Finally, the generalization to finite temperature. 



focusing mainly to applications in condensed matter physics has recently been reviewed by D. 
Ceperley 

Here we will mainly discuss VMC calculations. In order to explain the main ideas we will, as 



in the previous sub-section, use a simple system of bosons, obeying the hamiltonian of Eq. (93) 



with the scalar interaction of Eq. (|104|). First we consider a droplet of A ^He atoms. A trial 



function satisfying the symmetry requirements is given by, 

A A 

^(n, r2, r^) = n /(^u) n (155) 

i<j i=l 

i.e., the product of a Jastrow factor times a model wave function that assigns the same single 
particle wave function to all the atoms. In this case, for spinless particles and local interactions, 
the expectation value of the hamiltonian can easily be written in a form appropriate for Monte 
Carlo calculations. Let R represent the set of the 3^4 spacial coordinates of the atoms in the 
droplet. We define a local energy -El(R) by 

^l(R) = ^^^(R-) (156) 
and introduce a probability distribution function 

which by definition is normalized to 1. Notice also that for local potentials, the contribution 
of the potential to the local energy is just the value of the potential for the given configuration 
R. The energy is given by 

The Monte Carlo way of evaluating this integral is to consider p(R) as the importance sampling 
distribution and El(R) as the function to be sampled. The essence of the Monte Carlo method 
is to generate a set of statistically independent configurations (Rj) which are distributed ac- 
cording to the probability distribution p(R). Once this set of points (Rj) is obtained, the 
expectation value may be calculated through 

1 ^ 

E. = ^m j^T^Eim. (159) 

1=1 
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Actually, this way to organize the integral of Eq. ( p,58| ) is a special example of the more general 
expression 



F - ' ^^^^ " cifin^ 



where iy(R) is a probability distribution function. When iy(R) = \E'"I'(R)^(R) we recover 
Eq. (|158|) . Once we have a set of configurations in the 3 A dimensional space proportional to 
Vr(R) then the calculation of or the expectation value of any other operator O 

I O I ^) 

is straightforward 

^t(R)0^(R) 

O -- ^ (1Q2) 

^ W(R) 

where the sum runs over all the configurations (Rj). 

In order to get an estimate for the accuracy of the MC evaluation of the integrals, the 
samplings are grouped in blocks, each block with a sufficiently large number of points (Rj j = 
l,nmov), nmov being the number of configurations in each block. The central limit theorem 
guarantees that the average values of the energy obtained in each block i, 

1 nmov 

K = E M^^) (163) 

nmov , 

are distributed as a Gaussian centered around the true average value Ey 

-| nblock 

where nblock is the number of blocks and the statistical error is given by 

-1 / 1 nblock \ 1/2 

The configurations (Rj) following the importance sampling function can be generated through 
the so-called Metropolis algorithm ||135|| , which yields a sequence of random numbers, a Markov 



chain, distributed according to the required distribution probability. The way to generate a 
sequence of configurations (Rj) by means of the Metropolis algorithm is rather simple and may 
be described as follows: Assume that we are at a given configuration Roid, and the value of 
the wave function for this configuration is Woid = ^(Ro«d) then one tries to generate a new 
configuration by defining 3A new coordinates 

Rne«, = Rozd + (RAND() - 0.5) * step. (166) 

where RAND() produces random numbers uniformly distributed in [0, 1] and step defines the 
scale of the Monte Carlo moves. The value of the wave function at Rnew is Wmw 

Now we perform the Metropolis question to see if we are going to accept this new configu- 
ration Rnew in our Markov-chain: If the probability of the new configuration p{Rnew) is larger 
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than the older, i.e. if w"^^^ > w^i^, then Unew is accepted. If the probabihty is smaller, then 
Unew is accepted with probability w^^^/wli^. Both cases are included in the following condition: 

-^KwlJwli, < RANDO) (167) 

then Unew is accepted and is the starting configuration for the next move. If the condition 
( |167| ) is not fulfilled, the move is rejected and we try to generate another configuration from 



Hold- One keeps track of the number of acceptances, and a efficient calculation should keep the 
percentage of acceptances between 40% and 60%. Also one can allow for several moves between 
the different blocks in order to reduce the correlations between the different calculations of the 
energy. It is also convenient to perform several moves at the beginning of the process, in order to 
reach the region where the Metropolis algorithm works properly. In the Monte Carlo language 
this is known as thermalization. 

In the case of infinite homogenous systems one can of course not perform a simulation 
with an infinite number of particles. The usual way of dealing with these systems consists in 
representing the system in a simulation cell with periodic boundary conditions. If the size of 
the cells is L^, the density of the system determines the number of particles A that we should 
consider for each cell, p = A/Lr'. Metropolis moves are always recast into the main simulation 
cell, but in implementing the algorithm one must consider the particles in the periodic images 
as well. At a given density, the number of particles which can be treated in a calculation 
defines a clear limitation to the size of the simulation cell. Typical number of particles in the 
simulation box range between 50 and 100. 

There are no essential differences in the case of fermion systems with respect to what we have 
discussed for boson systems so far. There are, however, some technical problems, associated 
with the calculation of determinants (the model wave function should be antisymmetric) and the 
manipulation of spin and isospin degrees of freedom. When the interaction and the correlations 
do not depend on the discrete (spin and isospin) degrees of freedom, the trial wave function can 
be considered as the product of a Jastrow factor and as many Slater determinants as fermion 
species, i.e. one corresponding to spin up and the other to spin down. Once a configuration 
is defined, each determinant is a number and the acceptance criterium for a Metropolis move 
contains the squared quotient of two determinants times the contributions of the correlation 
factors. The number of operations involved in the computation of a determinant is proportional 
to A^ as a consequence, moving A particles involves A^ operations. There is an algorithm, 
based on the calculation of a determinant by using the matrix of cofactors which reduces the 
number of operations to A^, which is of course important for calculation with larger number of 
fermions||136| . 



In the variational method one aims to a minimization of the energy. This implies that 
one should perform VMC calculations for several variational parameters. This usually requires 
huge amounts of computing time. Besides, the statistical fiuctuations of the results may give 
an incorrect localization of the minimum. To avoid this problem, a method called reweighting 
of configurations has been introduced. This method tries to use configurations corresponding 
to the same importance sampling function to calculate the energy for wave functions which 



have small changes in their variational parameters | 136 |. 

In the nuclear case, the strong state dependence of the interaction and the correlation, 
produce several complications that limit the application of the method to a small number of 
nucleons. Actually all the applications with realistic potentials refer to light nuclei and only 
very recently attempts have been made to study homogeneous systems [|137| . 



For a finite number of nucleons, the trial wave function is a symmetrized product of non- 
commuting two-body operatorial correlations F2{i,j) (Eq. (|98[) ) times a Slater determinant of 
single particle functions (Eq. (|9TD). The trial wave function must be translationally invariant. 
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and the correlation function fulfils this condition automatically . Also the implementation of 
the translational invariant to the model wave function is rather straightforward in the Monte 
Carlo method. For instance, if the model function has some dependence on the single particle 
coordinates rj , one just substitutes that set of coordinates by fj with f j = — Rem and 
R-cm = lMZ]i"j- The fact that the trial wave function describes a localized bound state, is 
either refiected in the correlation operator 

Fij{rij) ^0 , ^ oo (168) 

or in the model wave function, by requiring that the model wave function tends to zero when 
one of the particles is moved far away. 

For very light nuclei one can take the model function (0) as a pure spin-isospin function 
with no spatial dependence, for instance for triton and the ^He nucleus the following functions 
have been used ||133| , |138|| 

\(l)CH))=A\in]nlp) (169) 

and 

m^He)) =A\in^ n Ip^l n) (170) 

where A is the antisymmetrizer and the up- and down-arrows indicate the spin of the different 
nucleons. In this case the localization is impose on the correlation functions. One of the 
problems is that the size of the spin-isospin space of A nucleons increases very rapidly. There 
are Ng = 2^ possible spin states ranging from all particles having their spins pointing down 
to all having their spins pointing up. If one fixes the charge of the system, i.e. the number of 
protons, then there are 

A\ 

^'^^^^ = Z\{A-Z)\ ^^^^^ 
states in the isospin space with the same total charge. The total dimension is therefore Ndim = 
Ns X Nchar For instance in the case of the "^He, A^^ = 12 and Nchar = 6 with Ndim = 96. 
While for the =^H, A", = 8 and Nchar = 3 while N^im = 24. However for the ^^O, A^, = 65536 
and Nchar = 12870 while Ndim = 843448320. Each correlation operator Fij is a very sparse 
matrix in the spin-isospin space. In evaluating the expectation value, this matrix acts on a 
vector in spin-isospin space to produce another vector. For instance the spin-isospin vector 
corresponding to the model wave function of ^He (Eq. ( p.70| )) has 24 non zero components, as a 
result of the antisymmetrization. When the correlations act on those initial vectors other states 
are generated. One can try a straightforward generalization of the VMC methods discussed 
above by writing 

M 

^ = E^P (172) 
p=i 

where corresponds to the term in in which the A(A-l)/2 correlation operators Fij operate 
in a specific order labeled by the index p. The energy expectation value is given by 

where the spin isospin degrees of freedom are treated explicitly. However, to sample the coordi- 
nates using the Metropolis algorithm would still be not practical, because it would be necessary 
to take into account all the (^(^4 — l)/2)^ possible terms, related to the ordering of the corre- 
lation operators in each side of the expectation value. Instead, one can sample randomly the 
coordinates R and the order labels p and q using for the importance sampling function 

W„{K) =1 i?e[vl/t(R)vl/^(R) I (174) 
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where the absolute magnitude is required in order to be sure that the weight function is never 
less than zero. The statistical error per sample in configuration space is increased by sampling 
the order of the operators. However, the time saved by evaluating each time only one term 
in the sum (Eq. ( |173| )) increases the efficiency of the calculation. The way to proceed is the 
following: From an initial configuration X = (R, p&g) a trial configuration X' = (R', p'hq') is 
generated. The Metropolis condition (|167|) is now applied to the sampling functions Wpq{R) 
and Wp'q'i^') to see if X' is accepted or rejected. 

Special codifications of the action of the operators on the spin-isospin vectors, based on 
binary representations of those vectors are useful. However, calculations performing the full 
spin-isospin summations, as the ones just described, are feasible with present computers only 
up to A=8. 

To study heavier nuclei one needs to sample also the spin-isospin states in a random walk 



which takes place in a combined spin, isospin, and coordinate space [|139|| . The main problem 
in sampling over spin-isospin states is the increase in the statistical error if the sampling is 
done blindly. In order to reduce the variance, the spin-isospin states must be treated in a 
manner similar to that commonly used for the spatial coordinates. In order to devise the 
proper strategy, we write the expectation value of the hamiltonian as 

/rfREa^i(R)^a(R) ' 

in this case a and h represent spin-isospin states. The basis of spin-isospin states is just the 
one discussed above, i.e. the states where each nucleon has a definite third component of spin 
and isospin. The Monte Carlo strategy consists in sampling values of R and a, while explicitly 
summing over all spin-isospin states 6, for each choice of R and a. The full summation over 6, 
although it is a lengthy calculation, reduces the statistical error of the calculation significantly, 
since for the exact wave function 

Y,H,,^, = E^, (176) 

b 

for each configuration (R, a) . 

The generalization of the standard Metropolis algorithm in order to sample R and a in the 
combined coordinate-spin-isospin space, must allow not only to move the particles in coordinate 
space but also for changes in the orientation of the spins and isospins of the nucleons. The way 
to do that is by fiipping the third component of the spin and/or by exchanging the spin-isospin 
of two nucleons, chosen in a random way. Of course these are not the only possibilities but 
these movements are easy to implement and enough to sample the full spin-isospin space. As 
before, the proposed move is accepted or rejected according to the ratio of the squares of the 
wave function evaluated at (R, a) and (R', a'), as in standard variational Monte Carlo. In this 
way the Monte Carlo method has been applied to ^^O, to calculate the energy, form factor and 
charge distribution, using the Reid V6 interaction, i.e. including up to the tensor channels in 
the general expression of the force (Eqs. ( p^j95| )). 

Another alternative is to perform a cluster expansion for the non-central correlations while 
the central correlations and the antisymmetry are treated to all orders. The different terms are 
exactly evaluated with Monte Carlo techniques. This procedure has allowed to study ^^O with 
the full Argonne Vi^^ two-nucleon interaction, supplemented with the Urbana VII three-nucleon 
potentials, using a trial wave function which contains pair- and triplet-correlation operators. 
Not only the binding energy but also other quantities like density distribution, momentum 



distributions, charge form factors and response functions have been studied |T^, 141 . 

To go beyond VMC one needs to invent an algorithm to generate the exact ground state 
wave function. These new algorithms are generally called as Quantum Monte Carlo and they 



48 



try to integrate the many-body Schrodinger equation by means of stochastic procedures. The 
simplest of these methods is the so called Diffusion Monte Carlo (DMC) which approximates 
the many-body time dependent Green's function by calculating the time evolution operator 
in small time steps. The starting point is to consider the many-body Schrodinger equation in 
imaginary time, 



dt 

The formal solution of Eq. (|177|) is 



(if-E)^(R,t) (177) 



^(R,t) =exp(-[i/-^]t)^(R,0) (178) 

where the propagator exp [—{H — E)t] is called Green's function, and E is a convenient energy 
shift. The imaginary time evolution of an arbitrary starting state \E'(R, 0), once expanded in 
the basis of stationary states (0j(R)) of the hamiltonian, is given by 

^(R,t) = 5]exp(-[E,-E]t)Q0i(R). (179) 

i 

The amplitudes of the different states change with time, increasing or decreasing depending on 
the sign of {Ei — E). Independently of the value of the energy shift E, the most important 
amplitude after a long time will be the one corresponding to the state with the lowest energy 
Eq. In other words the DMC method projects out the ground state wave function from the 
trial wave function using 

^0 = hm exp [-{H - ^o)i]^T , (180) 

t — ^oo 

provided that the overlap between the starting trial function and the true ground state is 
different than zero (cq 7^ 0). The eigenvalue Eq is exactly calculated during the process. The 
starting trial wave function is represented by a set of random vectors or walkers (Ri, R2, ....Riv^) 
distributed with probability proportional to the trial wave function, in such a form that the 
time evolution of the wave function is represented by the evolution of the set of walkers. Notice 
that the quantity interpreted as a probability distribution function is the wave function of 
the system, and not its squared as in VMC. Along the iteration process, going from a wave 
function to another wave function, means to change the random numbers representing the 
wave function. Each walker generates a number of descendants (none, one or more than one) 
which will finally accommodate to the exact wave function. The method is iterative and the 
ground state is asymptotically approached. In its simplest form, the DMC is only appropriate 
to describe the ground state of boson systems. In this case the wave function is positive and 
can be used as a probability distribution function. In the case of fermions, the ground state 
wave function is not positive definite, and one should be careful to interpret the wave function 
as a probability distribution function, giving rise to the so called " sign problem" . Importance 
sampling combined with the so called fixed node approximation allows to build a method also 
in this case. 

The actual computation of the time evolution is done in small time steps r, by writing the 
Green's function as a product over short time intervals 

n 

exp (- [H-E]t) = l[ exp (- [H - E]t), (181) 

i=l 

where r = t/n. These short time propagators are then approximated as the product of free 
particle propagators and a propagator involving the potential, 

exp {-Ht) ^ exp {-Vt) exp {-Kt) (182) 
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where V is the interaction and K the kinetic energy operator. Their coordinate representation 
is rather simple, 



GK{n',-R,t) = (R' I exp~Kt I R) = -^^^^^^^^ exp (-(R' - R)V4Dt) . (183) 

Where D = /2m is called the diffusion constant. If the potential is local the Green's function 
related to the potential is, 

Gy(R', R, t) = (R' I exp {-Vt) \ R) = exp (-1/(R)t)5(R - R'). (184) 

The coordinate representation of the total Green function for short times including the 
energy shift is given by 

^(R', R, r) = GkCR', R, T)exp[E - V{R)] + 0{t^). (185) 

Using the above expressions, the propagation in coordinate space reads 

(R I ^(r + 6t)) = J dR'{R \ G{6t) \ R')(R' | ^(r)). (186) 

This integral equation is solved by using Monte Carlo methods to sample the free particle 
propagator and including branching to take the potential into account, thus allowing for the 
process of replication of the walkers. The equation is iterated until convergence. In conclusion 
the scheme is well defined once one has found an appropriate approximation for the short time 
Green's function, and a collection of walkers describing the starting state has been chosen. 
These walkers are evolved in time, obtaining a new collection of walkers. If the number of 
time steps is large enough the final walkers will represent the ground state wave function. An 
important improvement is obtained by introducing an importance sampling function to increase 
the statistical accuracy of the method. 

The method is well established for systems of bosons interacting through central scalar 
interactions. For fermions one must deal with the so called sign problem, on which a lot 



of progress has been achieved in the last years, see |p.44| , p.43|| . Realistic applications are in 
the context of quantum liquids to study ^He and ^He. In the nuclear case, using realistic 
interactions there are only applications for a few number of nucleons. Recent calculations are 



able to arrive up to 7 nucleons, calculating the ground state and low-lying excited states |p.42 
Also new algorithms are being designed to treat the nuclear or neutron matter case with realistic 
interactions [|137| , |145 . 



3 Effects of Correlations derived from Realistic Interac- 
tions 

3.1 Models for the NN Interaction 

In our days there is a general agreement between physicists working on this field, that quantum 
chromo dynamics (QCD) provides the basic theory of the strong interaction. Therefore also 
the roots of the strong interaction between two nucleons must be hidden in QCD. For nuclear 
structure calculations, however, one needs to determine the NN interaction at low energies and 
momenta, a region in which one cannot treat QCD by means of perturbation theory. On the 
other hand, the system of two interacting nucleons is by far too complicate to be treated by 
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Figure 15: Nucleoli- Nucleoli interaction in the quark- and meson-exchange picture. The 
lines in the left part of the figure refer to propagating quarks, whereas those in the right 
part refer to nucleons and mesons. 



means of lattice QCD calculations. Therefore one has to consider phenomenological models for 
the NN interaction. 

Attempts have been made to develop models for the NN interaction, which consider the QCD 
degrees of freedom, quarks and gluons explicitly. As an example we mention the non-relativistic 
constituent quark models, which are very successful in describing the properties of baryons[ p.46| ]. 
Based on such constituent quark models, so-called quark cluster models have been developed for 
the NN interaction [0, p.47| , |148|| . These models include the effects of one-gluon-exchange terms 
and account for the exchange of mesons between the quarks. This exchange of mesons within 



the quark model is schematically displayed in the left part of Fig. |1^. The exchange of a quark- 
antiquark pair is substituted by the exchange of mesons of various kinds. This interpretation 
is in line with the arguments of t'Hooft and Witten [|149| , |15(]|] . They demonstrated that in the 
low-energy regime, the relevant degrees of freedom of QCD should be well described in terms of 
a meson theory. This means that also the quark-antiquark exchange processes displayed in that 
figure should be dominated by the exchange of the collective quark-antiquark modes, i.e. the 
exchange of mesons. Such constituent quark cluster models are very successful in describing 
the main features of the baryon-baryon interaction, however, the accuracy of the fits of these 
interactions to the empirical NN scattering phase shifts is not sufficient to use such interactions 
based on a quark model for nuclear structure calculations. 

If one considers the meson degrees of freedom as the most important degrees of freedom 
for QCD at low energies, it is quite natural to describe the NN interaction in terms of mesons, 
which are exchanged between nucleons. The internal structure of the nucleons can then be 
parameterized in terms of meson-nucleon form-factors. Detailed discussions of these meson 
exchange models have been presented, e.g. in B, p.51|, p.52|, [153||. The basic idea of this meson 



exchange model dates back to 1935, when Yukawa suggested that a new particle with "interme- 
diate" mass, eventually called a meson, should be responsible for the strong interaction between 



nucleons 154 



In modern one-meson or one-boson-exchange (OBE) model of the interaction one assumes 
that the basic NN interaction is described by the exchange of all mesons displayed in the right- 
hand part of Fig. |15|. The diagram of that figure can be interpreted as a Feynman diagram 
representing the simplest non-trivial contribution to the two-nucleon Green's function. Such 
contributions can be evaluated by using essentially the same rules as discussed in subsection 
273|. The basic building blocks for evaluating these contributions are the operators for the 
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meson-nucleon vertices, which are given by 



Fs = \'4:7igs for a scalar meson, 

Fps = z\/47r(7ps7^ for a pseudoscalar meson. 



r. = V47r 



gvl'^ + Y~i<J^^ki, for a vector meson. (187) 



The are the usual matrices using the conventions of Bjorken and Drell| |155[ |, the tensor 
operator is defined by the commutator cr^"^ = z[7^,7''] and A;^ refers to the 4-momentum of 
the exchanged meson. Matrix elements of these operators have to be calculated between Dirac 
spinors for plane waves. In the helicity representation these Dirac spinors are given as 



Ea+rn 



where A refers to the helicity, i.e. the projection of the nucleon spin on the direction of the 
momentum q and 



E^ = ^^ + m? (189) 

is the relativistic energy of the free nucleon. Note that the Dirac spinors of Eq. ( |188| ) are 
normalized due to 

S(g, X)u{q, A) = u\q, A)7^(g, A) = 1 . (190) 

Instead of using the pseudoscalar coupling Fpg in ( p.87| ) one also often employs the pseudovector 
coupling 

Tp, = zv^^7 V^/. • (191) 

^mes 

Both couplings yield equivalent results for on-shell nucleons, if one identifies mmes5'ps = 2m/pv. 
For nucleons described by Dirac spinors different from those of a free particle like in ( |188| ), how- 
ever, the pseudovector coupling suppresses the enhancements due to the antiparticle admixture 
as compared to pseudoscalar coupling. 

If now one considers the interaction of two nucleons in the center-of-mass frame with mo- 
menta q and —q before and the momenta q' and ~q' after the interaction, the matrix element 
for the exchange of a meson of the kind a is given by 

V,(g*', q) = {u{-(^)V^u{-q)) P^{k) {u{q')V ^u{q)) , (192) 

with the Dirac spinors u as defined in ( |188| ). Momentum conservation requires that the 4- 
momentum of the exchanged meson is k = q — q' and the meson propagators are given by 

-Ps = To for scalar and pseudoscalar mesons, 

Pv = — ^—7^ — ^— ^ for vector mesons. (193) 

— 

A very efficient way for the evaluation of the OBE matrix elements in ( |192|) using the helicity 
representation for the Dirac spinors has been presented in | 151 ]. The two-particle states can 
be expanded in terms of eigenstates with respect to the total angular momentum J. Since 
the OBE amplitudes are invariant under rotation the angular momentum is a good quantum 
number and one obtains matrix elements 

(A;A'2g'|V|AiA2g)j = 27r F dO sin^ di^,{e) (A'lA'^q'lVlAiAaq) , (194) 

Jo 
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where Aj denotes the hehcity of nucleon i, 6 is the angle between the momenta q and g', and 
the d-ly are the reduced rotation matrices with A = Ai — A2 and A' = A'^ — A'a- Inspection of 
the symmetries for the matrix elements in ( |194|) shows that there are six independent matrix 



elements between the various helicity states for each J. These matrix elements in the helicity 
representation can easily be transformed into the conventional basis of partial waves for two 
nucleons 

|AiA2g)j =^ r^'Ljq), (195) 

where S identifies the total spin of the nucleons, L is the orbital angular momentum of the 
relative motion, which is usually labeled by the letter S, P, D... for L = 0,1,2..., and 
J = L + S* is the total angular momentum. The Pauli principle for the interacting nucleons 
requires that the total isospin Tiso is related to these quantum numbers by the requirement that 
the sum L + S + Tiso is an odd number. 

As discussed before, one can account for the composite structure of the mesons and baryons 



by introducing form factors. This means that one may consider the coupling constants in (187) 
not as universal constants but as depending on the 4-momenta of the interacting hadrons. A 
simple choice for such a form factor, which is commonly used, is to assume that it depends only 
on the momentum transfer k, the momentum of the meson, and takes the form 

Uk)=gJ-^^\ , (196) 

with a cut-off parameter A and an exponent z/, which is 1 for the so-called monopole form 
factor. 

The operators defining the meson nucleon vertices in ( |187[ ) as well as the meson propagators 
in ( |193| ) refer to a relativistic description. This means that the two-body amplitudes which we 
have defined so far should be considered as an irreducible interaction kernel V to be used in a 
relativistic equation like Bethe-Salpeter equation, which can schematically be written 

TBs(g', g; P) = V(g', g; P) + / d^k-^V{q', k- P)gBs{k; P)%s{k, g; P) . (197) 

Note that this is an integral equation in the 4-dimensional space of momentum vectors. The 
momenta of the interacting nucleons are defined in terms of the center-of-mass momentum P, 
which is conserved, and the relative momenta q, q', and k in such a way that the momenta of 
the particles are, e.g., pi = 1/2P ± k. In the center of mass frame the total momentum P 
has a time like component, which is identical to the total energy y/s, with s referring to the 
corresponding Mandelstam variable (see e.g. |[L5(j||), whereas the space component of P is equal 



to zero. The uncorrelated two-particle Green's function ^bs occurring in (|197|) can be written 
as a product of relativistic single-particle Green's functions 

where the superscript (1) or (2) refers to the corresponding nucleon. It is common practice 
to ignore the propagation of the solutions of negative energy and furthermore reduce the 4- 
dimensional Bethe-Salpeter equation ( |1971 ) to an integral equation in three dimensions by 



fixing the time component of /c in a covariant way. One of the possible choices is the approach 



suggested by [157 



^ , i m2 A+«(fc)A+(2)(-fc) , , 
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with A+(*) referring to the projector on Dirac states with positive energy for particle i 



A+«(A;) 



^ , ^ (i) 
'y Ek — ^ ■ k + m 

2m 

Y,\u{k,\i)><u{k,\,)\ (200) 



with u the Dirac spinors of ( |188| ). The assumption of the Blankenbecler-Sugar propagator, 
that the time-hke component ko vanishes, means that for the propagation of the intermediate 
states both nucleons are considered to be off-shell by the same amount and the irreducible 
interaction terms do not transfer energy between the interacting nucleons. This implies that no 
energy transfer should be assumed for the meson propagators if one considers the contributions 
of OBE in the Blankenbecler-Sugar approximation. Replacing the Bethe-Salpeter propagator 
^Bs in ( |197| ) by the Blankenbecler-Sugar approximation one obtains 

/ X ./X /• , ,,m2A+«(A;)A+(2)(_A;)^^, ^ ^ ^ 



where we have replaced s by 4i?g, using Eg = i/g^ + m^. Assuming that the matrix elements of 
V are calculated between spinors, which correspond to the solution of the Dirac equation with 
positive energy (to account for the projectors A+), we may define 




This allows us to rewrite the Blankenbecler-Sugar equation (|201|) as 

Tscat(g", q) = V{^, q) + l -0^,Vi^, k)^-!^--T^,^,{k, , (203) 

which has the form of the non-relativistic Lippmann-Schwinger equation (^ for the scattering 
Tscat matrix. This means that if we evaluate the relativistic matrix elements of V and apply the 
so-called "minimal relativity" factors of (|202|) , the (relativistic) Blankenbecler-Sugar equation 
(|20 1|) becomes identical to the non-relativistic scattering equation. The states of this scattering 



equation can be rewritten in the usual partial-wave basis and the integral ( P03| ) can be solved 
with the techniques as described by 

Alternatives to the Blankenbecler-Sugar approach to reduce the Bethe-Salpeter equation 
to a three-dimensional integral equation have been developed. As examples we mention the 
approaches introduced by Kadychevsky[ p.58|] , Gross ||T59| , Thompson ||160|| , Schierholz ||161|| , and 



Erkelenz[ [151|] . A detailed discussion of the various approaches has been presented by Brown 
and Jackson[ [153|| . 

With the OBE ansatz one can now solve the Blankenbecler-Sugar or a corresponding scat- 
tering equation and adjust the parameter of the OBE model to reproduce the empirical NN 
scattering phase shifts as well as binding energy and other observables for the deuteron. Typical 
sets of parameters resulting from such fits are listed in table |I]. 

Some of the OBE parameters, such as the masses of the vr, t], p, u, and S mesons, are not 



free parameters but are taken from the mass table ||162||. Other parameters, such as the cut-off 
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Table 1: Parameters of the realistic OBE potentials Bonn A, B and C (see table A.l of 0). 
The second column displays the type of meson: pseudoscalar (ps), vector (v) and scalar (s) 
and the third its isospin Tiso- 











Bonn A 


Bonn B 


Bonn C 


Meson 




T- 

ISO 






K 




Aa 




A„ 








[Me V J 




[MeV] 


[MeV] 


[MeV] 






TT 


ps 


1 


138.03 


14.7 


1300 


14.4 


1700 


14.2 


3000 




ps 





548.8 


4 


1500 


3 


1500 







P 


V 


1 


769 


0.86^^ 


1950 


0.9^ 


1850 


1.0^ 


1700 


U 


V 





782.6 


25^^ 


1350 


24.5^ 


1850 


24a 


1400 


5 


s 


1 


983 


1.3 


2000 


2.488 


2000 


4.722 


2000 




s 





550"^ 


8.8 


2200 


8.9437 


1900 


8.6289 


1700 








(710-720)^ 


17.194 


2000 


18.3773 


2000 


17.5667 


2000 



^ The tensor coupling constants are /p=6.1 Qp and /^j = 0. 

^ The a parameters in the first line apply for NN channels with isospin 1, while those in the 
second line refer to isospin channels. In this case the masses for the a meson of 710 (Bonn 
A) and 720 MeV (Bonn B and C) were considered. 



parameters A and the contributions from 77 and 6 exchange, do not effect the fit very much but 
are used as a fine tuning. The coupling constant for the vr is very well constrained by the vrA^ 
scattering data. 

Also the coupling constants for the p meson, in particular the large tensor coupling fp 
are deduced from a dispersion analysis of ttN data in ||163||. As we will see below, this strong 



coupling for the p is of some significance for the nuclear structure calculation. A non-relativistic 
reduction for the p exchange, similar to the one performed for the vr in the preceding subsection, 
yields a tensor component for the NN interaction with a sign opposite to the one deduced from 
one-pion-exchange in ( p07D . Therefore a strong p exchange reduces the tensor force originating 
from the TT exchange significantly. 

The uj coupling constant used in the OBE potential displayed here but also in other realistic 
OBE models is rather large. A simple quark model with SU(3) flavor predicts the uj coupling to 
be nine times as large as that for the corresponding isovector vector meson, the p. The strong uj 
exchange contribution, however, is required to obtain sufficient repulsion for the NN interaction 
at short distances. A possible reason for this discrepancy may be that the uj exchange in the 
OBE model contains an effective parameterization of short-range repulsion originating from 
quark effectsfll]. Another explanation would be that this strong uo exchange simulates also 
repulsive ir — p exchange terms with intermediate isobar excitations [|164|| (see also discussion 
below) . 

The only part of the OBE model, that has a purely phenomenological origin is the a ex- 
change, which describes the medium-range attraction of the NN interaction. This exchange of 
the scalar a meson is used to describe various two-vr exchange processes, which are irreducible 
with respect to intermediate NN states and therefore not accounted for by summing ladder 
diagrams of one-pion-exchange in the Lippmann-Schwinger equation. Such two-vr exchange 



55 



N 



n, p 



a) 



b) 



Figure 16: Contributions to the NN interaction beyond the OBE approach due to interme- 
diate NA (a) and AA excitations (b) 



processes can be studied in a systematic way by means of dispersion relations using experi- 
mental information from vrN scattering processes. Such studies of the medium-range attraction 
employing dispersion relations have been performed e.g. by Vinh Mau and collaborators p7 
and the group in Stony Brook|2^. They are the basis of the medium range attraction, which 
is used in the so-called Paris potential ||165| . 

Another way to determine the effects of these irreducible two-vr exchange processes is to 
evaluate such contributions explicitly. A rather important contribution of this kind are the 
processes displayed in Fig. |T6|, in which the interacting nucleons are excited to intermediate 
A^A or AA excitations (with A representing the isobar excitation of the nucleon at 1232 MeV 
with spin and isospin 3/2). Since the A^ — A excitation requires a change of the spin and isospin 
of the baryon, such excitations could only be formed by the exchange of non-scalar, isovector 
mesons like the vr and p meson. The dominant contribution will be the iterated two-vr exchange 
which yields attraction. In a simple OBE model this attraction would be described in terms of 
a meson exchange. Note that in A^A^ channels with isospin T = only AA intermediate states 
can be reached because of isospin conservation. In T = 1 channels one obtains contributions 
from intermediate A^A and AA states. This is a plausible explanation for the feature that the 
masses of the a mesons listed in table |^ are larger for T = than for T = 1. 

Corresponding terms with vr and p exchange are repulsive and of shorter range [[L64|| . They 
might be simulated by a strong u exchange as discussed above [|^. 

During the past few years, considerable progress has been made in constructing realistic 
models for the NN interaction. In 1993, the Nijmegen group published a new phase-shift analysis 
including selected proton-proton and proton-neutron scattering data below a laboratory energy 
of 350 MeV with a P^r datum of 0.99 for 4301 data ||166| . Based on these data charge- 
dependent A^A^ potentials have been constructed by the Nijmegen group ||167|| , the Argonne 
group|8^ (Argonne Vig) and Machleidt et al. [|168|| (CD-Bonn) which reproduce the A^A^ data 
with a of 1.03, 1.09 and 1.03, respectively. In order to achieve fits with such a high accuracy 
one has to go beyond the OBE ansatz discussed so far. Machleidt et al. ||168|| obtain such an 
accurate fit in the so called CD-Bonn potential by adjusting the parameters of the a meson 
exchange in each partial wave, separately. 

One may say that these charge-dependent A^A^ interactions are essentially phase-shift equiv- 
alent, the on-shell matrix elements of the NN transition matrix T are almost identical. This 
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does, however, not imply that the models for the NN interaction underlying these descriptions 
are identical. Moreover, the off-shell properties of each potential may be rather different. All 
models for the NN interaction V include a one-pion exchange (OPE) term, using essentially 
the same irNN coupling constant. However, even this long range part of the NN interaction, 
which is beheved to be well understood, is treated quite differently in these models. 

The CD-Bonn potential is based on the relativistic meson field theory, which has been 



outlined above. Including the "minimal relativity" factors of ( p02| ) one obtains an expression 
for the One-Pion-Exchange contribution to the NN interaction of two nucleons in the ^5*1 
channel, in plane wave states with momenta k and k' for the initial and final state, respectively, 
which is of the form 

where denotes the momentum transfer, 

q^ = {k- k'f = k^ + k'^ - 2kk' cos e . (205) 



Note that the expression ( |204| ) contains a dependence on the momenta k and k' , which cannot 



be reduced to a dependence on the momentum transfer. This demonstrates that this expression 
for the One-Pion-Exchange yields a non-local interaction. Results for such matrix elements as 



a function of k, keeping k' = 95 MeV/c fixed, are displayed in Fig. p!7|||169 



If now we introduce in (|204D the nonrelativistic approximation for 



E,Ei -m'^ + h'\ (206) 



ignore the form-factor (A^ — m^) / (A^ + q^) and the "minimal relativity" factors, we obtain 

A-n 27rM2 



{ms\k)xoc^\ = '^7^2 I ,dcose [ - ^ ) . (207) 




This can be viewed as the local approximation to the OPE since the matrix element depends 
on the momentum transfer g^, only. It can easily be transformed into the configuration-space 
representation, resulting in a Yukawa term plus a 5 function, which originates from the Fourier 
transform of the constant 1/2 in Eq. ( p07| ). The various steps leading to this result are shown 
in Fig. It is obvious from this figure that all of the steps leading to the local expression 
(|207|) are not really justified for momenta k around and above 200 MeV, a region of relative 



momenta which is of importance in the deuteron wave function. This is true for the matrix 
elements V^^ as well as Vsd- 

Potentials like the Argonne but also the Nijmegen potentials contain the OPE contribution 
in the local approximation regularizing the limit of small r. In the case of the Argonne poten- 
tials the regularization is made in terms of a Gaussian function |p9||. It is remarkable that this 



regularization leads to matrix elements in the SD-channel which are close to those derived from 
the relativistic expression of the Bonn potential. This is not the case in the SS channel, where 
the removal of the 5 function term is very significant. This comparison of the various approxi- 
mations to the OPE part of the NN interaction demonstrates that even this long range part of 
the NN interaction is by no means settled. The local approximation and the regularization by 
form factors have a significant effect. 

The description of the short-range part is also different in these models. The NN potential 
Nijm-II ||167|| is a purely local potential in the sense that it uses the local form of the OPE 
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Figure 17: Plane-wave matrix elements of the one-pion-exchange potential (OPEP) using 
various approximations. As an example, the matrix elements in momentum space {k'\V\k) 
are shown as functions of k for a fixed value of k' = 95 MeV/c. The left part of the figure 
exhibits matrix elements for the partial waves ^Si-^Si, while the right part shows the tensor 
component in the ^Si-^Di channel with k referring to the momentum in the ^Di partial 
wave. 



potential for the long-range part and parameterizes the contributions of medium and short- 
range in terms of local functions (depending only on the relative displacement between the two 
interacting nucleons) multiplied by a set of spin-isospin operators. The same is true for the 



Argonne V^ig potential ^9[. The NN potential denoted by Nijm-I [|1671| uses also the local form 
of OPE but includes a term in the medium- and short-range central-force (see Eq. (13) of 



Ref. ||167|| ) which may be interpreted as a non-local contribution to the central force. The CD- 



Bonn potential is based consistently upon relativistic meson field theory 0. Meson-exchange 
Feynman diagrams are typically nonlocal expressions that are represented in momentum-space 
in analytic form. It has been shown ||168|| that ignoring the non-localities in the OPE part leads 



to a larger tensor component in the bare potential. 

By construction, all realistic NN potentials reproduce the experimental value for the energy 
of the deuteron of -2.224 MeV. However, the various contributions to the total deuteron energy 
originating from kinetic energy and potential energy in the ^Si and ^Di partial waves of relative 
motion, 

E = {^s\T\^s) + {"^dITI^d) + {"^slVl^s) + {"^dIVI^d) + 2{^s\V\^d) 

= Ts + Td + Vss + Vdd + Vsd, (208) 

exhibit quite different results. This can be seen from the numbers listed in Table ^ In this 
table, we display the various contributions to the deuteron binding energy employing the four 
potentials introduced above. 
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rot. 
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Td 
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Vdd 
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Pd 




[MeVJ 


[MeVJ 


[MeVJ 


[MeVJ 


[MeVJ 
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CD-Bonn 


9.79 


5.69 


-4.77 


1.34 


-14.27 


4.83 


Argon 1/18 


i i .Oi 


o.O ( 


-o.yo 


n 77 


1 8 Q/1 


^ 78 
0. ( o 


Nijm I 


9.66 


7.91 


-1.35 


2.37 


-20.82 


5.66 


Nijm II 


12.11 


8.10 


-5.40 


0.59 


-17.63 


5.64 



Table 2: Contributions to the kinetic and potential energy of the deuteron originating from 
the ^5*1 and ^Di parts of the wave function as defined in Eq. |20^ . Results are listed for 



the charge-dependent Bonn potential (CD-Bonn ||168|| ), the Argonne Vi% |^9[, the Nijmegen 
potentials Nijm I and Nijm II [|1671J . The last column of this Table shows the calculated 
D-state probabilities of the deuteron. 



The kinetic energies are significantly larger for the local potentials Vig and Nijm II than 
for the two interaction models CD-Bonn and Nijm I which contain non-local terms. The 



corresponding differences in the S-wave functions can be seen in Fig. |18[ The local potentials 
yield a stronger suppression of the ^5*1 wave function for small relative distances. This reflects 
stronger repulsive short-range components of the local interactions. These stronger short- 
range components are accompanied by larger high momentum components in the momentum 
distribution, which yields larger kinetic energies. 

Comparing the contributions to the potential energy, displayed in Table |^, one finds large 
differences particularly for the tensor contribution Vsd- The dominant part of this tensor 
contribution should originate from the tensor component of the one-pion-exchange potential 
which we discussed above. 

Although the modern NN potentials yield essentially the same NN scattering phase shifts 
and the same binding energy for the deuteron, there are significant differences in the contri- 
butions to both the kinetic and potential energy of the deuteron in the various partial waves. 
Speaking in general terms, these differences can be traced back to off-shell differences between 
the potentials. In particular it is the inclusion of non-local contributions in the long-range 
(tt exchange) as well as short-range part of the NN interaction, which is responsible for these 
differences. 

While the definition of a realistic two-body interaction between nucleons is a rather well 
defined subject with only little differences between the various models, the situation is much 
less clear for three-body and other many-body forces. As an example let us consider the process 
displayed in Fig. p!9|a) with one of the three nucleons being excited to the A resonance in the 
intermediate state. In a many-body theory which does not consider isobar degrees of freedom 
explicitly, this process should be included as two-meson exchange three-nucleon interaction. If 
one calculates the expectation value of this three-nucleon force with the uncorrelated ground 
state of the hole line expansion, the Goldstone diagrams displayed in Fig. |T^) and c) occur. 
For a nuclear system with total isospin T = the contribution in b) vanishes since it represents 
a A - hole (isospin T=l or 2) admixture to the ground state. The diagram of Fig. |T^) is 
closely related to the ground state expectation value of the two-body interaction term shown in 
Fig. |T^a). Fig. |l9|c corresponds to a correction of that two-body diagram which is due to the fact 
that the intermediate nucleon states below the Fermi level are blocked by the Pauli principle. 
Such corrections have been included in nuclear structure calculations, which employed NN 



interactions of the sort shown in Fig. |T6||H0|, |170|| . This demonstrates the model dependence of 
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Figure 18: Wave function for the deuteron (^^i and ^Di) calculated for different realistic 
interactions. 



such three-nucleon terms. 

The models for the three-nucleon interactions which are frequently used in nuclear many- 
body calculations are based on rather phenomenological grounds. Typically they include a 
local parameterization of the two-pion exchange terms with intermediate A excitations of the 



form 171 



Vijk = ^27rX! {{^ij^^ik} {n ' T^j,T^i ' ^k} + ^ [Xij,Xik\ [fj ■ fj-,fj ■ ffc] ) (209) 
cycl 



where 



-CTi ■ a,- + 



3 3 
1 + + 



(210) 
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a) 



b) 



c) 



Figure 19: Contributions to the three-body forces as discussed in the text 
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and the symbols [, ] and {, } denote commutator and anti-commutator, respectively. Further- 
more one introduces a repulsive short range term[p7[|. The parameters are then adjusted to 
obtain a good fit to the binding energies of few-body nuclei and nuclear matter. Various three- 
body forces, labeled e.g. UVII, UIX etc. have been introduced in this way by the Urbana 
group pF 



3.2 Ground state Properties of Nuclear Matter and Finite Nuclei 

In the first part of this section we would like to discuss the convergence of the many-body 
approaches and compare results for nuclear matter as obtained from various calculation schemes 
presented in section 2. The convergence of the hole-line expansion for nuclear matter has been 



investigated during the last few years in particular by the group in Catania|6^, Continuing 
the earlier work of Day||^ they investigated the effects of the three-hole-line contributions for 
various choices of the auxiliary potential U (see Eq. |lT]). In particular they considered the 
standard or conventional choice, which assumes a single-particle potential U = for single- 
particle states above the Fermi level, and the so-called "continuous choice", which has been 



advocated by the Liege group ||Tl||. This continuous choice supplements the definition of the 
auxiliary potential of the hole states in Eq. (|ll|) with a corresponding definition (real part of 
the BHF self-energy) also for the particle states with momenta above the Fermi momentum, 
k > kp. In this way one does not have any gap in the single-particle spectrum at k = kp. 

Fig. shows the results of BHF calculations for the Argonne Vi4^ interaction, assuming 
the standard choice (labeled BHF-s for standard) and the continuous choice (BHF-g) for the 
auxiliary potential. At kp = 1.4 fm~^, which corresponds to a density close to the empirical 
saturation density, the standard choice yields an energy of —10.9 MeV per nucleon while the 
continuous choice leads to —17.1 MeV. At first sight this dependence on the auxiliary potential, 
a difference of 6 MeV seems very large. One should keep in mind, however, that a Hartree-Fock 
calculation at the same density for the same interaction yields +41.8 MeV per nucleon. This 
means that the inclusion of two-hole line contributions yields an attractive contribution of 52.7 
MeV and 58.9 MeV for the standard and continuous choice, respectively. This means that the 
uncertainty, how to choose the auxiliary potential, leads to a 10 percent effect in this two-hole 
line contribution. 

Looking at the Bethe-Goldstone equation (^ it is easy to understand that the continuous 
choice leads to more binding energy. The absolute values for the attractive energy denominators 
in the Bethe-Goldstone equation, which correspond to the excitation energies of two-particle 
two-hole excitations are smaller if one accounts for an attractive single-particle potential also 
in the case of the particle states. This enhances the effects of correlations, leading to more 
attractive matrix elements for the G-matrix, and provides more binding energy. 

The Catania group also evaluated the three-hole line contributions to the energy of nuclear 
matter for both choices of the auxiliary potential U. Results of such calculations, including 
the sum of all Bethe-Fadeev ladders, are displayed by symbols in Fig. The contributions 
of the three-hole line terms are -3.1 MeV and +1.9 MeV per nucleon at kp = 1.4 fm~^ for 
the standard and continuous choice, respectively. This is about a factor 20 smaller than the 
two-hole line contributions discussed above. This indicates a very good convergence of the hole 
line expansion for nuclear matter around saturation density. This conclusion is supported by 
the fact that the results are rather independent on the choice of the auxiliary potential U after 
the effects of three-hole lines are included. Differences can be observed only for densities, which 
correspond to kp around 1.7 fm~^ and larger. At those higher densities the effects of four-hole 
line and higher order terms may get important. Also it is worth noting that the BHF results 
using the continuous choice are closer to the results with inclusion of three-hole line terms. 
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Figure 20: Energy per nucleon calculated for symmetric nuclear matter as a function of 
the Fermi momentum. Continuous lines represent the results of BHF calculations assuming 
the standard choice (BHF-s) and the continuous choice (BHF-c) for the auxiliary potential. 
The symbols (3BF) show results of calculations |^ with inclusion of the three-hole line 
contributions. The results were obtained using the Argonne Vi4^ interactional. 



How do such results obtained from the hole line expansion compare with those derived from 
variational calculations? In order to discuss this question we consider as a first step again the 
results obtained for the Argonne V14 interaction at a density of 0.185 fm~^, which corresponds 
to a Fermi momentum kp of 1.4 fm~^. As we have seen from the discussion above, the hole-line 
expansion yields -15. to -15.2 MeV per nucleon if the effects of three-hole line contributions are 
taken into account. The variational calculation of Wiringa et al.[^ predict an upper bound of 
-13.4 MeV using the FHNC/SOC approximation at this density. The perturbative corrections 
within the framework of the CBF theory (see Eq. ( p.43| ) and subsequent discussion) may provide 
another -1.5 to -2 MeV ||ll^ so that these very different approaches yield very similar results. 

In comparing such numbers one must keep in mind that the results of both approaches 
should be considered with a kind of "error bar" . This error bar should account for conceptual 
problems: On one side we do not know about the effects of higher order terms in the hole line 
expansion, on the other side we don't know about the relevance of e.g. elementary diagrams 
left out in the FHNC/SOC approach. However one should also be aware of uncertainties which 
are due to technical problems. For example, in performing self-consistent BHF calculations one 
typically solves the Bethe-Goldstone equation assuming an angle-averaged approximation for 
the Pauli operator]^, 172|| . Recent investigations show that an exact treatment of the Pauli 



operator may can modify the calculated energy up to around 0.5 MeV per nucleon |[T73| , p.74|| . 
Also details of the parameterization of the single-particle energies can lead to differences of 
this order of magnitude. Uncertainties due to technical reasons must also be considered for 
the variational calculations. As an example we mention the restriction to a specific variational 
ansatz for the wave function. As an indication for uncertainties of variational calculations 
we mention the calculation of the kinetic energy. The kinetic energy can be calculated using 
different expression. If all many-body clusters are calculated completely these expressions yield 
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Figure 21: Energy per nucleon calculated for symmetric nuclear matter and pure neutron 



matter as a function of the Fermi momentum. The results were obtained from BHF ||17(: 
and FHNC[^T^ using the Argonne Vis interaction 



the same results. So the difference in energies obtained using the so-called Jackson-Feenberg 
and the Pandharipande-Bethe expression is a measure of the error ||175| . This estimates yields 



is displayed in Fig. |2T 
potential 



error-bars of the order of 1 MeV per nucleon [p!6| at densities around 0.16 fm~^, the empirical 
saturation density of nuclear matter. 

A comparison between the energies for nuclear matter and neutrons at various densities 
calculated in the framework of the variational approach ||114|| and the BHF approximation |]T7^ ] 

For both calculations the charge-dependent version of the Argonne 
Vis has been used. One finds a rather good agreement between the two ap- 
proaches at densities up to 0.6 fm^'^, which is about 4 times the saturation density of nuclear 
matter. The predictions of the two approaches deviate from each other at higher densities. 
This differences at higher densities are to be expected since the convergence of the hole-line 
expansion should be good at low densities only. Also the variational approach should be less 
accurate at high densities, genuine n-body correlations should become important. The compar- 
ison also shows, however, that quite reliable estimates for the equation of state for symmetric 
nuclear matter as well as asymmetric nuclear systems can be deduced from such microscopic 



many-body calculations up to the densities shown in Fig. 21 



While the techniques developed for the variational calculations impose a restriction to the 
NN interaction to local potentials, more general interactions can be considered in the Brueckner 
hole-line approximation. Therefore in comparing the features in the predictions derived for 
various two-body interactions we will restrict ourselves to the discussion of results obtained 
within the framework of BHF. Results for the saturation points of nuclear matter, i.e. the 
minima of the energy per nucleon versus density curves, are displayed in Fig. The results 



presented include predictions from rather old versions of the NN interaction |^0| like the Reid 
soft-core potentialflT^ (E/A=-10.3 MeV, kF=lA fm'^), but also modern versions like the 
charge- dependent Bonn potential [|168|| (E/A=-19.7 MeV, kF=l-Q8 fm"^) and Argonne VispQl 
(E/A=-16.8 MeV, kF=1.59 fm~^). All these saturation points form a band, the well known 
Coester band, an observation which has been made already in 1970||178|: Some interactions. 
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Figure 22: Saturation points of nuclear matter calculated for various realistic NN interac- 
tions. Some of those results are taken from 



40 1 . Results were obtained using BHF with 



standard choice for the auxiliary potential, as well as the continuous choice. 



like the Reid soft-core, produce a minimum around the empirical saturation density, kp = 1.36, 
but predict an energy of only -10 to -11 MeV per nucleon, which is too small as compared to the 
value -16 MeV taken from the volume term of the empirical mass formula. Other interactions 
predict an energy, which is close to the empirical value but saturate at a density, which is 
around 60 percent too high (or even higher). 

Traditionally, the saturation properties of the different interactions have been related to 
the (i-state probability in the deuteron predicted by this interaction |H0| using the following 



argument: All interactions produce essentially the same phase-shifts for NN scattering, i.e. they 
yield the same T-matrix. One may distinguish "soft" interactions, in which the attraction of 
T originates predominantly from the bare potential V. The attractive terms of higher order, 
reflecting the importance of correlations are of little importance. Such soft interactions show 
a small amount of correlations, which is connected to a small (i-state probability, a direct 
indication for weak tensor correlations. The key-quantity of the BHF approximation is the 
G-matrix, which as we have discussed already in section ^]l|, corresponds to the T-matrix 
except that the terms of higher order in V are quenched because of the Pauli operator and 
dispersive corrections (the absolute values of the energy denominator are larger in the Bethe- 
Goldstone equation for the nuclear medium than in the Lippmann-Schwinger equation). For a 
"soft" interaction this quenching of the attraction in the higher order terms shows only a little 
effect: the G- matrix is almost as attractive as T therefore one obtains large binding energies 
at large densities. In contrast, a "stiff" interaction produces more correlations, the attractive 
contributions of the higher order terms in T are more important, the quenching mechanism is 
more important and one obtains saturation points for nuclear matter at smaller densities and 
smaller energies. 

In order to explore these saturation features a little bit more in detail, we will focus our 
attention to those NN interactions which have recently been fitted with high precision to NN 
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CDRonn 

V. 1 y jj Willi 


ArV18 


Nijixil 


Nijm2 


A 


B 


C 


Rpid 


< E > 


-17.11 


-15.85 


-15.82 


-13.93 


-16.32 


-15.32 


-14.40 


-12.47 


<V > 


-53.34 


-62.92 


-55.08 


-61.94 


-52.44 


-53.03 


-54.95 


-61.51 


<T > 


36.23 


47.07 


39.26 


48.01 


36.12 


37.71 


40.55 


49.04 


<K> 


-22.30 


-40.35 


-28.98 


-28.97 


-12.48 


-26.87 


-45.74 


-27.37 


Mr 


4.64 


30.34 


12.08 


36.871 


7.02 


10.07 


29.56 


176.25 


[% ] 


4.83 


5.78 


5.66 


5.64 


4.38 


4.99 


5.62 


6.47 



Table 3: Energies calculated for nuclear matter with Fermi momentum kp = 1.36 fm~^. 
Results are listed for the energy per nucleon calculated in BHF {< E >) and Hartree- 
Fock (< E >hf) approximation. Furthermore the expectation value for the NN interaction 
< V >, the kinetic energy < Tj^j^ > and the one-pion-exchange term < K- > are listed. For 
completeness we also give the D-state probabihty calculated for the deuteron Pd- Results 
are presented for the charge-dependent Bonn (CDBonn) ||168|| , the Argonne Vis (ArV18) 
and two Nijmegen (Nijml, Nijm2) ||167|| interactions. For a comparison results are also 
given for three older versions of the Bonn interaction (A,B,C) and the Reid soft core 
potential | |177[ , which is supplemented in partial waves in which it is not defined by the OBE 
C potential. All energies are given in MeV per nucleon. 



scattering phase shifts as we already discussed in the preceeding section. In particular we would 
like to explore the contributions of various components of these interactions to the calculated 
binding energy. The BHF approach yields the total energy of the system including effects of 
correlations. Since, however, it does not provide the correlated many-body wave function, one 
does not obtain any information about e.g. the expectation value for the kinetic energy using this 
correlated many-body state. To obtain such information one can use the Hellmann-Feynman 
theorem, which may be formulated as follows: Assume that one splits the total hamiltonian 
into 

H = Ho + AV (211) 
and defines a hamiltonian depending on a parameter A by 

H{X) = Ho + XAV. (212) 

If Ex defines the eigenvalue of 

i7(A)|*A>=^A|*A> (213) 

the expectation value of AV calculated for the eigenstates of the original hamiltonian H = H{1) 
is given as 

(214) 



The BHF approximation can be used to evaluate the energies Ex, which also leads to the 
expectation value < \E^|Ay|\E' > employing this Eq. ( pi4] ). In the present work we are going to 
apply the Hellmann-Feynman theorem to determine the expectation value of the kinetic energy 
and of the one-pion-exchange term AV = contained in the different interactions. 

First differences in the prediction of nuclear properties [|179|| obtained from the modern inter- 
actions are displayed in table ^ which contains various expectation values calculated for nuclear 
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Figure 23: Expectation value of the vr-exchange contribution to the energy of nuclear matter 
function of density 



matter at the empirical saturation density, which corresponds to a Fermi momentum kp of 1.3Q 
fm~^. The most striking indication for the importance of nuclear correlations beyond the mean 
field approximation may be obtained from the comparison of the energy per nucleon calculated 
in the mean-field or Hartree-Fock (HF) approximation. All energies per nucleon calculated in 
the (HF) approximation are positive, therefore far away from the empirical value of -16 MeV. 
Only after inclusion of NN correlations in the BHF approximation results are obtained which 
are close to the experiment. While the HF energies range from 4.6 MeV in the case of CDBonn 
to 36.9 MeV for Nijm2, rather similar results are obtained in the BHF approximations. This 
demonstrates that the effect of correlations is quite different for the different interactions con- 
sidered. However it is worth noting that all these modern interactions are much "softer" than 
e.g. the old Reid soft-core potential [p. 77|| in the sense that the HF result obtained for the Reid 
potential (176 MeV) is much more repulsive. 

Another measure for the correlations is the enhancement of the kinetic energy calculated 
for the correlated wave function as compared to the mean field result which is identical to Tpc, 
the energy per particle of the free Fermi gas. At the empirical density this value for TpQ is 
23 MeV per nucleon. One finds that correlations yield an enhancement for this by a factor 
which ranges from 1.57 in the case of CDBonn to 2.09 for Nijml. It is remarkable that the 
effects of correlations, measured in terms of the enhancement of the kinetic energy or looking at 
the difference between the HF and BHF energies, are significantly smaller for the interactions 
CDBonn and Nijml, which contain non-local terms. The non-local interactions tend to be 
"softer" in the sense discussed above and therefore lead to more binding energy in nuclear 
matter. 

The table |^ also lists the expectation value for the pion-exchange contribution \4 to the 
two-body interaction. Here one should note that the expectation value of \4 calculated in the 
HF approximation is about 15 MeV almost independent of the interaction considered. So it is 
repulsive and completely due to the Fock exchange term. If, however, the expectation value 
for 14 is evaluated for the correlated wave function, one obtains rather attractive contributions 
ranging from -22.30 MeV per nucleon (CDBonn) to -40.35 MeV (ArV18). This expectation 
value is correlated to the strength of the tensor force or the D-state probability Pd calculated 
for the deuteron (see table |^ as well). Interactions with larger P^, like the ArVlS, yield larger 



66 




_Q 25 I ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

r[fm] 



Figure 24: Correlated wave functions I'i/'jj > as defined in (0) as a function of the relative 
distance for the ^5*0 partial wave. Results are shown for a pair of nucleons in nuclear matter at 
empirical saturation density, which heal to an uncorrelated two-nucleon wave function with 
momentum q = 0.96 fm— 1 at larger distances. The curves are labeled by the interactions, 
which were considered. 



values for < \4 >. For a further support of this argument we also give the results for three 
different version of charge-independent Bonn potentials A, B and C, defined in 0. All this 
demonstrates that pionic and tensor correlations are very important to describe the binding 
properties of nuclei. In fact, the gain in binding energy due to correlations from alone is 
almost sufficient to explain the difference between the HF and BHF energies. 

The importance of pionic correlations has been emphasized by Akmal and Pandharipande [|T^ . 
They observe an enhancement of the pionic correlations in FHNC calculations for nuclear mat- 
ter at high densities and interprete this change in the wavefunction as an indicator for a phase 
transition to pion condensation ||1 80| , |181|| . We do not intend to go into a detailed discussion of 



pion condensation ||1 82|| , but simply show the expectation values for < > as a function of 
density in Fig. EBl Indeed one finds a smooth enhancement of this expectation value with den- 
sity if one considers the local representation of the vr-exchange as contained in the Argonne Vig 
interaction. It may be questionable if this should be called an indication for a phase transition. 
Quite a different behavior is obtained if the non-local components are taken into account as 
they appear e.g. in the CDBonn potential. 

Inspecting the expectation values for the kinetic energies we observe a feature very similar 
to the one observed for the deuteron (see section the local interactions, ArVlS and Nijm2, 
yield larger kinetic energies than CDBonn and Nijml, which contain nonlocal terms. This is 
independent of the density considered. 

A different point of view on nuclear correlations may be obtained from inspecting the rel- 
ative wave functions for a correlated pair lipij > defined in ([lOD . Results for such correlated 
wave functions for a pair of nucleons in nuclear matter at empirical saturation density are dis- 



played in Figs. ^ and As an example we consider wave functions which "heal" at large 



relative distances to an uncorrelated two-nucleon wave function with momentum q = 0.96 fm ^ 
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Figure 25: Correlated wave functions as a function of the relative distance for the ^Si and 
^Di partial waves. Further details see Fig. 24. 



calculated at a corresponding average value for the starting energy. 

Fig. ^ shows relative wave functions for the partial wave ^So. One observes the typical 
features: a reduction of the amplitude as compared to the uncorrelated wave function for 
relative distances smaller than 0.5 fm, reflecting the repulsive core of the NN interaction, 
an enhancement for distances between ^ 0.7 fm and 1.7 fm, which is due to the attractive 
components at medium range, and the healing to the uncorrelated wave function at large r. 
One finds that the reduction at short distances is much weaker for the interactions CDBonn and 
Nijml than for the other two. This is in agreement with the discussion of the kinetic energies 
and the difference between HF and BHF energies (see table |]). The nonlocal interactions 
CDBonn and Nijml are able to fit the NN scattering phase shifts with a softer central core 
than the local interactions. 

Very similar features are also observed in the '^5'i partial wave displayed in the left half of 
Fig. ^ For the ^Di partial wave, shown in the right part of Fig. ^ one observes a different 
behavior: All NN interactions yield an enhancement of the correlated wave function at r ^ 
1 fm. This enhancement is due to the tensor correlations, which couples the partial waves 
^Si and ^Di. This enhancement is stronger for the interactions ArV18, Nijml and Nijm2 
than for the CDBonn potential. Note that the former potential contains a pure nonrelativistic, 
local one-pion-exchange term, while the CDBonn contains a relativistic, nonlocal pion-exchange 
contribution. See also the discussion of the wave function for the deuteron in the preceeding 



section 3.1 



and Fig. [18. 

These most recent parameterizations of the nucleon-nucleon interaction not only predict 
different saturation properties of symmetric nuclear matter. Also one can observe differences in 
the calculated symmetry energy at high densities [p. 76|] . The differences between the predictions 
derived from these modern potentials, however, are significantly smaller than those derived 
from older models of the NN interaction. All NN interactions now predict a steady increase 
of the symmetry energy with density. This may be considered as an improvement in the 
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parameterization of the NN interactions. Nevertheless the remaining differences lead to non- 
negligible differences in nuclear astrophysics, in the studies of neutron star matter []TS3|, p.84| . 
The CD Bonn potential predicts a slightly larger symmetry energy at high densities. This leads 
to a larger proton fraction in /3-stable matter. As a consequence one finds that the so-called 
direct URCA process, a very efficient cooling mechanism for neutron stars by neutrino emission, 
may occur at densities around 0.88 fm~^, 1.05 fm~^, or 1.25 fm~^ according to the predictions 
of the CD Bonn, Argonne Vis, or Nijmegen I potential. 

These modern NN interactions account for a breaking of isospin symmetry to reproduce pp 
and pn phase shifts accurately. Even after the electromagnetic effects have been removed, the 
strong interaction between two protons is in general less attractive than the pn interaction. If 
included in the calculation, these isospin symmetry breaking effects lead to differences in the 
predicted energy of the order of 0.3 MeV per nucleon |185| |. Isospin symmetry breaking effects as 
well as charge-symmetry breaking effects (i.e. differences between proton-proton and neutron- 
neutron interactions which are partly due to the different masses of the nucleons) of these 
interactions have also been considered to explore their impact on the so-called Nolen-Schiffer 
anomaly 1 186]. 

Up to this point, we were mainly concerned with bulk properties of nuclear matter and 
correlations which are relevant for all nucleons in the Fermi sea. Special attention has also 
been paid to correlations around the Fermi level. This includes pairing correlations, which have 
been studied in the framework of the BCS approach by various groups [ p!8^ , |188|| . Traditionally, 



such studies concentrate on the 5* = 0, T = 1 pairing effects. Correlations, however, are even 
stronger for the S = 1, T = channel and corresponding pairing effects have been studied within 
the framework of the Green's function approach ||189| , |190|| . Instabilities of normal nuclear matter 
with respect to such pairing correlations are of particular importance at smaller densities, below 
the saturation density of nuclear matter. 

We are now going to discuss the situation of many-body calculations based on realistic 
interactions for finite nuclei. On one hand such calculations are in general much more difficult 
to perform since one has to determine also the appropriate single-particle wave functions within 
the self-consistent calculations. In nuclear matter the basis of single-particle wave functions, 
the plane waves, is determined by the symmetry of the problem. This additional complication 
shows up in all the many-body approaches, which we consider. For the hole-line or Green' 
function approach this implies that one has to determine the basis, in which the single-particle 
self energy is diagonal. The same basis must the also be used to define quantities like the Pauli 
operator in the Bethe-Goldstone equation. One meets a very similar problem in the coupled 
cluster method, for which one has to determine also the amplitude Si in a self-consistent way. 
For variational calculations the additional complication shows up in the fact that one also has 
to determine the long-range part of the wave function, e.g. in terms of an appropriate mean 
field wave function ^mf in the trial wave function of Eq. ( PT] ) from the variational calculation. 

On the other hand, however, the calculations for finite nuclei might lead to more reliable 
results as the average density is lower than for nuclear matter at saturation density. All the 
expansion techniques which we discussed should converge better for small densities, which 
means the results obtained at a given order might be more reliable for the finite system. 

Some features discussed for the infinite systems should also hold for finite nuclei. This 
should be true in particular for short-range correlations. Such short range correlations should 
not be affected by the long-range behavior of the wave functions. To repeat the argument 
from a different point of view: Short-range correlations are described in terms of particle-hole 
excitations of such high energies that the shell structure of the single-particle spectrum of finite 
nuclei is not relevant. The information on such short-range correlations in finite nuclei might 
therefore be deduced from the corresponding information in nuclear matter by means of a 
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Figure 26: Ground state properties for ^^O as predicted in BHF calculations using differ- 
ent OBE potentials from 0. In terms of open symbols we also present the results after 
corrections due to general ring diagrams are included. 



local density approximation [191, 192 , |193| , 194 |. The situation is different for the long-range 



correlations, which are built in terms of low energy configurations for which the shell structure 
is relevant. Such features of the correlations may also vary, depending on the specific nucleus 
under consideration. 

Many-body calculations using the hole-line expansion, the coupled cluster method, the 
Green's function approach or variational calculations have been done for various different nuclei. 
In the following discussion we will mainly focus our attention to results on ^^O, since for this 
double magic nucleus results from the different approaches are available. 

The complication in solving the BHF equations for finite nuclei have been discussed e.g. in 
5B| , where also a computer program for solving the Bethe-Goldstone equations has been made 
available. Different techniques to solve the Bethe-Goldstone equation have been described by 
Barrett et al.[ 195 and Kuo and coworkers] 196|| . Also in BHF calculations for finite nuclei one 
has to choose the spectrum of the intermediate particle states in the Bethe-Goldstone equation. 
The standard choice is the same as for nuclear matter, use pure kinetic energies. Old studies 
of three-body terms in the coupled cluster method led to the suggestion that a shift of the 
kinetic energy spectrum by a value around -8 MeV, would minimize the effects of three-body 
correlations] 197 |. Such results indicate that the effects of three-body correlations may even 
be smaller than for nuclear matter (see the convergence argument given above), however, a 
systematic study of three-hole line contributions for various interactions should be performed 
and is still missing. 

Results for the energy and the radius of the charge distribution of ^^O obtained from BHF 
calculations ||198|| using various meson exchange potentials are displayed in Fig. The plot 



of the results in terms of the energy per nucleon versus the inverse of the radius of the charge 
distribution has been chosen to have a presentation of results which is similar to the plot 
of saturation point for nuclear matter in Fig. ^ The results for nuclear matter and finite 



nuclei are very similar in so far that also the calculated ground state properties for ^^O show 
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a behavior like the Coester band discussed for nuclear matter. There are interactions like the 
OBEP A (defined in Table A.l of 0), which yield about the correct binding energy per nucleon 
of -7.98 MeV but predict a radius of the charge distribution, which is much smaller than the 
experimental one. Therefore we also have in this case a density which is too large. Other OBE 
interaction, which show a larger D-state probability in the deuteron, like OBEP B and OBEB 
C, are "stiffer" and lead to slightly larger radii but less binding energy. Note that all these 
modern OBE potential yield more binding than the old Reid soft-core potential ||177| . 

The BHF approximation can be considered as a very efficient tool to account for the ef- 
fects of two-nucleon short range correlations. This approach, however, might be too simple to 
account for long-range correlations corresponding to low energy configurations in the nuclear 
shell-model. Therefore attempts have been made to consider a model space, in the sense as we 
introduced this term in the discussion around Eq. (^, in which these long range are treated ex- 
plicitly, and use the lowest order hole-line expansion to account for the short-range correlations 
outside this model space. This is similar to the concept of model space for nuclear matter as it 



has been discussed by Kuo et al. One possible way along this line is to consider a 

solution of the Bethe-Goldstone equation (|^) using a Pauli operator which restricts intermediate 
two-nucleon states to configurations which are outside the model-space considered. One may 
then use the resulting G-matrix in shell-model configuration mixing calculation, in which all 
nucleons are considered as active particles which are allowed to form all configurations within 
the model space. Such kind of no-core shell-model calculations have recently been done by Bar- 



rett and coworkers ||199| , |200|| . Due to the dramatic increase of configurations with increasing 
number of active nucleons, however, such studies are restricted to very light nuclei. 

For heavier nuclei one can consider a model space defined in terms of configurations in e.g. an 
oscillator basis and consider in a systematic way the effects of correlations defined within this 
model space. One important class of such long-range configurations are the RPA correlations 
represented by the so-called particle-hole ring diagrams. They correspond to an iteration of 
the particle-hole ladders to any order. Within such a limited model space the number of 
particle-particle configurations is of the same order as the number of hole-hole configurations. 
Therefore one may like to treat particle-particle and hole-hole ladders on the same footing as 
it is done in the Green's function formalism. This is achieved by summing the contribution of 
all particle-particle hole-hole ring diagrams. A technique has been developed which allows the 
consistent summation of all particle-hole and particle-particle hole-hole ring diagrams leading 
to a so-called Super RPA (SRPA)||70[|. It is remarkable that the resulting SRPA equations yield 
a stable solution only if the self-energy, used to calculate the single-particle Green's function. 



is determined in a self-consistent way 1 201 



Such long-range correlations lead to an additional binding energy for closed shell nuclei 
like ^^O of around 1 MeV per nucleon. They might be characterized by the depletion of the 
hole-state occupation. Again for the example of ^^O, the SRPA correlations in a model space 
which includes the IpOf oscillator states leads to an occupation probability for the Op hole 
states which is of the order of 0.85. Note that this depletion of the hole state occupation due 
to the long-range correlations should be added to the depletion which is due to the short-range 
correlations (see discussion in the next subsection). The long-range correlations also have an 
effect on the calculated radii. The effects of these SRPA correlations on the ground-state 
properties of ^^O are also shown in Fig. |26| for the interactions OBEP A to C. The SRPA 
correlations yield a small improvement of the results obtained within the BHF approach. This 
improvement, however, is not sufficient to meet the experimental data. 

Calculations for finite nuclei using the "Exp(S)" or Coupled Cluster Method have been 
performed already 25 years ago by Zabolitzky|^, Using the Reid soft-core potential he 
finds that the SUB2 approximation yields results which are very similar to those obtained 
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in the BHF approximation (see result for the Reid potential in Fig. ^) and that the effects 
of three-body correlations, included in the SUB3 approach, leads to an additional binding 
energy of around 0.5 MeV per nucleon connected with an enhancement of the wound integral 
or depletion of the hole-state occupation by around 0.05 percent. This indicates a very good 
convergence of the Coupled Cluster Method as well as the hole-line expansion. 

More recent calculations using the language of the Coupled Cluster Approach have been 
reported by Heisenberg and Mihaila[]73[ for the Argonne Vis interaction [^]. The calculations 
are performed in a configuration space defined in terms of harmonic oscillator functions. They 
use a G-matrix approximation for the two-body interaction and calculate a mean field which 
they correct to account for 3p3h and 4p4h correlations. They obtain for the case of ^^O a binding 
energy of -5.9 MeV per nucleon and a radius of the charge distribution of 2.81 which is even 
larger than the experimental value. Unfortunately, there are no BHF calculations available for 
the same interaction. However, looking at the results obtained for nuclear matter, one would 
expect that BHF calculations for the interaction might produce results with a binding 
energy around -6 MeV but a smaller radius. It is not clear whether these possible differences 
are due to the 3p3h and 4p4h corrections included in |7^, due to the restricted oscillator space 
or due to a lack of self-consistency in solving the Bethe-Goldstone equation. 

Variational calculations using the FHNC summation techniques which we discussed in sec- 
tion |2]^ have recently been performed for our reference nucleus ^^O as well ||113|| . As we have 
discussed before the techniques of these variational calculation restrict their applications to 
local interactions. In the work of Fabrocini et al. | |1 13[ various versions of the Argonne po- 
tentials Vi4 and Vg have been used supplemented by three-nucleon potentials of the Urbana 
group. One of the main problems for such variational calculations for finite systems is to ob- 
tain a reliable description for the long range structure of the wave function, i.e. the shape of 
the single-particle wave functions which form the mean field part of the trial wave function 
([91|). In ref. ||113|| two different models for these single-particle wave functions are considered: a 
harmonic oscillator model and a Woods Saxon parameterization. The optimal parameters are 
determined from the variational calculation. It turns out, however, that the functional exhibits 
a rather fiat minimum as a function of the parameters characterizing the single-particle waves. 
This is mainly due to the balance between kinetic energy and potential energy. In particular 
for wave functions with small radii one observes a cancellation between these contributions, 
which is rather sensitive to the details of the variational form of the wave function. 

In fact, one finds minima of the variational calculations, which are quite different for the 
two parameterization of the single-particle wave functions considered. Therefore also the ex- 
pectation values of single-particle operators like the radius are rather model-dependent. The 
Woods Saxon parameterization yields a radius for ^^O, which is 0.2 fm smaller than the one 
derived from the harmonic oscillator parameterization. This certainly a drawback of variational 
calculations for finite systems. On the other hand, however, the variational calculation yield 
upper bounds for the energy of the order of -5.2 MeV per nucleon in the case of ^^O which is 
a very good benchmark for other many-body calculations. This is supported by the fact that 
the FHNC/SOC calculations yield results very close to corresponding Variational Monte-Carlo 
calculation in cases for which such a comparison is available [[L13|| . 

Very extensive extensive Variational Monte Carlo (VMC) as well as Green's function Monte 
Carlo (GFMC) calculations have recently been performed for nuclei with particle number up to 
A = 8 ||134| , |142|| . Results of such calculations for the ground state of these nuclei are presented 
in table The calculations leading to the results displayed in that figure used the Argonne 
Vi8 potential for the two-body interaction supplemented by the Urbana IX three-nucleon force. 
This version IX of the Urbana three nucleon forces has been adjusted to reproduce together 
with the Vis NN interaction the binding energy of ^H and to give a reasonable saturation point 
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Energy [MeV] 


Radius [fm] 




VMC 


GFMC 


Exp. 


VMC 


Exp. 




-8.32 


-8.47 


-8.48 


1.59 


1.60 


4tt 


-27.78 


-28.30 


-28.30 


1.47 


1.47 


°He 


-24.87 


-27.64 


-29.27 


1.95 




6t • 


-28.09 


-31.25 


-31.99 


2.46 


2.43 


'He 


-24.43 


-25.16 


-28.82 








-32.78 


-37.44 


-39.24 


2.26 


2.27 


rie 


-19.71 


-25.77 


-31.41 








-29.70 


-38.26 


-41.28 






8Be 


-48.06 


-54.66 


-56.50 







Table 4: Results for ground state properties of nuclei with nucleon number up to A = 8 
as obtained in VMC and GFMC calculations. These results for the total energy and rms 
radii for the protons have been obtained using the Argonne Vig for the two-body interaction 
supplemented by the Urbana IX three-nucleon force. The results displayed in this table have 
been copied from [|134|| and [142 . 



for nuclear matter. Because of these adjustments the calculated energies cannot directly be 
compared to the results which we discussed before in which no three-nucleon force has been 
employed. The expectation value of the three nucleon force alone yields a contribution to the 
energy of ^Li of -8.9 MeV. 

Despite the use of this adjusted three-nucleon force it is quite remarkable to see how well the 
GFMC calculations reproduce the experimental values for the binding energies. It is also very 
satisfactory to see that he variational calculations (VMC) yield predictions for the energies 
which are above those obtained in GFMC by values which are typically less than 1 MeV 
per nucleon. Here one must keep in mind, however, that for these very light nuclei more 
sophisticated shapes could be considered for the trial wave function than it was possible for 
heavier nuclei. In particular one does not start assuming a set of single-particle wave function 
in the mean field part of the trial function, but determines also the long range part of the radial 
shape in terms of correlation functions. 

Also the results for the radii are in very good agreement with the empirical data. Again one 
should keep in mind, however, that this success is at least partly due to the use of the adjusted 
three-nucleon force, which has been determined to remove the problem of the Coester band in 
nuclear matter. This indicates that sophisticated variational calculations using a hamiltonian 
which has been adjusted to reproduce the saturation point in nuclear matter tend to give proper 
saturation mechanisms also for finite nuclei. This leads to the expectation that hamiltonians, 
which lead to a saturation of nuclear matter at a too high density shall predict radii for finite 
nuclei, which tend to be too small. 



3.3 Single Particle Properties in Nuclear Matter and Finite Nuclei 

Single particle properties are most conveniently described in terms of the Green's function 
which has been introduced in section |2.3| . All single particle properties can be calculated from 
the single-particle Green's function, in the sense that the expectation value of any one-body 
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operator can be obtained by using the hole spectral function (Eq. (|75D). In this section we will 
mainly discuss the properties of the hole spectral function in finite nuclei. However, we will 
also present some results in nuclear matter which show some features of the spectral functions. 
A recent review on the effects of correlations in the independent particle motion in fermion 
systems, with special emphasis in the nuclear case, can be found in Ref. p02[ . 

The physical meaning of the hole spectral function can be read from Eq. ([72|) : it gives the 
probability of removing a particle with momentum k from the system of A particles leaving the 
resulting (A-1) system with an energy E^~^ = Eq — u, where Eq is the ground state energy of 
the A particle system. Analogously, the particle spectral function Sp{k,uj) is the probability of 
adding a particle with momentum k and leaving the resulting {A + l)-system with an energy 
E^~^^ = to + Eq. (see also Fig. |l^ and discussion in section 



In the Self- Consistent Green's Function (SCGF) formalism, the self-energy is the key quan- 
tity to determine the one-body Green's function. The self-energy takes into account the strong 
interactions that a nucleon in the nuclear medium has with the other nucleons. Since the 
self-energy is determined from the two-body effective interaction between two dressed particles 
(Fig.(|llj)), and this requires in turn the knowledge of the propagator, one needs to deal with a 



coupled problem which must be solved self-consistently |Tj]. In the case of an infinite system, 
the self-energy is diagonal in k and the formal solution of the Dyson's equation (Eq.([78|)) is 
particularly simple 

gik,uj) = ^ . . (215) 



By combining the Lehmann representation ([7TD, which expresses the single particle prop- 
agator in terms of the spectral functions, with the Dyson equation ([78|) , which relates the 
propagator to the self-energy, the following expressions for the spectral functions are obtained 
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1 Imll(k,uj) 

Due to the short range repulsion present in any reahstic interaction, a meaningful approx- 
imation to the effective interaction, used in the calculation of the self-energy, should sum up 
ladder diagrams which in this approach include also hole-hole propagation, either to all orders 
189, n] or to a second order in the propagation of holes | [2Uij| , 194|. The inclusion of hole- hole 



propagation is required since at the single particle level it yields the coupling to the excited 
states of the (A-1) particle system. The self-energy is used in the Dyson equation and the re- 
sulting single particle propagator is plugged again in the calculation of the effective interaction, 
until self-consistency is achieved. In the calculations presented here, self-consistency has been 
established only for the quasi-particle energies, i.e. at the level of the real on-shell part of the 
self-energy 

^ipik) = ^ + Re^ik, e,pik)). (217) 
Zm 

When this self-consistency is achieved, the complete energy dependence of the self-energy can 
be studied [[189|| . The imaginary part of the self-energy is different from zero over a wide range 



of energies, being positive below ep and negative above. This yields a spreading of the single- 
particle strength, but ioi k < kp Sh will still contain a peak at the quasi-particle energy. Since 
Sp contains some fraction of the total strength the occupation probability (Eq. (^) ) ioi k < kp 
will be depleted to a value below one. 

The energy dependence of the hole spectral function is shown in Fig. ^for several momenta. 
These spectral functions have been obtained for the Reid soft core potential at kp = 1.36 fm~^ 



in the framework of SCGF ||189[ . The hole spectral function in nuclear matter has been also 
calculated for the Urbana V14 interaction |^ in the framework of CBF theory, by including 
one hole and 2hlp intermediate correlated states p04| . The agreement between both methods 



is rather satisfactory. For k < kp the hole spectral function shows a sharp peak (quasi- hole 
peak) located at the quasi-particle energy that concentrates most of the strength. When k 
approaches kp the quasi-hole peak becomes sharper. In the limit k = kp, the peak is just a 
(5-function, the strength of which defines the discontinuity of the momentum distribution at 
kp. The typical occupation probability for k < kp is about 80%. Since the number of particles 
must be conserved, it is necessary to have partial occupation of states k > kp, which were 
unoccupied in the mean field description. A smooth distribution of strength is observed for 
k > kp indicating the possibility to find a nucleon in a momentum that would be empty in 
absence of correlations. 

A very important question concerns the location of the high momentum components in the 



nucleus as a function of the energy. To illustrate this point, Fig.^ shows the occupations for 
k > kp obtained by integrating Sh{k,uj) over excitation energies up to 100 Mev, 200 Mev and 
up to infinity. One concludes that it is necessary to integrate over high excitation energies of 
the A — 1 system in order to catch the high-momentum components. As a result of correlations 
there are particles promoted above the Fermi sea and n{k) is characterized by a depletion below 
kp and a long tail (occupation above kp) which gives ~ a 62% of the total kinetic energy. The 
distribution, when the nuclear matter is considered as a normal Fermi fluid, is discontinuous at 
kp and the discontinuity is given by Zkp- Basically all calculations, independently of the method 
and of the interaction used, provide rather similar results for the depletion. The result for the 
occupation at A; = using the Reid soft core at kp = 1.36 fm~^ within SCGF is n(0) = 0.83 



189|| . CBF gives the same result for Urbana V14, however the calculation was performed at 



slightly smaller density kp = 1.33 fm ^ [|123||. The calculations in the framework of Green 
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Figure 28: Momentum distribution for k > kp obtained with various energy cut-offs in the 
integration of the spectral functions. 



function formahsm but using a self-energy with up to second order in the propagation of holes, 
give also very similar results | |194| , pu^| . However, the discontinuity at kp is more sensitive to 
both the interaction and the method. Under the same conditions as before, Z^^ = 0.72 for the 
Reid and 0.7 for the Urbana V14. We will have further discussions on n{k) in connection with 
the momentum distribution in ^^O. 

The spectral function fulfils energy weighted sum-rules, which besides providing useful infor- 
mation on the spreading of the strength can also be used as a way to control the approximations 
employed in the calculation. The two lowest sum-rules have a simple expression. mo{k) is ob- 
tained by using the completeness relation in the energy integration of the spectral functions, 

mo{k) = jXo dujSh{k, u) + duSp{k, u) 

(vl/^ I {ak,al} I vl/^) = l, (218) 

due to the fermion character of the nucleons, the anti-commutator jok, aj^j is equal to the unit 
operator. 

A similar procedure, leads to the first-order energy weighted sum rule mi{k): 
mi{k) = JX^ujSh{k,uj)duj + ujSp{k,uj)duj 

= {^i\{[a^,H],ai}\^^). (219) 

In order to evaluate the right hand side of mi{k) it is necessary to assume a Hamiltonian, 
which in the present case is taken to be nonrelativistic and with only two-body forces. Under 
these assumptions, 

« I {[a^,H],ai} I O = |!- + ^ J d'k'n{k'){Kk' \ V \ k,k%, (220) 

where (k, k' | | k, k')^ is the anti-symmetrized two-body matrix element of the bare nucleon- 
nucleon interaction. The second term on the right hand side of Eq. (|220| ) can be identified with 
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Figure 29: Momentum dependence of the removal energy and the quasi-particle energies, 
in nuclear matter at kp = 1.36 fm~^. 



the energy independent part of the self-energy which is obtained as the high frequency limit of 
Tjlkyu) [|203|| . These sum-rules have been investigated in Ref. ||205|| . There, it was found that 
the SCGF approach yields spectral functions that fulfil both sum-rules rather accurately. Due 
to the short range repulsion of the potential and to the positive character of the kinetic energy 
contribution, the right hand side of Eq. ( |220| ) is usually positive definite and rather large. For 
instance using the Reid potential it is around 300 MeV at saturation density for k = 0. For 
momenta below kp , most of the strength 80%) is below tp which is a negative quantity. 
Therefore the contribution of the hole part is negative and it must be the positive contribution 
of the high energy tail of the particle part of the spectral function which, besides compensating 
the negative contribution of the hole part, brings the left hand side to fulfil the equality. This 
confirms the need for the appearance of single-particle strength at very high energies. 
Another useful check for the spectral functions is the density sum-rule 



P 



deg 
(27r)3 



(fk n{k) 



deg 
(27r)3 



(fk 



Sh{k, uj)duj . 



(221) 



with deg denoting the degeneracy of the single-particle level, which is 4 for nuclear matter. 
This sum-rule is well respected by the different approaches. In the case of the SCGF formalism 
it is fulfilled within 1%. The integration up to kp provides 81% of the sum-rule and one gets a 
97% if the integration is carried up to "ikp. 

In the case of two body interactions, the hole spectral function gives access, through the 
Koltun sum- rule ([76|) , to the binding energy per particle. 
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to the kinetic energy per particle. 
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and to the potential energy per particle v{p) = e{p) — t(p). Introducing the removal energy 
er{k), 

Mfc) = ■^-^,'^"/f-r' (224) 



one can express the binding energy as 



In Brueckner-Hartree-Fock, er{k) = k"^ /2m + u{k), i.e. it coincides with the quasiparticle 
energy, and n{k) = 9{kp — k), therefore 

Due to correlations, the spreading of the spectral function implies that e,.(/c) 7^ ^qp{k). Obviously 
this comparison between er{k) and eqp has sense only for k < kp. 

The momentum dependence of the removal energy is shown in Fig. ^ together with the 
quasi-particle energies (below kp) for the Reid soft core within SCGF. er{k) is an increasing 
function of k for k < kp and decreases quite fast for k > kp having a large discontinuity at 
kp. For k = 0.027 , er = —121 MeV, while eqp = —70.38 MeV, and the occupation probability 
n(0.027) = 0.83. 

Averaging ( p24|) over all momenta one can define a mean removal energy, e^, by 

. Jd'kZl,du;u;S,{k,u;) . ^ , . . 

In the case that we are considering as an example, i.e. nuclear matter with the Reid soft core, 
at kp = 1.36 fm~^ in the framework of SCGF, ^ —86 MeV. On the other hand, the binding 
energy is ~ — 18MeV, which is too large. However, what we want to show here is the balance of 
kinetic and potential energy and the contribution to the Koltun sum-rule from momenta above 
kp. 

This binding energy is the result of a cancellation between a kinetic energy, t{p) = 48.4 
MeV, and a potential energy ,f (p) = —66.7 MeV. The contributions to t(p) and v{p) of k's 
below kp are the 38% and the 47% respectively. This supports the crucial role of the high 
momentum components also in the energetics of the system, when the kinetic and the potential 
energy are calculated separately. One should keep in mind that the contribution of the low 
momenta {k < kp) to the density sum-rule was ~ 80% and their contribution to the total 
energy is around 70%. The details depend on the interaction used. Note that the expectation 
value for the kinetic energy obtained in SCGF agrees rather well with the result listed in table 
^ for the Reid potential. The comparison in that table also shows that the non-local OBE 
interactions are softer than the Reid potential leading to smaller kinetic energies and smaller 
contributions from high momenta. 

To carry on the full self-consistent scheme, and calculate the effective interaction with 
dressed particles, i.e. considering the spectral functions in the intermediate states in calculating 
the G-matrix is a very difficult task. There are however, recent attempts in that direction 



2061, |207|] and also efforts to analyze the meaning of using dressed particles when one considers 



the scattering in the medium |P08 



The calculation of the single particle spectral function for finite nuclei is much more compli- 
cated. Both the self-energy and the Green's function are not anymore diagonal in momentum 
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and furthermore, one should deal with discrete (bound) and continuum single particle basis 
states. Microscopic calculations have mainly been performed for very light nuclei |f^09[ . For 



heavier nuclei, several procedures based on local density approximation have been devised 
191, |192| , 19ij|| . Here we will pay more attention to the explicit calculations of the spectral 



function directly in finite nuclei. The calculations have been mainly dealing with the hole part 
of the spectral function for ^^O nucleus for which also systematic exclusive electron scattering 
measurements exist, which will be discussed in the next section. A convenient way to take ad- 
vantage of the spherical symmetry is to introduce a partial wave decomposition of the spectral 
function and work in the single-particle basis characterized by the orbital angular momentum 
/, total angular momentum j, isospin third component r and momentum k, 

Si.Ak, ^) = E I i^n~' I ^Mr I P S{u - {E^ - Ef)) (228) 

n 

The momentum distribution 

ni.Ak) = « I al^^a^i.r \ ^t) = E I I ^ki,r \ ^t) T, (229) 

n 

is obtained by integrating Sijr{k,u!) over the excitation energies of the A — 1 system. In the 
independent particle model (IPM), the sum in this equation is typically reduced to one term, if 
(/,j, r) refer to a single-particle orbit occupied in \1/q. In that case, Eq.( p29| ) yields the square 
of the momentum space wave function for this single particle state. The contribution nijr{k) 
vanishes in the IPM if no state with this quantum numbers is occupied. If correlations are 
present beyond the IPM approach this simple picture is no longer true and one can study the 
effects of correlations on n{k). Taking into account the degeneracy factors of each orbital one 
obtains the total momentum distribution 

n{k) = J2{2j + l)nijrik). (230) 

l,j,T 

The spectral function for the various partial waves, Sijr{k,uj) can be obtained from the 
imaginary part of the corresponding single-particle Green's function gijT-{k,uj) which can be 
evaluated by solving the Dyson equation (Eq. (^)) 



gij{ki,k2;uj) = gij'{ki,k2;uj) + J dk^ J dk4glj'{ki,k3;uj)AJ:ij{ks,k4;uj)gij{k4,k2;uj), (231) 

where g^'^^ refers to a Hartree-Fock propagator and AS^j represents contributions to the real and 
imaginary part of the irreducible self-energy , which go beyond the Hartree-Fock approximation 
of the nucleon self-energy used to derive g^^\ Notice that here and in the following we have 
dropped the isospin quantum number r. Ignoring the Coulomb interaction between the protons 
the Green function are identical for N = Z nuclei and therefore independent of the quantum 
number r. 

As in the nuclear matter case, the key point of this approach is the calculation of the self- 
energy, its calculation as well as the solution of the Dyson equation has been discussed in detail 



in Ref. However we include here a brief summary of the relevant aspects of the method. 

The self-energy is evaluated in terms of a G matrix which is obtained as a solution of the 
Bethe-Goldstone equation for nuclear-matter. The Bethe-Goldstone equation is solved at a 
given density and starting energy. In the calculations reported here, the chosen interaction is 
the one-boson- exchange (OBE) potential B p| for kp = 1.4/m~^ and starting energy -10 MeV. 
These choices are rather arbitrary. It turns out, however, that the final result is not sensitive 
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Figure 30: Momentum distribution (solid line) for different partial waves in ^^O. The 
momentum distribution is the sum of the quasihole contribution (short-dashed line) and the 
continuum contribution (long-dashed line). 



to this choice. The self-energy contains a Hartree-Fock contribution plus terms of second order 
in G with intermediate two-particle-one-hole (2plh) and two-hole-one-particle states, assuming 
harmonic oscillator states for the hole states and plane waves appropriately orthogonalized for 
the particle states. The techniques, which are needed to evaluate matrix elements in such a 
mixed representation are described in | |22| ] . 



Vif'\k, k'- u) = - ; ' ' ' du', (232) 



Using such matrix elements one can determine the imaginary part of the 2plh and 2hlp 
contributions, W^J'^^ and Wf^^^^ in a straightforward way. 

The real parts of the self-energy are calculated by means of dispersion relations like, 

'if K]uj} = — / 

where V indicates a principal value integral and /c, k' refer to the modulus of the single-particle 
momentum k for the nucleon under consideration with angular momentum quantum numbers 
I and j. The imaginary part W"^"^^^ is different from zero only for positive energies. Since the 
diagonal matrix elements of W"^^^^ are negative, the dispersion relation (Eq. ( P32|) ) implies that 
the diagonal elements of V"^^^'* will be attractive for negative energies. They will decrease and 
change sign only for large positive values for the energy of the interacting nucleon. 

Since the Hartree-Fock contribution is calculated in terms of a nuclear G-matrix, it 
already contains 2plh terms. In order to avoid such an over-counting of particle-particle ladder 
terms, a correction term (Vc), which contains the 2plh contribution calculated in nuclear matter 
with the same starting energy and Pauli operator as used in the Bethe-Goldstone equation, is 
subtracted from ^^pih ^ dispersion relation similar to Eq.( p32| ) holds for v"^^^^ and W"^^^^. 

Summing up the different contributions, the self-energy is expressed as 

S = S^^ + AS = S^^ + (\/2pi/^ _v^ + ^ i{W^P^^ + W^'^^P). (233) 
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The resulting self-energy is non-local and energy dependent. It represents a microscopic 
derivation of the optical potential which can be used to analyze the data of nucleon-nucleus 
scattering. Careful studies of the non-locality of this self-energy have been presented in Ref. 
21U| . In this reference, a parameterization of the self-energy, for both the real and the imaginary 



part, was given in terms of local energy dependent potential of the Woods-Saxon form. 

Once the self-energy is calculated one must solve the Dyson equation (Eq. ( p31| )) for the 
single particle propagator. To this aim, it is useful to discretize the integrals involved in this 
equation by considering a complete basis within a spherical box of a radius Rbox- The calculated 
observables are independent of the choice of Rbox, if it is chosen to be around 15 fm or larger. 
A complete and orthonormal set of regular basis functions within this box is given by 

^ujmir) = (r I hljm) = Najiihr)yijmiO, 0). (234) 

In this equation yijm represent the spherical harmonics including the spin degrees of freedom. 
Nil is a normalization constant and ji denote the spherical Bessel functions for the discrete 
momenta ki which fulfil 

MhRbox) = 0. (235) 

These basis functions, defined for discrete values of the momentum ki within the box, differ 
from the plane wave states defined in the continuum with the same momenta just by the 



normalization constant, which is y2/7r for the latter. 

As a first step, one constructs the Hartree-Fock approximation for the single-particle Green's 
function in the "box basis". To this end, the Hartree-Fock Hamiltonian is diagonalized, 

E I ir^in + T.^/ I fc„)(fcn I = eSf {h I (236) 

n=l 

The set of basis states in the box is truncated by assuming an appropriate N^ax (Usually 
Nmax = 20 is enough). In the basis of Hartree-Fock states | a), the Hartree-Fock propagator is 
diagonal and given by 

where the sign in front of the infinitesimal imaginary quantity irj is positive (negative) if e^^^ 
is above (below) the Fermi energy. With these ingredients one can solve the Dyson equation 
(Eq. ( [^31| )). An efficient procedure is to determine first the reducible self-energy, 

(a I S[f (^) I (3) = {a I ASy(^) | /3) + I AS,,(^) | 7) x gl^\r,u;){j \ S[f (^) | f3) (238) 

7 

and obtain the propagator from 

gi,{a,P;u;) = 6^pgl^\a;u;) + gl^\a;u;){a \ ^f{u;) \ (3)gf^\p-,uj). (239) 

Using this representation of the Green's function one can calculate the spectral function in 
the "box basis" from 

Sij{km,kn;uj) = ^Im (^{km I a)ijgij{a,p;uj){p \ kn)ij^ . (240) 

For energies below the lowest sp energy of a given Hartree-Fock state (with Ij) this spectral 
function is different from zero only due to the imaginary part of S^'^'^. This contribution involves 
the coupling to the continuum of 2hlp states and is therefore 

2m' 



E (^^ I TTT^- + + ASz,(a; = e%) \ kr.){k„ \ T),, = e%{h \ T),,. (241) 



n=l 
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Figure 31: The total momentum distribution of ^^O. The quasihole contributions also 
shown together with the results obtained with various energy cutoffs in the integration of 
the spectral functions. 



Since, in this approach AS contains a sizable imaginary part only for energies below e^^*, the 
energies of the quasihole states are real and the continuum contribution is separated in energy 
from the quasihole contribution. The quasihole contribution to the hole spectral function is 
given by 

Sfijikm, kn] uj) = Zrij{km I T)ij(T | kn)ijS{u - efij), (242) 



with the spectroscopic factor for the quasihole state given by 

Finally, the continuum contribution of Eq. ( |240D and the quasihole parts of Eq. ( p^43| ) can be 
added and renormalized to obtain the spectral function in the continuum representation at the 
momenta defined by Eq. 



SlJik^,uJ) = (^Sf^{h,u)+Y.Sti^{h,uj)j . (244) 

The quasi-hole states with a probability defined by the spectroscopic factor Zij (in this 
example 0.78, 0.91 and 0.90 in the case of Si/2,P3/2 and Pi/2) are located at the corresponding 
quasi-hole energies -34.30 MeV, -17.90 MeV and -14.14 MeV, respectively. Some strength has 
been moved to more negative values of u in the (2hlp) continuum. 

The momentum distribution nij{k) is given by the energy-integrated spectral function. Here 
it is also convenient to separate the contribution in two pieces: the continuum and the quasi-hole 
part 

(k) = r duo ISI^ {k, u) + ~Sl'^{k, uo)] . (245) 

This separation into the two parts is displayed in Fig. ^ for various partial waves. This figure 
shows quite clearly that the momentum distribution at small momenta is dominated by the 
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Figure 32: Momentum distribution of ^^O, obtained in the Green function approach em- 
ploying the OBE potential, compared with n{k) obtained in nuclear matter. Also reported 
are the n{k) of nuclear matter within the SCGF approach for the Reid potential. 



quasihole contribution (for those partial waves for which exists) whereas the high momentum 
components are given by the continuum part. In order to show the importance of the continuum 
part of the spectral functions as compared to the quasihole contribution, the particle numbers 
for each partial wave including the degeneracy of the states. 



/ep roc 
du / dkk'^Sij{k,uj), (246) 
-oo Jo 



separating the quasihole and the continuum contributions are reported in Table ^. In this 
approach, there are only 14.025 out of 16 nucleons in the quasi-hole states. Another 1.13 
"nucleons" are found in the 2hlp continuum with partial wave quantum numbers of the s and 
p shells, while an additional 0.79 "nucleons" are obtained from the continuum with orbital 
quantum numbers of the d,f and g shells. The quasi-hole strength can be identified with the 
experimental spectroscopic factor. In this approach, for the pi/2 it turns to be 90% which is 
too large compared with the experimentally determined value, 63%. This discrepancy , can be 
associated to the fact that this approach is treating mainly the short-range correlations. Long- 
range (low-energy) correlations, will typically yield another 10% reduction of the quasihole 
strength [Q, |211| , |212|| . Similar results are obtained with variational Monte Carlo techniques, 
calculating the overlap between a correlated ground state wave function for ^^O and the quasi- 
hole states of the {A — 1) residual system (^^A^) [|213|| . However, an appropriate treatment of 



center of mass motion can enhance the spectroscopic factor by up to 7% [21^] increasing again 
the discrepancy with the experimental data. The absolute value for the spectroscopic factors 
seems not to be understood yet. 



83 



In 


^qh 
ft 




n 


Sl/2 


3.120 


0.624 


3.744 


P3/2 


7 "iT A 
1 .Oi4 


U.OOZ 


( .040 


Pl/2 




u. i ( o 


O. ( 04 


"5/2 


n n 
u.u 


U.Z04I: 


n 9'?/! 


"3/2 


n n 
u.u 


u. iyo 


u. iyo 


/7/2 


0.0 


0.117 


0.117 


/5/2 


n n 

u.u 


n 1 AO 


u . x^u 


99/2 


0.0 


0.040 


0.040 


97/2 


0.0 


0.064 


0.064 


E 


14.025 


1.920 


15.945 



Table 5: Distribution of nucleons in ^^O. Listed are the total occupation number n for 
various partial waves and also the contributions from the quasihole {n'^^) and the continuum 
part (n"^) of the spectral function separately. The last line gives the sum of particle numbers 
for all partial wave listed. 



The sum of the particle numbers listed in Table |^ is slightly smaller (15.945) than the 
particle number corresponding to ^^O. There are several possible sources for this discrepancy: 
First of all the analysis only accounts for momenta below 3.3 fm~^ and the partial waves with 
/ > 4 have not been considered. The restriction in k is determined by the choice of Nmax in 
truncating the "box basis". Inspecting the decrease of the occupation numbers listed in Table 
I with increasing / one can expect that the "missing" nucleons may be found in partial waves 
with / > 4. Furthermore, however, one must also keep in mind that this approach to the single- 
particle Green function is not number conserving, as the Green functions used to evaluate the 
self-energy are not determined in a self-consistent way [|I^ . 

Similarly to Fig.^ for nuclear matter, the momentum distribution including the quasi-hole 
states obtained with various energy cutoffs is shown in Fig. RTI. The quasihole part reflects the 



momentum components that one can detect in knock-out reactions with small energy transfer. 
As in the nuclear matter case, the high momentum components due to short-range correlations 
can be detected only for large excitation energies. 

The total momentum distribution is displayed again in Fig.^ and compared to predictions 
for nuclear matter. In order to enable the comparison with nuclear matter results, the mo- 
mentum distributions have been normalized such that / d^kn{k) yields 1. In order to show 
the sensitivity of the calculated momentum distribution to the nucleon-nucleon interaction, the 
results obtained for Reid soft-core in SCGF are also shown |T89[| . The nuclear matter results 
for the OBE-B potential have been obtained by using a self-energy which contains only second 
order terms on the propagation of holes ||194]| . The OBE-B is consider to be softer than the 
older Reid potential. This is reflected by the fact that the momentum distribution obtained for 
the Reid potential yields larger occupation values for k > kp than those obtained for the OBE 
potential. The comparison between nuclear matter and finite nuclei seems to indicate that the 
enhancement of the momentum distribution predicted in this approach for high momenta is 
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Figure 33: Comparison of different approaches to calculate the momentum distribution in 
^^O. The results obtained from the spectral function are given separately including both the 
quasihole and the continuum part (total). The curve labeled LDA corresponds to the results 
from the local density approximation (LDA) obtained in Ref. ||193|| , while the FHNC belongs 



to the Fermi hypernetted chain approach of Ref. ||109|| and VMC stands for the variational 
Monte Carlo calculation of Ref. ||141|| . 



below the corresponding prediction derived for nuclear matter. These problem has been clar- 
ified in Ref. [r^] were an LDA estimation of the contribution to n{k) of the partial wave not 



included in the calculation brought the tail of n{k) in ^^O in close agreement with the nuclear 
matter results, when calculated with the same interaction. 

In Fig. ^ we compare the n{k) in ^^O obtained in different approaches. The results obtained 



in the Green function approach (total) are compared with those obtained in LDA from Ref. 



193|| , which were obtained using the Urbana V14 interaction. Also shown are the results, 
obtained by a direct FHNC calculation of the expectation value of the number occupation 
operator in a correlated wave function of the ground state of ^^O, where the correlations were 
optimize by minimizing the ground state energy calculated with the semi-realistic Afnan and 
Tang interaction ||109|| and finally the variational Monte Carlo calculations of Ref. [ |141|| obtained 



with a ground state wave function containing two- and three- operatorial correlations which 
minimize the ground state energy for the Argonne V14 interaction. 

Now we want to illustrate the effects of correlations, which are taken into account in the 
Green's function approach beyond the BHF approximation, on the ground state properties of 
^^O . To calculate these properties one needs also the non-diagonal part of the density matrix 
which is given by 

nij{ki,kn) = f du-Im f V(/ci | a)ijgij{a, p;uj){P \ kn)ij] . (247) 
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and contains as in the case of the spectral function, a continuous contribution and a part 
originating from the quasihole states, 



~ qh 



n 



{ki, kn) = E Zrijiki I T);,(T | A;„);,. (248) 



With this density matrix, the expectation value for the square of the radius can be calculated 
according to 

^max 

« I I O = E2(2j + 1) E {k^ I | A;„);n,,(fc„ A;„). (249) 

Ij i,n=l 

The matrix elements of in the basis states are given by 

{ki I I kn)i = NuNni / drr%{kir)ji{knr). (250) 

Jo 

The total energy of the ground state is obtained from the Koltun sum rule ( [76D 



£;o^ = E2(2j + l) E / rfc^;. Ui + ^ + • (251) 



Ij i=l """" " V^"" / \ X 



Using the previous equations we calculate the contributions of the different partial waves 
to the binding energy, the results are reported in Table As a first step we consider the 
Hartree-Fock (HF) approximation. The resulting binding energy per nucleon (-1.93 MeV) is 
quite small. This is probably due to the use of a G matrix calculated at the saturation density 
of nuclear matter, which overestimates the Pauli effects as compared to a BHF calculation 
directly for ^^O. 

The treatment of the Pauli operator is improved by adding the 2plh part ,with the cor- 
responding correction term, to the self-energy. This approximation can be considered as an 
approximation to a full BHF and we will labeled "BHF". This correction increases the bind- 
ing energy to -4.01 MeV. This number is in reasonable agreement with self-consistent BHF 
calculations performed for ^^O using the same interaction ||198|| . However, as the single parti- 
cle states are more bound than the single particle states obtained in the HF calculations, the 
gain in binding energy is accompanied by a reduction of the calculated radius of the nucleon 
distribution. 



The inclusion of the 2hlp contributions to the self-energy reduces the absolute values of the 
quasi-hole energies. Despite of this reduction, the total binding energy is increased compared 
to BHF calculations. This increase of the binding energy is mainly due to the continuum part 
of the spectral function. Comparing the various contributions, one finds that only the 37 % of 
the total energy is due to the quasi-hole part in Eq. ( p51| ). The dominant part (63 %) results 
from the continuum part of the spectral functions although this continuum part only represents 
11 % of the nucleons. The inclusion of 2hlp terms increases also the radius, moving the results 
for the ground state off the Coester band. 

The relative importance of the various contribution to the single-particle strength can also be 
seen from inspecting Fig. ^ which shows the density profile of ^^O and some of its components. 



This figure also demonstrates that the strength located in single-particle states with / > 1, 
which would not be occupied in the mean field approach, provide a small but non-negligible 
contribution. 

Another way of describing correlation effects in a single-particle basis is the representation 
of the single-particle density matrix in terms of natural orbits [pT5| , |216| , |217| , |218 . 
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2.76 
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17.19 
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-103.62 


35.94 


-5.84 


Z > 1 c 














-98.87 


63.17 


-12.27 


E/A 
(r) 


-1.93 
2.59 


-4.01 
2.49 


-5.12 
2.55 



Table 6: Groundstate properties of ^^O. Listed are the energies e and kinetic energies t of the 
quasihole states (qh) and the corresponding mean values for the continuum contribution (c), 
normalized to 1, for the various partial waves. Multiplying the sum: l/2(t + e) of these mean 
values with the corresponding particle numbers of Tab. I, one obtains the contribution AE to 
the energy of the ground state. Results are presented for the Hartree-Fock (HF), Brueckner- 
Hartree-Fock (BHF) and the Green function approach (Total). The particle numbers for the 
qh states in HF and BHF are equal to the degeneracy of the states, all other occupation 
numbers are zero. The results for the radii are given in fm, all other entries in MeV. 
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Figure 34: Density distribution of nucleons in ^^O as a function of the distance from the 
center of the nucleus. 



3.4 Correlations in Nucleon Knock-Out Experiments 

The uncorrelated Hartree-Fock state of nuclear matter is given as a Slater determinant of plane 
waves, in which all states with momenta k smaller than the Fermi momentum kp are occupied, 
while all others are completely unoccupied. Correlations in the wave function beyond the mean 
field approach will lead to occupation of states with k larger than kp. Therefore correlations 
should be reflected in an enhancement of the momentum distribution at high momenta. Indeed, 
as we have already discussed above, microscopic calculations exhibit such an enhancement for 
nuclear matter as well as for finite nuclei |T^, 12^. At first sight one feels encouraged to measure 



this momentum distribution by means of exclusive (e, e'p) reactions at low missing energies, 
such that residual nucleus remains in the ground state or other well defined bound state. 
From the momentum transfer q of the scattered electron and the momentum p of the outgoing 
nucleon one can calculate the momentum of the nucleus before the absorption of the photon 
and therefore obtain direct information on the momentum distribution of the nucleons inside 
the nucleus. 

This idea, however, suffers from a little inaccuracy. In such exclusive (e, e'p) experiments 
one does not measure the whole momentum distribution but rather the spectral function, at 
the energy which correspond to the specific final state. As we have seen in the discussion of 
the preceeding section the spectral function at low energies does not show this enhancement of 
the high momentum components, which shows up only if one integrates the spectral function 
up to high excitation energies (see Fig. ^ and discussion there). A general review on nucleon 
knock-out by means of electromagnetic probes is presented in the book of Boffi et al.[ [^19|| 

A second problem is related to the fact that the nucleon, which is knocked out in such an 
(e, e'p) process, feels the interaction with the remaining nucleus. It will be retarded and might 
also get reabsorbed. These effects of the final state interaction (FSI) can be taken into account 



88 



10' j J 




Figure 35: Reduced cross section for the ^^0{e,e'py^Ngs reaction in parallel kinematics. 
Results from |24| with and without (plane wave impulse approximation, PWIA) inclusion of 
FSI are compared to experimental data taken by Leuschner et al.||220|| 



by means of an optical potential. As an example for the importance of FSI effects we show in 
Fig. ^the reduced cross section for the ^^0{e, e'pY^Ng^ reaction with and without inclusion of 



FSI effects. The reduced cross section is defined as the cross section divided by the elementary 
electron - nucleon cross section times a kinematical factor. Data were taken |220l in so-called 



parallel kinematics. This means that the missing momentum of the proton, the momentum of 
emerging proton minus the momentum transfer is parallel or antiparallel (negative values) to the 
momentum transfer. Note that the missing momentum roughly corresponds to the momentum 
of the nucleon before absorption of the virtual photon. The retardation effects yield a reduction 
of the momentum for the outgoing proton. Therefore FSI effects lead to a shift of the cross 



section to smaller missing momenta, which can nicely be seen from Fig. p5 . 

The data displayed in Fig. ^ exhibit results only up to moderate values of the missing 
momenta. Experiments trying to explore high missing momenta were performed at MAMI in 
Mainz and NIKHEF in Amsterdam | [^22| |. These data can be well reproduced over a range 
of missing momenta going up to 600 MeV/c by calculations which account for the FSI using 



a relativistic model for the optical potential ||223|| . The shape of the calculated cross section 



however, is rather insensitive to the use of spectral functions, which are derived either from 
mean field or more general single-particle Green's function (see Fig. 36). The only difference is 



the global spectroscopic factor which is adjusted to reproduce the total cross section. 

This demonstrates that exclusive one-nucleon knock-out experiments only yield limited 
information on NN correlations. Therefore one tries to investigate exclusive (e, e'NN) reactions, 
i.e. triple coincidence experiments in which the energies of the two outgoing nucleons and the 
energy of the scattered electron guarantee that the rest of the target nucleus remains in the 
ground state or a well defined excited state. The idea is that processes in which the virtual 
photon, produced by the scattered electron, is absorbed by a pair of nucleons should be sensitive 
to the correlations between these two nucleons. 

Unfortunately, however, this process which is represented by the diagram in Fig. |37|a, com- 
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Figure 36: Reduced cross section for the ^^0(e, e'p)^^N reaction leading to the ground state 
(1/2^) of ^^N in the kinematical conditions considered in the experiment at MAMI (Mainz) 
2l[| . Results for the mean-field description (HF) and the correlated spectral function (Full) 



are presented. 



petes with the other processes described by the diagrams of Fig. |37[ b and c. These last two 
diagrams refer to the effects of final-state-interaction (FSI) and contributions of two-body cur- 
rents. Here we denote by final state interaction not just the feature that each of the outgoing 
nucleons feels the remaining nucleus in terms of an optical potential. Here we call FSI the 
effect, that one of the nucleons absorbs the photon, propagates (on or off-shell) and then shares 
the momentum end energy of the photon by interacting with the second nucleon which is also 



knocked out of the target. The processes described in Figs. |37| a and |37[b, correlations and 
FSI, are rather similar, they differ only by the time ordering of NN interaction and photon 
absorption. Therefore it seems evident that one must consider both effects in an equivalent 
way. Nevertheless, most studies up to now have ignored this equivalency but just included the 
correlation effects in terms of a correlated two-body wave function. For the sake of consis- 
tency one should assume the same interaction to be responsible for the correlations and this 
two-body FSI. Correlations can be evaluated in terms of the Brueckner G-matrix while the 
T-matrix derived from the very same interaction should be used to determine FSI. 

The nine-fold differential cross section of such an (e, e'2N) reaction can be written as a 
product of a matrix element of the leptonic current times the matrix elements of the corre- 
sponding hadronic current < > calculated between the initial and final nuclear state 

K\3, <r>\^ (252) 



dEld^]ldE2dfi2d^M^]' 



with K a kinematical factor. This expression can be rewritten into the form [224, 225] 



dV 1 1 ~ ~ r . 

■ ~ , - ,^ ^F/^n/ = 7 El E2 cTMott {vcWl + vtWt + vsWtt + viWlt \ 

X (27r) 5{Ef - Ei) (253) 
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Figure 37: Diagrams for the different processes contributing to the (e, e'2N) reaction. 
Diagram (a) and (b) show the absorption of the photon by a single nucleon. The nucleon- 
nucleon correlations are described by the G matrix. Diagram (c) depicts photon absorption 
via meson exchange (MEC) or isobaric currents (IC) 



where Ei, E2 and pi,p2 denote the energies and momenta of the outgoing nucleons, respectively. 
The virtual photon created in the electron scattering process carries a transfered momentum q 
and energy uj. The leptonic structure functions Vi {i = C, T, S, I) are defined by 



(M. 



A*. 2 



Vt = tan 



2V q 



- If (254) 

Here, 9e is the angle of the scattered electron with respect to the incident electron beam and 
Eg, E'^ are the energies of the incident and the scattered electron, respectively. The nuclear 
structure functions Wi {i = L,T,TT, LT stands for longitudinal, transverse, etc.) contain the 
matrix elements of the nuclear current operator for a given photon polarization A. These matrix 
elements have to be calculated for the various processes displayed in Fig. 

The relative importance of the various contributions under different kinematical setups, 
can be explored in calculations of the structure function Wi for nuclear matter at saturation 
density P26| ]. As a first example we consider the longitudinal structure function for the knockout 



of a proton-proton pair. In this case one can ignore the effects of meson-exchange currents 
(MEC) and isobar currents (IC) displayed in Fig. since two protons do not exchange charged 
mesons. One of the protons is assumed to be emitted parallel to the momentum of the virtual 
photon with an energy of Tp 1 = 156 MeV, while the second is emitted antiparallel to the photon 
momentum with an energy of Tp_2 = 33 MeV (see left part of Fig. [38| ). This is called the 'super- 
parallel kinematic', which should be appropriate for a separation of longitudinal and transverse 
structure functions. In this situation the dominant contribution to the longitudinal response 
function is due to correlation effects (dashed line). But also the FSI effects contribute in a non- 
negligible way to the cross section (curve with circle symbols), although the two protons are 
emitted in opposite directions. The effects of FSI are much more important, if we request that 
the two protons are emitted in a more symmetric way. As an example we show the longitudinal 
structure function for (e, e'pp), requesting that each of the protons carries away an energy of 70 
MeV and is emitted with an angle of 30° or —30° with respect to the momentum transfer q of 



the virtual photon. Corresponding results are displayed in in the right part of Fig. p8|. For this 
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Figure 38: Longitudinal structure function for the knockout of a proton-proton pair as a 
function of the photon momentum g, keeping the photon energy constant at u = 215 MeV. 
The left figure shows results for a 'super parallel' kinematical situation with angles 9p i = 0° 
and 0'p2 = 180° of the two protons with respect to the direction of the photon momentum, 
while the figure on the right-hand side displays results assuming 9p- = ±30°. The dashed 
line represents the contribution from correlations, the line with circles those of the FSI and 
the solid line the coherent sum of these two 



kinematical situation the FSI contribution is much more important than the correlation effect. 

The situation is even more complicated in the case of the (e, e'pn) reaction, since in this 
case also MEC effects need to be considered. Results for the longitudinal structure function for 
(e, e'pn) assuming the same kinematical setup as in the right part of Fig. |38| are displayed in 
Fig. The resulting structure function for (e, e'pn) is almost an order of magnitude larger than 
for the corresponding (e, e'pp) case. The dominating contribution to the longitudinal response 
is again the correlation part. Comparison with Fig. ^demonstrates that the pn correlations are 
significantly larger than those for the pp pairs. This supports our conclusion from discussing the 
results of table ^ that the pionic or tensor correlations which are different for isospin T = and 
T = 1 pairs play an important role and are even more important than the central correlations, 
which are independent of the isospin. Note that the MEC contribution for this 'superparallel 
kinematic' is smaller than the correlation effect. This is due to a strong cancellation between 
the pion seagull and the pion in flight contributions to the MEC. In fact, the dominant MEC 
contribution in this kinematical setup is due to the coupling of the photon to the p meson. This 
does not hold for other situations, in which the pion contributions to MEC are dominant. 

The calculation of (e, e'NN) reactions for finite nuclei typically account for the FSI effects 
only on the level of the mean field approximation. This means one considers an optical potential 
for the outgoing nucleons but ignores the FSI effects which are due to the residual interaction 
between the two ejected nucleons. As a typical example we will consider again the closed shell 
nucleus ^^O as the target nucleus. For an exclusive experiment leading to a discrete final state 
of the daughter nucleus with well defined angular momentum J and isospin T the initial state 
$j , for which the matrix elements < > (see Eq. ( p52| ) of the current operator have to be 
calculated can be expanded in terms of correlated two-hole wave functions 

$r(ri,f2) = <^.2 < n,r2\^2Wii^2JT > (255) 

where we have used the nomenclature of the coupled cluster method introduced in subsec- 
tion 2.2. The expansion coefficients a^j^j determined from a configuration mixing calcu- 



92 



1.5 




250 300 350 400 450 

Photonimpuls [MeV/c] 



Figure 39: Longitudinal structure function for (e, e'pn) for the kinematical condition as 
displayed in the right part of Fig. The contribution of MEC is shown by the solid line 
with box symbol 



lation of the two- hole states in ^^O, which can be coupled to the angular momentum and 
parity of the requested state within the model space considered for the treatment of long range 
correlations ||22 7| . 



As an example, exhibiting the effects of correlations, we display in Fig. ^ the two-body 
densities, resulting from the knock-out of two nucleons from the pi/2 shell in ^^O. In Fig. ^ such 



two-body densities are displayed for a fixed fi = {xi = 0,yi = 0, zi = 2 fm) as a function of r2, 
restricting the presentation to the X2, Z2 half-plane with (x2 > 0, ?/2 = 0). The upper left part of 
this figure displays the two-body density without correlations (S'2 = 0). One observes that the 
two-body density, displayed as a function of the position of the second particle r2 is not affected 
by the position of the first one fi. Actually, the two-body density displayed is equivalent to the 
one-body density. This just reflects the feature of independent particle motion. If correlation 



effects are included, as it is done in the upper right part of Fig. one finds a drastic reduction 
of the two-body density at r2 = fi accompanied by a slight enhancement at medium separation 
between n and r2- 

In order to amplify the effect of correlations, the lower part of Fig. ^ displays the corre- 
sponding correlation densities (i.e. the corresponding amplitudes S'2 squared). While the left 
part shows the correlation density for the removal of a proton-proton pair, the corresponding 
density for a proton-neutron pair is displayed in the right part. Comparing these figures one 
sees that that the pn correlations are significantly stronger than the pp correlations. This is 
mainly due to the presence of pionic or tensor correlations in the case of the pn pair. This part 
of Fig. ^ also exhibits quite nicely the range of the correlations. This range is short compared 
to the size of the nucleus even in the case of the pn correlations. All results displayed in this 
figure have been obtained using the Argonne V14 potential for the NN interaction]^]. 

Results for the cross section of exclusive (e, e'pn) reactions on ^^O leading to the ground 



state of ^^N are displayed in Fig. The calculations have been performed in the super-parallel 
kinematic, which we already introduced before. The kinematical parameters correspond to 
those adopted in a recent ^^0(e,e'pp)^^C experiment at MAMI ||228|| . In order to allow a direct 



comparison of (e, e'pp) with (e, e'pn) experiments, the same setup has been proposed for the 
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Figure 40: Two-body density for the removal of two nucleons from shell in ^^O. The 
upper part displays the results for the total density \E'2 without (left) and with inclusion 
(right) of NN correlations. The lower part shows the defect function squared for pp (left) 
and pn knock-out (right). See text for further description. 



first experimental study of the ^^0(e,e'pn)^^N reaction ||229||. This means that we assume an 



energy of the incoming electron Eq = 855 MeV, electron scattering angle 9 = 18°, uo = 215 
MeV and q = 316 MeV/c. The proton is emitted parallel and the neutron antiparallel to the 
momentum transfer q. 

Separate contributions of the different terms of the nuclear current are shown in the figure 
and compared with the total cross section The contribution of the one-body current, 



entirely due to correlations, is large. It is of the same size as that of the pion seagull current. 
The contribution of the A-current is much smaller at lower values of pb, whereas for values 
of Pb larger than 100 MeV/c it becomes comparable with that of the other components. It is 
worth noting the the total cross section is about an order of magnitude larger than the one 
evaluated for the corresponding {e,e'pp) experiment P30|| . This confirms our finding about the 
relative cross sections for pp and pn knock out, which we have discussed above for the study in 
nuclear matter. 



The right-hand part of Fig. ^ shows the quantities as the left part but calculated with the 
simpler prescription of correlations, i.e. by the product of the pair function of the shell model, 
described for l^*" as a pure (^1/2)"^ hole, and of a Jastrow type central and state independent 
correlation function. The large differences between the cross sections in the two parts of Fig. ^ 
indicate that a refined description of the two-nucleon overlap, involving a careful treatment 
of both aspects related to nuclear structure and NN correlations, is needed to give rehable 
predictions of the size and the shape of the (e, e'pn) cross section. 
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Figure 41: The differential cross section of the ^^0(e,e'pn) reaction as a function of the recoil 
momentum pe for the transition to the 1^ ground state of ^^N (£'2m = 22.96 MeV), in the 
super-parallel kinematics with Eq = 855 MeV, and cu = 215 MeV q — 316 MeV/c. The recoil- 
momentum distribution is obtained changing the kinetic energies of the outgoing nucleons. 
Separate contributions of the one-body, seagull, pion-in-flight and A-current are shown by 
the dotted, short-dashed, dot-dashed and long-dashed lines, respectively. Positive (negative) 
values of pe refer to situations where is parallel (antiparallel) to q. The calculations 
leading to the left part used correlated wave functions calculated in the frameworkof the 
CCM, while for the right part a simple Jastrow correlation function has been considered. 



4 Conclusion 

The main aim of this review has been to demonstrate that nuclear systems are very intriguing 
many-body systems. They are non-trivial systems in the sense that they require the treatment 
of correlations beyond the mean field or Hartrcc-Fock approximation. Therefore, from the point 
of view of many-body theory, they can be compared to other quantum many-body systems like 
liquid He, electron gas, clusters of atoms etc. A huge amount of experimental data is available 
for real nuclei with finite number of particles as well as for the infinite limit of nuclear matter 
or the matter of a neutron star. 

It is a challenge for theoretical physics to develop many-body theories for a description of 
these data, which are based on a realistic model for the NN interaction and yield predictions 
of the many-body data which are free of any parameters. A lot of progress has been made 
during the last years to develop such tools of many-body theory as the Brueckner hole-line 
expansion, the coupled cluster or "exponential S" method, the self-consistent calculation of 
Green's function and variational calculations, which are based on cluster expansion techniques 
as well as Monte Carlo methods. These techniques are very useful for the study of nuclear 
structure but are also used in investigationis of other quantum-many-body sytems. 

The different techniques yield predictions for the bulk properties of nuclear systems at 
normal densities which tend to agree with each other. This is certainly true for the treatment 
of short-range correlations. Additional effort may be required for a reliable treatment of long- 
range correlations in finite nuclei. Those correlations arc related to the low energy configurations 
within a shell-model basis. Shell-model configuration mixing calculations using Monte-Carlo 
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techniques ||231|| or other sophisticated methods to deal with nuclear structure calculations in 
large model spaces |23^ might be appropriate and should be combined with the formalisms to 
handle short-range correlations. 

Many-body theories have reached a level of sophistication that their predictions can be 
considered as a reliable test of models for the NN interaction. Characteristic differences can 
be observed if either a complete meson exchange model is used, which leads to a non-local NN 
interaction, or if the local approximation is considered. In any case, however, one finds that a 
two-nucleon interaction alone does not lead to the empirical saturation point of nuclear matter 
and also fails to reproduce binding energy and density of finite nuclei. The calculations tend 
to predict binding energies which are too small and/or densities which are too large. 

This can be corrected by introducing empirical three-nucleon forces. It is not clear, however, 
whether such three-nucleon forces simulate sub-nucleonic degrees of freedom (A excitations of 
the nucleons) or relativistic features as they are contained in the Walecka model[^]. Further 
studies are needed to clarify this point. One must try to inspect special observables like e.g. the 
spin-orbit splitting in the single-particle potential, the energy dependence of the optical po- 
tential or the response functions for (e, e'p) experiments |P33| , P34|| , which are sensitive to the 
relativistic features. 

There is also a challenge for a cooperation between experimental and theoretical physics to 
search for observables, which test the significance of the short-range NN correlations and thereby 
the short-range structure of the underlying NN interaction. The study of exclusive (e, e'NN) 
experiments seems to be an appropriate tool to explore the proton-proton and proton-neutron 
correlations. 
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